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ABSTRACT 


The theory of optimal contro! presents a powerful method of controller design. 
Ihe basis of the method 1s selection of a “performance index” which compares the 
actual to desired controller performance for such items as output response and energy 
consumption. In most previous work, the design process has radically simplified the 
performance index for analytic feasibility and practical utiiity. [n order to provide a 
more accurate and versatile method, state of the art numerical optimization methods, 
using the Automated Design Synthesis Program, are applied to numerical modeling of 
muitivariable controllers, using the Dynamic Simulation Language. The resulting 
Optimum values stemming from the analysis are the controller gains which minimize 


the desired performance index for a specified set of system constraints. 
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I. INTRODUCTION 


The design of controllers has been an engineering endeavor since the 19205. Тһе 
road to success was significantly changed by the introduction of the International 
Business Machine Model 360 computer in the 1960's, which allowed the control 
designer to develop numerical models for various systems and efficiently test his design 
before prototyping. The development of FORTRAN programs such as he 
Continuous System Modeling Program (CSMP), and Dynamic Simulation Language 
(DSL) greatly simplified the modeling proeess and allowed the designer to have highly 
accurate numerical tools at hand, rather than having to relv on elosed form or 
reduced-matrix type solutions such as the Riccatti equation. 

Along similar lines the engineering discipline of optimization has developed. 
Numerical optimuzation, using mainframe and micro computers, has rapidly developed 
over the past ten years with the work of G. N. Vanderplaats and others. His 
FORTRAN programs Control Program for Engineering Synthesis (COPES) and most 
recently Automated Design Synthesis (ADS) have provided routines that implement all 
of the known methods available to do numerical optimization. 

The goal of optimization has been always to provide the best for the least, and 
this credo ean be applied to the design of all controllers. One method of aceomplishing 
this task will be demonstrated by the author. It involves the modeling of single input, 
single output (SISO) or multiple input, multiple output (MIMO) systems using DSL, 
specifying an applicable performance index by which to judge the system, defining 
constraints which specify the desired system performance characteristics as well as 
verifying system stablity, and then calling upon ADS. Through iterations the 
optimization routine finds the necessary values for system gains (defined as design 
variables) whieh will minimize the performance index (defined as the objective 


function). The final result being the “best” controller design. 


il. MODERN CONTROLLER DESIGN 


A wtNALYSIS TEGHNIQUES 

Two general methods of controller design have developed, each having certain 
advantages, depending upon the application of the controller involved. The first is 
frequency domain analysis, in which the response of the desired system to a sinusoidal 
input is examined. The other method is time domain analysis, where the system 
response to various time dependent inputs such as impulse, step, or ramp functions are 
investigated using a state-space representation of the system. This second method 
lends itself more easily to numerical modeling techniques that can be programmed ona 


computer, and will be the method of choice in this investigation. 


ЕЕ NUMERICAL ANALYSIS TECHNIQUES 
Р 

The CSMP- III program was developed бу International Business Machine 
Corporation (IBM) to provide users with a simplified method for numerically solving a 
system of ordinary differential equations and related mathematical problems. [Ref. Гір. 
2]. Today it is a widely used and respected tool for the design of control systems. 

Some drawbacks to CSMP were found when the initial considerations were 
made on how to perform numerical controller optimization. The most significant areas 
of concern were that CSMP uses the FORTRAN/HX complier, it 1s only available 
locally in the Multiple Virtual System (MVS) environment, it uses single precision 
variables throughout, and that graphics are available only after processing through the 
DISSPLA graphics system. These facts led to a decision to select a different method to 
do controller optimization. 

2,/Фв51.ГУМ 

DSL/VM 1s a recent product from IBM [Ref. 2:p. in]. The program is similar 
to CSMP in that simplified source code is translated into standard FORTRAN 
Statements (as output to the file FORT FORTRAN), compiled and then run in order 
to perform the desired system simulation. Following the successful run, tabled output 
IS directed to the user along with any desired graphic displays. Fig. 2.1 1s a flowchart 


of this procedure. 


Significant improvements, from CSMP to DSL inelude (80 of The 
FORTRAN/VS compiler, Virtual Machine (VM) or MVS environment availability, on 
line graphics capability, double precision calculations throughout, and automatic 
linkage to user defined libraries. Additionally, a one time job initialization segment 1s 
available which is separate from the model initial segment. Thus the programmer can 
easily rerun the time dependent simulation for the the same model, after changing any 
desired parameters, within the same Job specification. These facts were significant 


reasons (or selecting DSL as the modeling tool to use for controller optimization. 
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Figure 2.1 DSL/VS FloweChart 


HT. OPTIMIZATION 


a ee BOTS ТЕСТ СЕХ 
l. Basic Viethcds 

The determination of the “best” value for any design requires that several 
items be considered. One is the resource that can be varied, made up of what are called 
the design variables. Also included are the mathematical reiationships that describe 
how the resource relates to the design environment. These relauoaships define what 
are called the objective function and constraint functions. Many diverse methods have 
been developed to determine the values of the design variables that cause the obiective 
Поп о 06 41 15 оргшайга value. With the advent of the compuier, numerical 
computation methods have become the most practical ones for treatment of COMPEX 
ystems, supplanting grapnical and analyuc determinations. 

2. Numerical Methods 

Many different approaches nave been developed for solving the general 
optimization problem by numerical methods. Numerical Optimization Techniques for 
Engineering Design [Ref. 3] describes most of the common analytic techniques and tne 
corresponding algorithms to perform this task using a computer. A major extension of 
this work was the development of the Control Program for Engineering Synthesis 
(COPES) and Constrained Minimization (CONMIN) FORTRAN routines that 
accomplish the various numerical tasks necessary to successfully determine an optimal 
solution to problems that can be appropriately defined. 

a. COPES and CONMIN Programs 

COPES is a FORTRAN program that was developed as a metnod for 

easily coding the parameters required to define an optimization problem for solution by 
a computer. [ts use requires that a FORTRAN subroutine called ANALIZ de written 
which evaluates the objective function and all constraint functions. Previouslv, the 
FORTRAN subroutines grouped under the name CONMIN were then called to 
conduct the desired optimization analvsis. Recent improvements to the CONMIN 
routines have been generated and the package has been renamed Automaied Design 
Synthesis Program (ADS) [Ref. 4]. To use COPES with ADS, the routine COPESA 
should be called. 


A major drawback to the use of COPESA s the specific formatting and 
subroutine requirements. This hinders interactive optimization, and has resulted in the 
selection of an alternative method for performing controller optimization. 

b. ADS Program 

ADS can be called directly from any FORTRAN program (Ref. 4:pp. 3-31. 
Prior to making the call, all input parameters must be defined, and a looping logic 
must be provided in the calling program that causes the objective and constraint 
functions to be evaluated upon the completion of each optimization iteration, as 
depicted in Fig. 3.1. | 

With this simple logic many engineering problems can benefit from the 
advanced capabilities offered by the ADS program. The specific area that has served 
as a focus for this study is the interfacing of controller design, done with DSL, with 


optimization using ADS. 
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Figure 3.1 Program Logic Calling ADS Directly 
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IV. AUTOMATED DESIGN SYNTHESIS/DYNAMIC SIMULATION 
LANGUAGE 


A. OVERVIEW 

Linking of the Automated Design Synthesis Program and the Dynamic 
Simulation Language has been successfully demonstrated, and the acronym ADSL 
represents the close coupling of these two FORTRAN programs. In the following 


discussion the methodology of how thcse two programs interface will be explained. 


B. PROGRAM SEGMENTS 

Figure 4.1 is a flowchart of how the individual segments of a DSL program are 
utilized with the ADS optimization code. The DSL Language Reference Manual 
(Ref. 2] is the source document on the syntax and logic of the DSL language. Chapter 
3 of the manual explains in detail the function of each phase of DSL program 
processing. For the proposed method of controller design each phase 1s utilized. 

1. Program Initialization (INIT). 

This portion of a DSL program is the first executed segment of code сн 
arrays are dimensioned, non-standard variable types (fixed, complex, etc.) are defined, 
and initial values of variables are assigned. Additionally, DSL parameters which 
specify the output variables, the desired integration method, and the technique for 
plotting can be specified. 

The important fact about this section is that it 1s only executed one time. 
Thus the initial values that ADS requires about the status of the optimization, denoted 


by the variable INFO, and the array dimensioning can be done in this section. 


2. Run Initialization (INITIAL) 

The beginning of each system response simulation starts at this point with a 
call to ADS. On the initial call, ADS is provided with best-guess values for the design 
variables, the upper and lower bounds for each design variable, along with the number 
and types of constraints and their allowable ranges. ADS will return new values for 
the design variables to be used in the controller design. On subsequent calls to ADS 
updated values of the objective and constraint functions will be provided to ADS for 


its determination of the optimal values of the design variables. 
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The remainder of this program section sets various constants at the initial 

values that will be needed later in the dynamic section. | 
3. Derivative Section (DERIVATIVE) 

This section is the heart of the simulation. The inputs to the system are 
defined, the system transfer function is specified, and various outputs, states and the 
objective function are evaluated. Great latitude is possible in each of these 
specifications, and for this reason DSL is an outstanding method for analysis of 
controller design. 

This section also determines the integration time step that the simulation will 
take, based upon two criteria. The first is how much time is remaining before a 
specified data output time is reached, and the second 15 how big an integration step can 
be taken while maintaining the specified relative and absolute errors. The output time 
ПИ Ооа through the use of the DELPRT and DELPLI commands. These 
control parameters indicate the desired time intervals for saving intermediate values of 
specified variables for both plotting and listing purposes. 

4. Dynamic Section (DYNAMIC) 

Once the derivative section finalizes what the size of the next integration 
interval will be, it then calculates the magnitude of the various system states, outputs, 
and objective function at the end of the integration interval. This information is then 
passed onto the dynamic section where it can be evaluated to determine the various 
constraint values and whether or not the simulation is complete (specific information 
on the evaluation of these parameters will be left to subsequent discussions of the 
single input, single output (SISO) example). If the simulation is incomplete, 
computational flow ts 15 directed back to the derivative secton. 

When the completion of the simulation run is determined, by either reaching a 
Seeciied slope Of the Objective function, Or a maximum time limit, the DSL routine 
ENDRUN is called. This procedure completes all outstanding calculations and then 
transfers control to the terminal segment. 

5. Terminal Section 

In this last program segment the ADS status variable INFO is checked to 
determine whether or not the optimization is complete. If it is, the DSL program is 
terminated by calling the routine ENDJOB. If further optimization is required, control 


is transferred back to the initial segment. 
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V. SINGLE INPUT SINGLE OUTPUT PROBLEM ANALYSIS 


Но SECOND ORDER SYSTENI 

In order to obtain a baseline for the ADSL program the well behaved, general, 
second order system was taken for examination. Ogata [Ref. 5] completely describes 
this simple system shown in Fig. 5.1 and its responses to various inputs. For simplicity 
this system will be called Case 1. In the present study, the input is taken to be the unit 
step function, R(s)= 1.0 for time greater than zero, and the magnitude © was set equal | 
080. 









2 





с) 


s2 + s(20@) 


Figure 5.1 Case | Block Diagram 


Various analytic optinuzation methods have also been developed [Ref. S:pp. 296-306] 
for this system that can be used to verify the ADSL program. 
Objectives of this baseline investigation were as follows: 
e Determine the most efficient and correct integration method. 
о Determine the most efficient optimization strategy. 


° Эше the variations to be expected from using different performance 
indices. 


e e the most overall-efficient method to obtain optimum design variable 
values. 


A sample program to accomplish these tasks 1s listed in Appendix C. 
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B. PROBLEM FORMULATION 
1. Objective Function 
The fundamental part of any optimization problem is definition of an 
objective function. For this problem the performance index 1s the logical relationship 
that detines how “good” a system is. Ogata [Ref. 5:p. 296], defines the various 
performance indices that have been used in this investigation, and these will be 
described in the following discussion. 
2. Design variable 
The design variable selected for this investigation was the damping ratio ( GB 
[t is a reasonable parameter because it directly effects how the system responds to any 
input over time. Thus if Č is increased, the system output response will be slower and a 
larger accumulated error between the actual to desired output may be generated with 
respect to time. This larger error could then increase the value of the performance 
index/objecuve function. On the other hand, insufficient damping will lead to excessive 


output oscillations which aiso tend to increase the objective function. 


C. APPLICATION OF CONSTRAINTS : 

Ogata [Ref. 5:p. 232] defines transient response specifications for the typus 
second order system. In almost all cases, these standard definitions can be directly 
applied to the response of a DSL model. A brief review of these specifications and 
how they are applied are discussed below. 

1, Maximum Overshoot 

How much the output response of an under-damped system goes beyond that 
which ts desired, is the maximum overshoot. The larger this value 15, (һе more energy 
1s Wasted, as a restoring force must be applied to return the system to the desired state. 
In addition, physical systems may be limited in the amount of amplitude overshoot 
that can be tolerated. A simple constraint of the maximum allowable amplitude, as a 
fraction of the intended steadv state, can be defined. 

2. Delay Time 

The time the system requires to reach fifty percent of its Intended resp omms 
the system delay time. It indicates how fast the system begins to respond, and a 


constraint of a maximum time can be established. 
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3. Rise Time 
When the system initially reaches 100 percent of the desired response, the rise 
time has been reached. Like the delay time, this parameter is also an indication of how 
quickly the system will respond. For under-damped svstems, a constraint of maximum 
rise time can be programmed. 
4. Peak Time 
The time of maximum response 1s called the peak time. It also shows how 
fast the system responds. In cases where a peak time occurs, a maximum peak time 
constraint can be specified. 
5. Settling Time/Completion Time 
Settling time is defined [Ref. S:p. 236] as when the response falls and stays 
within a prescribed percentage (typically two or five percent) of the total response. For 
computational considerations, the settling time was redefined for the two possible 
response conditions: 
1)  Underdamped case: that time when the absolute maximum of a response 
oscillation peak 1s within one percent of the total response transient. 
2)  Overdamped case: that time, greater than five time constants, at which the 
absolute system response is within one percent of the total response. A tune 
constant is defined as that time required for the system output response to 
reach sixty three percent of the final response value. 
In order to keep the above definition separate from settling time, the new consideration 
will be called completion time. 

The different constraint criteria are illustrated in Fig. 5.2, with the constraints 
of maximum amplitude, rise time, delay time, peak time, and completion time, being 


designated by (1) to (5) repectively. 


D. INTEGRATION CONSIDERATIONS 
DSL provides nine optional integration methods which can be grouped into the 
following types: 
о Fixed-step 
e Variable-step 
е Variable-step, Variable-order. 
The default method is the Runge-Kutta Fifth Order variable-step method. Selection of 


the method which is best suited to the controller design problem is discussed below. 
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Figure 5.2 Case 1 Demonstration of Constraints 
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1. Limits of Integration Determination 

For this analysis the objective function has been defined as the system 
performance index. The particular performance indices to be examined are weighted 
time integral errors between the actual and desired output response of the system from 
the time the input is applied to infinity. Noting that infinity is an impractical number 
for a computer (or anything else, for that matter), a reasonable approximation must be 
found. Two criteria were considered for determining what cut-off time to use to mark 
the end of calculations, or time equivalent to infinity. 

One cut-off time criteria would be to simply select the completion time, but 
this has a significant drawback. Different values of the design variabie will cause 
different cut-off times, and thus this method does not provide a consistent manner for 
evaluaung tne obiective function for different designs. A second method is to evaluate 
the rate of growth of the objective function with time. For a stable system the 
BEIeenye funcauon initially increases rapidly due to the mis-match of input and output. 
йе рговге55с5 this error approaches Zero and the rate of change of the objective 
junction slows down. Once the completion time is reached, if the rate of change of the 
objective function is evaluated to be essentially zero, a consistent method of 
determining the cut-off time can be established for all values of the design variable. 

2. Integration Method Analysis 

Typical output and objective functions are shown in Figs. 5.2 and 5.3. These 
rapidly increase from zero and approach a constant value asymptotically. This type of 
behavior is described by Speckhard [Ref. l:p. 268], and is characteristic of “stiff” 
systems. In considering which integration method would best deal with this stiff 
system problem, fixed-step methods were ruled out due to their inability to accurately 
deal with rapidly changing functions. The Runge-Kutta Fifth Order (RK5), Adams 
variable-order (ADAMS), and Gear Full Jacobian (STIFF) variable-step methods were 
compared to determine which would result in the maximum computational efficiency. 

The computational efficiency in this problem can be considered in two parts: 

1) The number of individual integration steps 
2) The number of objective function evaluations requested by the optimization 
routines for gradient determination. 
Ihe conclusion 1s then to use an integration method that maximizes the step size and 
an optimization routine that minimizes the number of function evaluations. The DSL 


Reference Manual [Ref. 2:p. 3.16] describes the reasons for differences in execution time 
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between the Integration methods, and points out that for oscillatory solutions, such as 
those seen for the second order underdamped system response, all three of the methods 
are essentially equal. Consequently, the best method to use 1s the one that generates 
the least number of objective function evaluation requests by the optimizer. 

Three ADSL programs were generated, each using a different integration 
method, but using the same input, system, optimization method, and objective function 
(performance index). Discussion of optimization methods and performance indices will 
follow. In the problem at hand the optimization method was to use an exterior 
penalty function strategy with the Broydon - Fletcher - Golfarb - Shanno optimizer 
and polynomial interpolation for the one dimensional search. For the performance 
index the integral square error (ISE) was employed. The results are compared in Table 
Ihe applied constraints, designated by (1) to (5) repectively in Fig. 5.2 and Table 1, 
Were: | 

(1) maximum amplitude 

ШІ xise time 

(3) delay time 

(4) peak time 

(5) completion time 

Interesting points from this Investigation Include: 

e All three integration methods found the same result within the default tolerance 
of ADS. 

° Each method, starting with the same initial design variable value, calculated a 
different initial objective function value. 

ө The number of objective function calls for each method was almost the same, 


ПБ iene ioc аз опе тпеһе паус predicted, required the fewest. 


E. OPTIMIZATION METHODS 

In order to use ADS, an optimizer and one-dimensional search must be selected. 
Optionally, a strategy may be invoked. 

1, Optimization Strategy 

An optimization strategy is simply the algorithm by which the overall 

optimization problem is transformed into one in which constraints are weighted or 
eliminated and the objective function may be redefined by use of a “pseudo-objective” 
function. Although use of a strategy is not required, some problems have shown 


significant computational improvement as a consequence. 
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TABBE I 
| | CASE 1 INTEGRA TION METHODERESUELTs 


Method RK5 ADAMS ШЕЕ 
Initial Design 1.0 1.0 1.0 
| Initial Objective 4008.9 12334. 5073.9 
| Final Design DET 0.718 0.720 
Final Objective 1.968 IORS 1.965 
Constraint 1 -.361 -.361 -.362 
| Constraint 2 -.157 -.190 -.157 
| Constraint 3 -.12| -.103 -.077 
| Constraint 4 170 ИЮ -.016 | 
| Constraint 3 -2.24 -2.55 -2.69 | 
Number of Runs 49 46 42 | 





a. Sequential Unconstrained Minimization Techniques (SUM T) 


SUMT involves the creation of a pseudo-objective function that is the sum 
of the original objective function and a penalty function. This penalty function may be 
one of four types: 

1) | Exterior Penalty 

2). -Linear Extended Inte noi 

3) Quadratic Extended Interior 

4) Cubic Extended Interior 
All these methods cause the objective function to have its actual value when the design 
variable is within the feasible region, but when the design variable is outside the 
feasible region, the penalty function significantly increases the objective function value. 
The penalty functions listed above vary in the value that will be added to the objective 
and how the transition between the feasible and infeasible regions 15 conducted 
[Re SEE 

5. Augmented Lagrange Multiplier (ALM) 
In this method a pseudo-objective function 1s generated based upon the 


exterior penalty function, but with an additional term in order to create the 
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Lagrangian. It can be shown that the solution to the Lagrangian 1s the solution of the 
original optimization problem [Ref. 3:p. 141]. 
c. Sequential Linear Programming (SLP) 

The basis of this strategy is the simple linearization of the objective 
function and constraints. Rectangular move limits are established and then linear 
programming techniques are applied to optimize the objective function [Ref. 3:p. 155]. 

d. Method of Centers 

Tne SLP method can produce infeasible designs with each iteration. To 
correct this problem the Method of Centers generates an inscribed hypersphere based 
upon the linearized objective function and constraints, then uses linear programming 
techniques to determine the optimum. 

e. Sequential Programming 
Sequential Pregrammung can be of either of two types: 
1) Sequential Quadratic 
ee sequential Convex 
Both of -hese :nethods refine the method for determining the search direction for the 
objective function minimum. In the quadratic method the objective function 15 
transformed into a quadratic function, while the convex method generates the 
Mciproeal of the  objeetive function. These methods create conservative 
approximations to the actual optimization problem which can then be easily evaluated 
(Ref. 4:p. 2]. 
2. Optimizers 
The method by which new design variables are chosen is a result of the 
optimizer selected. Each of the possible methods predicts a new set of design variables 
which should be closer to the optimum design, based upon evaluations of the objective 
function at previous values of the design variables. 
a. Conjugate Direction Method 

This optimizer, developed by Fletcher and Reeves [Ref. 6:pp. 149-154], 1s 
based upon finding the negative gradient of the objective function from some initial 
point and then performing a one-dimensional search in that direction for the minimum 
value of the objective function. This process is repeated from the newly found point, 
and then a new direction is determined based on the previous two points and 


directions. 
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TASE? 
ADS OPTIOT > 


Sirateazy (15 207 

0 - Мопе 

і ЭШИ 7 Хо епПОог ела I! muon 

2 - S'IMT, Linear Extended Interior 

3 - SUMT, Quadratic Extended Interior 

4 - Cubic Extended Interior 

5 - Augmented Lagrange Multiplier Method 
6 - Sequential Linear Programming 

T- Method of Gentens 

8 - Sequential Quadratic Programming 


9 - Sequential Convex Programming 


ОИЕ ТОРТ) 

үн ессе есте 

2 - Davidon-Fletcher-Powell (DFP) 

3 - Broydon-Fletcher-Golfarb-Shanno (BFGS) 
4 - Method of Feasible Directions 

5 - Modified Method of Feasible Directions 


One dimensional Search (IONED) 

| - Golden Section Method 

2 - Golden Section and Polynomial 

3 - Polynomial Interpolation, bounded 

4 - Polynomial Extrapolation 

5 - Golden Section Method 

6 - Golden Section and Polynomial 

7 - Polynomial Interpolation, bounded 

8 - Polynomial Extrapolation | 

Note: One dimensional searches | through 4 are for use with 


Optimizers | through 3, and searches 5 through 8 are for 
use with Optimizers 4 and 5. 
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Description Manual [Ref. 4: pp.3-5]. In ADS, each possible combination is identified 
ру use of a three digit code (XYZ) with the X digit representing the optimization 
strategv, the Y digit the optimizer, and the Z digit the one-dimensional search 
employed. 

Vanderplaats [Ref. 4] points out that although there is a multitude of different 
combinations, eight combinations have, as a result of ENS SEN eno One 
recommended for general applications. These eight combinations were tested using 


otherwise identical ADSL programs. The results are presented in Table 3. 


ПОВІВ З 
ен eS STRATEGY RESULTS | 


Strategy 047 057 857 957 553 235 155 656 
Objective 1.966 1.966 1.985 1.959 1.964 1.9646 1.964 1.972 


Design Variable 0.720 0.719 0.7023 0.729 0.722 0.722 0.722 0.714 


Constraint 1 -.362 -.161 -,155 -,165 -.362 -.162 -.162 -.559 
Constraint 2 т.А45/ -,157 29250 2.054  -.158  -.158 -.158 -.154 
Constraint 3 тоог 5 гэ 7 2-0208 -.057 -.069 -.069 eco 
Constraint 4 -.009 0.000 -.036 +.048 -.018 -.014 -.014 +.033 
Constraint 5 -2.69 -2.84 -2.95 -2.49 -2.6% -2.61 -2.61 -2.89 
Number of Runs 10 27 11 10 G7 51 51 29 


Observations that can be made are: 

• А1 the combinations found essentially the same optimum value. 

ә [here are significant differences in the efficiency of the various methods. If the 
criteria of minimum number of objective function evaluations is used, the 047 
and 957 strategies are the best. 

9 All the combinations found solutions that met all the constraints within the 
ADS default tolerance band. 

9 More objective function evaluations did not necessarily lead to reduction in 
objective function value. In fact, one of the quickest methods (957), found the 


smallest value. 


F. PERFORMANCE INDEX SELECTION 
Ogata [Ref. 5: pp. 296-301]. presents four performance indices that are based 

upon the time-integrated weighted error between the desired output and the actual 
svstem response. These indices are: 

* Integral square-error (ISE) 

*  Integral-of-time-multiplied square-error (ITSE) 

e Integral absolute-error (IAE) 

о  [ntegral-of-time-multipled absolute-error (ITAE) 
These performance indices vary in complexity for analytical evaluation over the 
integration period of zero to infinity. For two of the cases, namely [SE and ITSE, 
exact analytical answers can be found as indicated in Table 4. Bv using ADSL the 


analytic problems are greatly simplified for the other indices. 


TABLE 4 | 
CASE 1 PERFORMANCE INDEX RESULTS 


Index ТАЕ ISE ТЕЛДЕ ЖЕЛЕ 
Strategy 047 047 047 047 
ADSL Objective 1.605 1.014 . 1.966 OT 
Вх cce n.a. 1.000 n.a. ОЛ 
Design Variable ( Ü ) .6605 .5909 7202 5928 
Ес | n.a. .5000 n.a. .5946 
Constraint l -.337 -.300 -.362 -.301 
Constraint 2 22] -.152 -. 157 -.207 
Corsirairt з -.447 -.858 „ШШ ] ЗІ 
Constraint 4 2299 -.717 -.009 -.646 
Constraint: 5 -3.32 -3.64 -2.69 -3.94 
Number of Runs 20 40 10 45 


The four performance indices were placed into identical ADSL programs as the 
objective function for optimization, and the results are shown in Table 4. Significant 


observations that can be drawn from this table are: | 
° As expected, each of the methods found a different optimum value for the 


design variable. 
e All of the methods found an optimum value that did not violate any constraint. 


5 


e Within numerical tolerance, the results are the same as those found by Ogata 
722171 

° The apparent discrepancy in optimum damping ratio for the ISE performance 
index is due to the flat minimum of the objective function in this case (a 


sA 
SET 


Gticremce Of 20% In | gives a performance index variation of less than 29^). 

[t should be noted that the magnitude of the resulting objective functions cannot 
be compared directly to determine which method is "better, This would be like 
comparing apples and oranges, due to the significant mathematical difference in the 


objective functions/performance indices. 


VI. INTEGRAL FEEDBACK CONTROL DESIG 


А. 5ТКТЕ5РХСЕАМРТ 5 

ine time and frequency domain methods discussed in the previous chapter are 
most useful in the analysis of single input, single output (SISO) systems. State-space 
methods give added flexibility in controller design and allow the treatment of complex, 
multiple input, multiple output (MIMO) systems. The concept of state, which 1s the 
basis of modern control theory, embodies the following kev elements: 

e State Variables: the smallest set of n variables, with their initial quantities 
known, that completely describes the behavior of a dynamic svstem for a given 
input over time. [t should be noted that these states do not necessarily have to 
be phvsicallv accessible. 

з State: the condition of a dynamic system at a given time (t) that 1s uniquely 
determined by the initial state of the system and the inputs to the system after 
the initial time (t,). 

, State Vector: the vector made up of the n state variables needed to describe the 
system. The state vector 1s usually designated x(t). 

е Input Vector: the vector made up of the m variables that are inputs to the 
system. The input vector is usually designated u(t). 

ө State Space: is the n dimensional space made. up of the coordinate єє. 
corresponding the individual states. 

Kuo [Ref. 7] describes two types of basic systems that are typical of modern 
control systems: 

D State feedback control 

2) State feedback with integral control 

These systems have been modeled using the combined Dynamic Simulation Language 
(DSL) and Automated Design Synthesis (ADS) optimization code ADSL. 


B. CASE 2 STATE FEEDBACK WITH INTEGRAL CONTROL 

Case 2 1s the next problem to be investigated. In this system control is attained 
by feedback of the state variables through constant gains and an integrator. A signal 
flow graph representative of state feedback with integral control is shown in Fig. 6.1. 


As a matter of interest, if the integrator is removed the resulting system is called a state 
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feedback controller. When the reference for any state feedback controller is set to Zero, 
the resulting system is classified as a regulator. It can further be shown that the 
common proportional-integral-derivative (PID) and rate-feedback control systems are 
special cases of state feedback control [Ref. 7:p. 522]. 

Case 2 systems are of interest because the state feedback systems are generally 
only useful as regulators that have no noise input. Integral control suppresses the 
fluctuations caused by undesirable noise and thus provides for a much more effective 
system. 

1. Model Description 

The system selected for optimization was a dc electric motor speed control 
svstem proposed by Kuo [Ref. 7:pp. 532-536]. In this example the motor shaft angular 
velocity and armature current were selected as the state variables x, and x» 
respectively. Two inputs were considered: 

° wy, asa direct, torque loading acting on the motor shaft 

° wy as the desired constant speed set point value. 
The control and output variables can respectively be represented by the following 
relationships where gj, g5, and g3 are the feedback gains, and the state variables are x4 
and X»: 

2-1 11-25" 53| c(t)dt 

Фи e(t)nc Wess X1. 
The corresponding signal flow graph showing the various relations in the control 
system is shown in Fig. 6.1. The initial conditions and values assumed for the design 
are shown in Tables 5 through 8. The resulting state responses are displayed in Fig. 
бий» 

The goal of the optimization is to determine the values of the three gains -- 
81» £9, and g4 -- necessary to minimize a stated objective function without violating 
specified constraints for controller performance and stability. 

2. ADSL Program 

Appendix D contains the source code listing of the ADSL program which 
models the state feedback system with integral control for Case 2. It follows the same 
form of construction as used for Case І. The initialization of parameters was expanded 
into three distinct groupings: 

1) Array dimensioning and parameter initialization for ADS 


2) DSL control parameter setting 
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STATE FEEDBACK WITH INTEGRAL CONTROL 


Figure 6.1 Case 2 Signal Flow Graph 


3) Model description, constraint evaluation, and Integration completion control 
parameter implementation. | 

Using the results from Case 1, selection of the stiff integration method was made. An 

additional change was to reduce both the relative and absolute error criteria for the 

integration routine from 1.0E-05 to 1.0E-04. This was done In order to increase 

computational efficiency while still providing for an order of magnitude greater 

precision than that expected for the objective function optimization done by the ADS 
program, using the default parameter settings. 
3. Performance Index/Objective Function 

For optimizing, the classical quadratic performance index (QPI) was adopted. 

This performance index: J = о Ох dr u  Ru)dt allows the user to specify 

individual weights for the Q and R matrices. Thus the desired relative importance 

between the error seen in each of the actual state values to the desired state values, and 


for the energy use associated with the control variable over time can be specified. For 


Er 





TABLE 5 
CASE 2 INITIAL CONDITIONS 


Motor Shaft Angular Velocity (x, ) = 0.0 rad/s 
Armature Current (xo ) = 0.0 Amps 

| Feedback Control Signal (u(t)) = 0.0 Volts 

x Input Disturbance Torque (w, ) ^ 0.0 N-m 
| 


Input Reference Set Point Speed (уә ) = 1.0 rad/s 


TABLE 6 
CSE 2 INIMMAL POLE SELECTION 


Initial Pole Selection 
-10.0+ 10.0; -10.0-10.0; -300.0 


Exact Resulting Initial Gains 
81 = -0.383 82 - 0.60 83 = -6.00 


Case 2 each state was selected to be independent of the other, and thus Q and Ë were 
identity matrices. Table 8 shows each of the matrices used in the Case 2 performance 
index. 
4. Exact Initial Gain Determination 
ADS requires that an initial guess for the design variables be provided in order 


to start the optimization. In conjunction with the the determination of these values a 


TABLE 7 
CASE 2 COP СЕТЬ 


CSP = Maximum Value of the Output c(t) = 140% of initial value 


TSP = Maximum Completion Time = 1.5 sec 
USP = Maximum Value of the Control u(t) = 10.0 V 





| TABLE S 
| CASE 2 PERFORMANCE INDEX MATRICES | 


0.9 R = 1.0 
Ev 


Hu cg l Х 1:82Х2-8: | (541 





reasonable assumption has been made that the controller designer be able to specify 
stable poles for the system. Kuo [Ref. 7:pp. 527-528] outlines an analytic, general 
method for the determination of initial state feedback gains without transformation to 
the phase variable canonical form of the system characteristic equation. The method 
requires defining the following terms: 


e Open loop characteristic equation: 


шоо: 4 


e Closed loop characteristic equation with state feedback: 


ДЕ ет Ав 


e Relation Matrix: 


k(s) = [Adj(sI-A)]B 


S 


STATE VALUES 


Figure 6.2 





0.2 0.4 
TIME IN ЗЕС 


Case 2 Initial Condition State Responses 
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Where the indicated matrices are as follows for n states: 

A- the coefficient matrix for the state variables (n x n). 

B- the coefficient matrix for the control input (n x 1). 

G- the siate feedback eain matrix (1 x n). 

I- the identity mairix. 

The relationship СК (5) 7 A((s) - A (s) can then be solved for G using the 
specified initial design system poles as the eigenvalues for the closed loop characteristic 
equation. In the Case 2 ADSL program, shown in Appendix D, the necessary matmx 
operations have been mult:plied-out and the like terms relating the various powers of s 
have been grouped together to yield a system of three equations and three unknowns. 
Matrix solver routines from the LINPAC library (Ref. $8}, DGEFA (401016 precision 


eneral matrix conditioning) and DGESL (double precision general matrix soiver). were 


(JQ 


them called to respe ман decompose and solve the system of equations. and thus 
generate the exact initial state feedback gains. 

. Approxiniate Initial Gain Determination 

The matrix method iust described above works well, and is an efficient method 
for systems of order three or less. With higher order systems the determination of 
adjoint matrices is difficult and ume consuming. To counter this problem, an 
alternative method was devised. The determination of the gains was considered as an 
optimizauon problem, where the constraints and objective function were formulated 
such that the there would be a minimum difference between the poles corresponding to 
the gains determined by optimization using ADS, and the original poles specified by 
the designer. Two different forms of the optimization problem were developed. 

a. Determinant Method 

For the first method, the complex matrix: ( s] - A ^ BG ) was formeatey 

substituting in one of the desired poles for s and setting the initial values for the gains 
to 1.0. Tire determinant of this matrix was then found using the LINPACK routines 
ZGEFA (complex general matrix conditioning) and ZGEDI (complex general matrix 
derterminant). If the correct value of the gain, corresponding to the desired poie was 
substituted into the matrix equation, the resulting determinant would be equal to zero. 
Thus an equality constraint can be established equal to the absolute vaiue of tne 
resulting determinant. This process was then repeated for each of the other desired 
poles. The objective function was then defined as the sum of all the constraint values. 


For an exact solution the objective function wouid be equal to the value of zero. ADS 
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was then called to vary the gain values to find the minimum value of the objective 
function. The final determined gains thus should represent an approximate match to 
the desired initial poles. 

а авт Wasmetremmpted, the results were not as expected. 
Өс 2 эоттеһетезен5 fer variousunitialpoles using the Case 2 model. Looking at 
Dossibie reasons for these outcomes it can be found that the optimization problem, as 
formulated above, was very ill conditioned. When the exact answers were substituted 
as initial guess for the gains, the determinant values for each gain resulted a machine 
zero of {.0E-08. A very small deviation of 0.01 from the exact gains was then tried and 
resulted in a determinant value of greater than 1.0E+06, as did the determinant for the 
initial guess gain values of 1.0. These radically different determinant values produced 
very large gradients for ADS to work with, and in turn, ADS could not find a solution 
that closelv approximated the desired initial poles. 

à. E:igenvalue Method 

in order to overcome the il! conditioning problem, a more robust method 
was devised. ihe real matmx: (A - BG) was formed, using initial guess vaiues for the 
gains of -1.0. The eigenvalues of this matrix are the corresponding poles of the system, 
and were then found using the Eigensystem Routines (EISPACK) subroutine RG (real 
general matrix solver) [Ref. 9]. The resulting largest and smallest, complex eigenvalues 
Were then compared to the largest and smallest, desired complex poles using the 
FORTRAN min and max functions for both the real and imaginary parts. Four 
equality constraints were then defined as the absolute differences between the 
corresponding min or max, and real or imaginary parts of the desired poles and 
calculated eigenvalues. The objective function was then set equal to the sum of the 
constraint values. ADS was then called to minimize the objective, and thus find the 
gains corresponding to the desired poles. Table 10 shows the results for this method 
using the same initial conditions as in the determinant method above. It can be seen 
that this method finds reasonably accurate values for the gains, compared to the exact 
method. 

It should be pointed out that in order to prevent the the optimizer from 
finding a local minimum corresponding to a positive real pole soluticn the following 
mechanisms were programmed into the algorithm: 


өе Strategy 133 was selected. 
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TABLED 
| ‘CASE 2 INITIAL GAIN APPROXIMATION USING DETE RACINAS 


Desired Pole Derived Pole Derived Gains Exact Gains 


| 
| -10.0 + 10.0) -14.33 + 9.66j -125 -.380 
-10,0-10.0} -14:33-9,66) Ж +600 
| 1300.0 -294.8 -8.81 -6.00 
| 
| -10.04-10.0j + 257 4 414 -.920 
K 4332 5151 450 
| 300 -30:1 З 7367 690 
| 210.0 + 10.0; -24.13 + 39] + 5.02 
22 1106-1020) 5529 2200 +44] 
243009. 20001 НӨ -60.0 
| 
100.9 + 10.0} -99.00 + 167 + 1.50 
^ :100.0-10.0j :274.7 + 863 + 6.0! 
300.( + 1.08 4- 2.94 -303.0 
-10.0 + 100.0} E 4.362 + 610 
10:0-100.0) -10:59 + 343 + 600 
-300.0 -201.4 14.13 -303.0 


ADS Strategy 957 


° The default ADS scaling of the optimization problem was disabled (to be 
discussed in a latter section). 

* A check for positive real parts of the smallest cigenvalue was done. If a 

positive eigenvalue was found, then the sign of the first gain value was reversed. 

An observation that should be pointed out, as a consequence of these 

determinations, is how sensitive optimization is to problem definition. Mathematically, 

both of the methods discussed above are equally good for establishing the gains 

corresponding to the desired poles, but the second method results in always dealing 

with numcrical values of essentially the same order. of magnitude, and a much better 

conditioned problem for ADS to solve. The source code for both of the approximate 


gain determination methods has been included in Appendix E. 
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TABLE 10 
INITIAL GALN DETERMINATION USING EIGENVALUES 


Desired Pole Derived Pole Derived Gains Exact Gains 
-10.0 + 10.0] -9.99-- 10.05) -.384 | -.380 

-10.0- 10.0] -9,99-10. 05) — 9909 + .600 
-300.0 -300.0 -5.98 -6.00 

-10.0--10.0) -7.84 + 9.48] -.939 -.920 

-10.0-10.0) -7,84-9,48) 32773 2150 

zou 2212 52272 -.600 

ОКЕ 00] -118.1 + 37.6 1 5.02 

-10.0- 10.0 EO. + 14.6 -- 14,1 

25000-0 -3000.0 2010 -60.0 x 
-100.0 + 10.0} s 75 ТЭ + 1.50 | 
-100.0-10.0) -05.3-8.76)| + 3.25 + 6.01 | 
-200,0 -292.0 212211 2303.0 | 
-10.0 4- 100.0j -10.0+ 99.7} г 2005 --.610 | 
-10.0-100.0) -10.0-99.7| + 598 + .600 
220070 -299.9 250870 250570 


ADS Strategy 133 


6. Constraint Definitions | 

The maximum overshoot and completion time constraints, as previously 
described, were applied directly to this problem. Using a manner similar to the 
maximum overshoot, a new constraint was defined for the maximum allowable value of 
the feedback control signal u(t). 

An additional feature seen in the state feedback problem is that as ADS 
begins its optimization it has no way of knowing whether or not the design variable 
values proposed for the feedback gains will cause the system to become unstable. To 
account for this, two options are available: 

1) Ina reverse manner to the initial gain computations the ADS proposed gains 
could be used to compute the new corresponding system poles. These poles 


could then be evaluated for stability. 
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2) An extreme upper limit could be set for the allowable system output in order 
to monitor for stability. When this condition was reached the run could be 
terminated and an unstable solution equality constraint could be invoked. 

This second method was selected as more numerically efficient as it can be assumed 
that most proposed gains would be reasonably stable, and the cost of pole evaluation 
for every set of ADS proposed gains is high due to the need for finding eigenvalues to 
the closed loop characteristic equation. 

A stability constraint was proposed, in private discussions with Vanderplaats, 
as a way to eliminate unstable solutions. A simple way to look at this method is to 
consider the total design space. The constraints related to specific performance (such 
as completion time) define the borders of the space, and the resulting interior region 1s 
the domain of feasible solutions. Addition of the stability constraint essentially 
generates a “hole” in the feasible region by eliminating possible solutions at specific 
points, while allowing solutions at points close by. Normally the stability constraint 
would be set to a constant, negative value of say minus one, and thus it would never 
define an active boundary of the design space. However, when an invalid design 15 
proposed by ADS, through a request for evaluation of controller gains that result in 
violating the maximum allowable output (an unstable condition), the stability 
constraint would be set to a positive one, and the previous value of the objective 
function would be maintained. When ADS sees a step change in the unstable solution 
constraint, it will adjust the design variables to return this constraint negative. The 
result will be that only valid values of the design variables are considered for 
evaluation. 

Nye and Tits [Ref. 10:p. 1695] proposed the definition of hard and soft 
constraints. All the previously defined constraints for system performance (maximum 
overshoot, completion time, etc.) fall into the category of soft constraints, because they 
are set by the designer to meet flexible goals through the optimization. One type of 
hard constraint they consider is that of negative response of the input. This constraint 
simply checks to see that the response curve takes off in the same direction that the 
Input goes in, otherwise the solution will be divergent. This constraint was 
implemented in a similar manner to the stability constraint check above. 

7. Integration Termination | 
With the introduction of integral feedback and multiple inputs, the 


performance index/objective function can become a ramp, rather than having a plateau 
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as in Case 1. To account for this effect, an additional check for a constant second 
derivative of the objective function was added to the model as a possible criterion for 
the simulation integration termination. 
9. Final Pole Determinations 
Ones ADS completes the optimization, the final system gains are known. In 
order to give the designer information on how much the optimization caused the 
systems poles to be changed from the originally specified ones, the program calculates 
the final system poles. This requires a two step process: 
i) Calculation of the the final closed loop system matrix equation: A- BG 
2) Determination of the resulting matrix eigenvalues by the same procedure 
described in the approximate initial pole to gain determination using 


eigenvalues. 


С. SCALING 

in work bv Cnow [Ret. 11] and Gordon [Ref. 12], where ADS was used in the 
anaiysis ot poie placement techniques, scaling was "turned off" by manipulation of the 
ADS integer work matrix (IWK). The technique to perform this operation is described 
in the ADS Manual [Ref. 4:pp. 12-16]. One reason for not scaling is to increase 
computational efficiency, as the scaling calculations are not required and thus directly 
reduce the time required to run the problem. It was felt that a determination should be 
made as to how scaling effects the optimization problem when using ADSL. 

l. Scaling the General Optimization Problem 

Scaling of any numerical problem is generally considered good engineering 

practice. Vanderplaats (Ref. 3:pp. 136-137] describes in detail the theory of scaling the 
optimization problem by two different methods. First, the initial problem can be 
scaled to provide relatively the same order of magnitude to the design variables, 
constraints, and objective function. In Case 2 this was.done in the initial problem 
шэг Піс Second Scaling is then done to the gradients of the individual 
parameters, so that they also have the same order of magnitude. Scaling of the 
gradients, in addition to the values, is critical because the basic optimization problem 
rests on finding a minimum, or in other words, finding where the gradient goes to zero. 

2. Scaling in ADS 

Scaling within ADS is done by applying separate scaling factors to the 


objective function, and each design variable and constraint. From examination of the 
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ADS source code [Ref. 13], the design variable and constraint scaling factors were 


found to be applied as follows: 


I) 


4) 
5) 
0) 
7) 
5) 
2) 


Gradients are found using finite differences based upon the initial ADS call 
values. 

The absolute value of the gradient is taken. 

If the gradient is less than the ADS numerical value for zero, the gradient 
value is set to 0.1. 

If the gradient is greater than 1000.0, the gradient value is set to 1000.0. 

The scale factor is set equal to the inverse of the gradient value. 

The product of each design variable and its respective scale factor is taken. 
The product of each constraint and its respective scale factor is taken. 

All scale factors are maintained unchanged throughout the optimization. 
Upon compietion of optimization, each design variable and constraint 15 


áivided by its respective scale factor, to provide a final value that 1s not scaled. 


ror the objective function the scale factor is based upon the square root of the sum of 


ine squares of the design variable scale factors. When scaling is not used, ADS sets 


the scale factors to a value of 1.0. 


3. Scaling in Case 2 


In order to investigate how scaling would affect the problem posed in Case 2, 


the initial design noted in Tables 6 through 8 was assumed to provided a reasonably 


scaled problem. Various strategy combinations were then selected and optimized with 


and without scaling applied. Table 11 shows the results of these runs, and Fig. 6.3 is a 


plot of the state responses for the optimized controller, using initial poles of 
-10.0+ 10.0}, -10.0-10.0;, and -300.0 and the 047 strategy for optimization. | 


Trends that can be seen from these results are: 


In all cases the strategy that had scaling applied found a lower value of the 
objective function. 

All strategies found essentially the same optimum if there was no scaling. 

By scaling, and using the most costly strategy (133) with respect to total 


simulations required, the lowest objective function value was found. 


In Appendix F the output listing for the 057 strategy combination with and 


without scaling can be found. Close inspection of the iteration to iteration 


optimization results reveals that the first iteration is the most critical to the problem 


solution. In the run without scaling, the optimization proceeds slowly, and relatively 
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TABLE 11 


ЕЛ ING RESULTS 





Total | 

Strategy Scaling Objective Final Poles Simulations 

057 11 0492 -20.6 + 66.2] 69 
-20.6-66.2] 
-31.93 

057 No .0658 SR 922] 67 
-11.5-52.2] 
-21.75 

047 Wes 0492 -20.6 + 66.3] 69 
-20.6-66.3] 
ZB с 

047 NO .0640 SL 79 
-11.1-52.8] 
-21.77 

DIT Yes . .0432 - ОВК 15 5 76 
-22.5-75.5| 
-42.57 

957 Мо 0543 -15.4+ 61.9; 80 
-15.4-61.9| 

-25.7 

133 Yes — 0250 -346.4-383. 136 
-346.-383.| 
-366.2 

133 Мо 0645 -9.35 + 54.4) 101 
Бо 


Initial Poles 
-10.0+ 10.0] -10.0-10.0) -300.0 
Initial Objective = 0.1396 


constantly towards the minimum. The gradients are small and thus provide less 
significant data for ADS to process in the determination of search direction and step 
size. When scaling is applied, a significant initial step is made. In fact, the initial step 
“Jumps over” the optimum found in the case without scaling, and then simple 


refinement is done. Thus the classical case of a relative minimum being found is seen. 
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STATE VALUES 





0.00 0.02 0.04 0.06 0.08 0.10 
TIME IN SEC 


| ин 


Figure 6.3 Case 2 Optimal Response 
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It can be concluded that the computational efficiency of turning scaling off is not 
beneficial when using ADSL. 

A word of precaution about the use of ADS print control parameter 
(IPRINT) [Ref. 4:0. 10] is in order. The value assigned to this parameter should be 
selected :0 provide 2 printout that reveals, at least, the optimizer termination criteria. 
Use of [PRINT= 1000, for example, prints only the initial and final optimization 
results. The maximum iteration termination could then be invoked and the user would 
not know that the result obtained from ADS was questionable. While IPRINT = 3050 
produces more output than is really useful, it is valuable for checking what is really 
happening during the optimization process, the type of scaling that occurs, and to 
reveal the complete status of the solution. In general, use of IPRINT= 1010 for cases 
when no strategy is used (ISTRAT=047 or 057) or IPRINT=11!100 for cases where a 
strategy is used (ISTRAT= 133 or 957), proved out to be the most beneficial print 


metro! “aiues. 


ШІ ЖЧП PORE SELECTION 

Different initial poles were selected to determine how initial selection effects the 
optimization process. Intuition may lead one to assume that the optimizer should find 
the same result no matter where the initial poles are placed. This estimate turns out to 
be incorrect on several counts. First, the optimizer 1s likely to find the minimum value 
of the objective function closest to where it initially starts. This is know as a relative 
or local minimum. Without further investigating the whole spectrum of possible values 
for the design variables, the optimizer has no way of differentiating between a relative 
minimum, and the global minimum, which represents the minimum of all relative 
minimums of the objective function [Ref. 3:pp. 12-13]. In Case 2, changing the initial 
poles results in corresponding gains being sent to ADS that can vary drastically. A 
significant change in the initial design then leads to solving essentially a completely 
different optimization problem, that could have its own distinct relative minimum. 
Whether or not this minimum is the global minimum is not known. Secondly, the 
poles selected may result in an initial solution in which one or more constraints are not 
met. This is know as starting in the infeasible domain. ADS will attempt to find a 
feasible solution that will meet all the constraints, but in some cases it 1s unsuccessful. 
In the cases where ADS cannot meet all the constraints, it will still attempt to find the 


best solution, but with one or more violations present. This means that although ADS 
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found a solution, the designer should still try a different set of initial poles in an 
attempt to start ADS in the feasible domain. A shift in the initial pole placement can 
generate a solution that is significantly better, and one that meets all constraints. 
Table 12, shows the results for both small and large changes in the initial guess for the 
system poles using the Case 2 model. Figures 6.4a and 6.4b shows the pole placement 


resuits grapnically (not to scale). 


TABLE 12 
CASE 2 POLE SELEC ORR ESUE TS 


Initial Initial Initial Final Final | Final 


Poles Gains Objective Poles Gains Objective 
-10.04-10.0) ..380 -344.2+ 342.0} + 25.3 
-10.0-10.0| +.600 .139 -344.2-342.0| +264 .024 
Ç £ -300.0 -6.00 -40.88 -962.4 
о -10.0-10.0j -.920 -3.17 + 25.9; 2225 
-10.0- 10.0] 750  .189 -3.17:2588) O17 А 
-30.0 -.600 здано -. 699 
-10.0+ 10.0) +5.02 -75.2 + 150.8 + 2.54 
-10.0-10.0; +14.1 .133 -15.2-150.8| ON 1028 
-3000.0 -60.0 -46.7 159 
1271: Е? -386.5+ 165.3; 4- 20.9 
-100.0-10.0}° +6.01 .026 -386.5-165.3| +3.14 .023 
-300.0 -303.0 -60.4 -642.4 
-10.0+ 100.0} +.610 -248.8 + 210.7} + 12.6 
-10,0-100.0}°  —.600 181 -248 221071 +1.79 .024 
-300.0 -303.0 -60.4 -624.4 


ADS Strategy 047 


Analysis of the results shows that: 
e The assumed initial poles of: -10.0+10.0j, -10.0-10.0j, and -300.0 provided a 
fairly reasonable value for the optimized objective function. 
9 If the complex poles are dominant, nearly the same value for the objective 


function was found for every case of selecting different initial poles. 


5l 





nia Final 


Figure 6.4a Case 2 Initial Pole Selection Results 
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Figure 6.46 Case 2 Initial Pole Selection Results (cont.) 


38 


e There was a wide variety of initial pole combinations that resulted in essentially 
the same value for the final objective function. 


° The average values for the resulting poles are -50, -250+ 150), and -250-150). 


E. CONSPRAINIS 

ADS considers two general types of constraint violations: The first, side 
constraints, are based upon the upper and lower bounds applied to each design 
variable (VUB and VLB). The second, defined constraints, are those relationships of 
the design variables to specific limits, that are developed by the designer to specify the 
system. 

1. Side Constraints 

Througi side constraints, ADS prevents the optimization process from using a 
value less than the specified lower bound or greater than the desired upper bound for 
sacn design variable. The values of these bonds are specified based upon physical 
restrictions ic the system. For Case 2, these restrictions would be derived from the 
physical circuitry used to provide the feedback gains. If at the beginning of the 
optimization, the initial value of any design variable is set outside the allowable 
bounds, the design variable is reset to the closest bound, and a warning message 1s 
generated. It should be noted that this situation can easily occur using the Case 2 
ADSL model. The designer will specify a stable pole combination, and in general will 
not have a good idea of what the corresponding gain values are. In reality the gain 
values will be physically limited to finite values of usually only several orders of 
magnitude. To explore how the optimization is affected by restrictive side constraints 
multiple runs of the Case 2 model, described by Tables 5 through 8, were conducted 
with varying upper and lower bound values. Table 13 shows the results for order of 
magnitude changes Im both the upper and lower bounds. 

Examination of these results shows that by making the side constraints more 
restrictive, a direct effect can be seen on the objective function value. As the 
constraints are tightened, the objective function grows in value, indicating that the 
system has a greater amount of error associated with it over time, and that a larger 
feedback control signal is required to achieve the desired response. | 

As a word of caution, it should be noted that in order to suppress any design 
variables from the model, the programmer must directly reduce the number of design 


variables (NDV). Simply putting the upper and lower bounds (VUB and VLB) to the 
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x 


"әтек Lower Active Final 
,  Dound Bound Constraint Objective 
x + 1000.0 -1000.0 none 0.0233 
| 
| 
É -1000 в; 0.0305 
| +10.0 -10.0 е; 0.0561 
| 
Initial Poles 
-10.0+ 10.0] -10.0- 10.0] 


same value may seem to eliminate the variable, but it will still be taken into account. 


ADS will calculate the corresponding gradients and include them in the optimization 
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with unpredictable results. 


influence the optimization process. It should first be pointed out that ADS keeps 


2. Defined Constraints 


The next area to be considered was how does reducing individual constraints 


tracks of constraints in two ways: 


l) 


2) 


Violated constraints: 


exceeds 0.0 + CITMIN, where CTMIN is defaulted to a value of 0.0] for 


Final 
Poles 


-381 42247 
-38 1-247} 
177 


-300.0 


non-linear constraints, and 0.001 for linear constraints. 


Active constraints: a constraint is considered as active when its value exceeds 


0.0 + CT, where CT is defaulted to -0.03 for non-linear constraints mum 


-0.005 for linear constraints. 
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Final 
Gains 


a constraint is considered as violated when its value 


CTMIN and CT may be modified by changing their corresponding WK array values 
[Ref. 4:pp. 12-14]. 

Recalling that all equality constraints, whether linear or non-linear, must be 
equal to zero, and all inequality constraints must be equal to, or less than zero. the CT 
and CTMIN values can then be considered to generate a tolerance band for ADS to 
work within. Thus ADS never exactly meets an equality constraint of zero, but rather 
finds a value in the band between CT and CTMIN. In the the CONMIN program this 
band was known as the "constraint thickness." 

Constraint thickness comes into play for the Case 2 model with the 
specifications for maximum response (CSP), completion time (TSP), and maximum 
feedback control (USP). If these specifications are made too small, they will essentially 
become lost within the constraint thickness boundaries. In order to look at this effect 
with Case 2 an additional constraint was added after noting the findings about the х» 
(motor current) state in the initial pole investigation. In a manner similar to that used 
for the maximum feedback control signal constraint, a maximum value for state X» 
(X2SP) was established. This constraint was then tightened in order to observe its 
effect on the objective function and the other constraints. For this investigation, the 
more for, and USP valuesiwere reduced to where they resulted in almost active 
constraints for the initial conditions and design described by Tables 5 through 8. The 
exact values are usted in Table 14, and Table 15 documents the 1. of the 


Simulations conducted for reduced values of X2SP. 


TABLE 14 
See CONSTR | SPECIFICATIONS FOR X2SP REDUGTION 


Constraint Initial Constraint Reduction 
Specification Model Model 

CST 140% 22 

TJE |55336 roo see 

US? 10 У MOS 
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TABLE 15 
САЗЕ2”Х2 УР” КМ ОСҮИ МЛ”  РРЭ, 


Final Final Final 
X2SP . Objective Gains Poles Constraints 
2155 -772.0 175Р active 
0775 0.0235 4.23 -215.0 lower bound 
- 1000.0 -60.2 g4 active 
27.99 -Tiu 05р, Сорас 
0.70 0.0235 4 -259.9 lower bound 
-1000.0 -53.3 g4 active 
25.83 -1089.0 A2SP active 
0.60 0.0248 5.53 -189.7 
-996.0 -48.2 
| 24.64 -2146.0 X2SP active | 
0.40 0.0350 10,32 4103 CSP active | 


[nitial Poles 
-10.0 + 10.0; -10.0- 10.0} 
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From the results in Table 15, it appears that reducing the maximum allowable 
value for the x4 state will have varying effects on the other constraints, but in genera] 
will cause the objective function to increase in value. This is as would be expected, 
because with a smaller allowable maximum motor current, the controller will see more 
error ove: time, and require a larger feedback signal to produce the desired output 
response. An additional note of interest is that for every case investigated, the final 
gains resulted in an optimum solution with all poles being real. 

When the other constraints were tightened, individually, similar results were 
seen. [t should be pointed out that once a constraint is pushed too far the model 
collapses. This collapse can be explained if one looks at the design variable space. 


With over restrictive constraints, the feasible design region disappears, and an optimum 


5 


design is impossible. An example of this phenomenon can be seen by looking at 
reducing the completion time constraint (TSP). Tables 16 and 17, respectively, show 


the initial conditions that were assumed and typical results. - 













ЧЕН 16 
то тү 5РЕСТЕТСАТТОМ5 FOR TSP REDUCTION 


Constraint Initial Constraint Reduction 
Specification Model Model 

ESP 140% 5% 

ро Ша. 10А 





Сэ? 10 У 1.05 У 
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TABELCE I7 
CASE 2 BSEC GRON RESULTS 


Final | Final Final 
Objective Gains Poles Constraints 


ДЕС USP active 
lower bound 


(0255 4.23 -215. 
-1000.0 010) g4 active 


-.38 | 
0.6 J TSP violated 
-6.0 i 


Initial Poles 
-10.0+ 10.0) -10.0-10.0) -200.0 
ADS Strategy 047 
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From the results of Table 17 it is apparent that if a constraint is tightened to 
the point where it becomes violated, the optimization will fail to find an improved 
solution. In most cases, ADS will stop after three iterations and return the initial 
solution as the final solution. This fact can be used as an aid when the design process 
is being conducted. The designer would see no improvement for the run and could 
then immediately use this result to relax the constraint specification and rerun the 
model. Typical process time for this cycle, when used with the Case 2 example, was 


one minute (including graphic display time for the IBM 3279 terminal). 


F. INITIAL POLE SELECTION WITH ACTIVE CONSTRAINTS 

In the previous discussion of initial pole selection, none of the proposed 
constraints were active due to the large values of respective specifications that were 
initially proposed. In order to compare the essentially unconstrained results of Table 
12 to a more realistic problem, the Case 2 initial pole investigation was run again, out 
with the constraint specifications of Table 18 being applied. The results of these runs 


are shown in Table 19 and graphically in Figs. 6.5a and 6.5b (not to scale). 
srap y g£ 


TABLE 18 
CASE 2 CONSTRAINT SPECI EI@A ONS Ore NI TIAL ewe 
SEES TION 


Constraint Initial Constraint Reduction 
Specification "^ Model Mode 

CSE 140% 5% 

Х25Р no 0.6 A 

кер 0.85 sec 


USPE 1.05 V 





From comparing the results of Tables 12 and 19, one finds that the тебе 5 


very similar with respect final objective function values. In most cases the constrained 
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TABLE 19 
oe eo) oe ep lA POLE SELECTION RESULTS 


Initial Final Final Final Constraint 
Poles Poles Gains Objective Status 
-10.0+ 10.0; -1089 25.9 22258 
-10.0-10.0) -189.7 3:30 .0248 Active 
-300.0 -48.21 -996 

-10.0+ 10.0] -4.08 + 29.31 -.903 None 
-10.0-10.0) -4.08-29.3] -.906 1225 Active 
-30.0 -10.6 20 

-10.0+ 10.0} -12112:12022) 3.20 Х25Р 
-10.0-10.0) -121.7-130.3j 0.439 „0299 Active 
-3000.0 -44.6 -141.8 

-100.0+ 10.0] -197.6 26:1 X2SP and 
-100.0-10.0} -1070 De .0249 83 lower 
-300.0 -47.27 -1000 ctive 
ШОШ 100.0) -61.7 +39.3) 19.57 21221 
-10.0-100.0) -61.7-39.3] no .0264 Active 
-300.0 -1624 -869.5 
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case values were slightly higher, and this would be expected due to the optimum being 
constrained. In terms of dynamic response, however, the tightening of constraints 1s 
seen to have a considerable effect upon the resultant poles of the optimum designs. 

These results indicate that the sensitivity of the final design to initial pole 
selection remains, in spite of the tightened constraints. It is appparent that for the 
Case 2 system local minima abound in the design spaces that have been investigated 
thus far. Restated in another way, there are many system responses, with or without 
overshoot, that lead to a relatively small value of the objective function. [Factors that 
must be futher explored in an effort to find the global minimum, or at least more 
consistent results would be: 


ө investigation of more exact integration methods 
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Figure 6.5a Case 2 Constrained Initial Pole Selection Results 
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Figure 6.5b Case 2 Constrained Initial Pole Selection Results (cont.) 
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о increased precision of the optimization process through reduced tolerances of 
the ADS control parameters 
° conduct a sensitivity analysis of the various optimization parameters. 
Additionally, selection of a more practical application, together with known parameters 
associated «with component hardware, would no doubt eliminate much of the design 


Пехота charactenimesse acca 


G. NODE INPUT 

The next inquiry into the Case 2 model, as describe by Tables 5 through 8 and 
with the maximum x» state constraint equal to 10.0 Amps, was to extend the initial 
problem to multiple inputs. This simply required changing the w, input from Zero to a 
positive torque loading value, eliminating the negative response constraint, and 
establishing new initial and final, response and feedback control signal values. The 
resulting state response curves, for the assumed initial conditions and a constant torque 
‘oading, are shown in Fig. 6.6. 

In order to see how the final system gains responded to different noise inputs, the 
constant noise torque load (w,) of 1.0 N-m was modified to simulate a more realistic 
load condition. Two different signals were superimposed onto the constant load: 

(1) 1000 Hz sinusoid of amplitude 0.1 N-m 

(2) A random variable from a normal distribution with a standard deviation of 0.1 

N-m 
These inputs were’ generated by single calls using DSL functions and can be found in 
the program listed in Appendix D. 

The results of these three runs are presented in Table 21. Figure 6.7 shows the 
final state responses for the optimized system with a constant torque loading, using 
strategy 133. Observations that can be drawn are: 

° When the constant torque noise loading is modified, with either of the 
superimposed load signals, the resulting imaginary and real poles are reduced in 
magnitude and thus the system response will be slower. 

e The final objective functions are higher in value. This would be expected as the 
system is further in error, due to the applied noise input, and additional energy 
is needed to return the system to the desired state. 

Overall, the application of a noise input demonstrated the robustness of the 


ADSL program. By the addition of one. or two simple DSL function calls, a much 
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Figure 6.6 Case 2 Initial Response with Constant Torque Loading 
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Figure 6.7 Case 2 Optimal Response with Constant Torque Loading 
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Тэ 21 
шог ЁГОГ КЕБОТ.15 


| oise Final Final а Total | 
x Type Poles Gains Objective Simulations 
Constant  -17.75 + 82.6} -.209 ‚399 109 
| Torque -17.75-82.6| -.716 
-21.29 52 | 
iNormal -9.002+ 25.41] -.279 484 23 
21151. -9.602-25.41| -.001 
| -180.6 -13 
| Sinusoid -9.93+3.79; 2379 1.19 19 
| -9.98-3.79] 4-,503 
| -300.8 -6.02 


[nitial Poles 
-10.0+ 10.0) -10.0-10.0; -300.0 


Initial Objective = 1.920 


ADS Strategy 133 


more complicated, and realistic noise input could be investigated by the control 
designer. Any attempt at this tvpe of analysis, using conventional methods, would be 


virtually impossible. 
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VII. MICROCOMPUTER IMPLEMENTATION OF OPTIMIZATION FOR 
CONTROLLER DESIGN 


si 


OVERVIEW 
The ADSL program, discussed in detail in the previous chapters, was run on a 
mainframe computer (see Appendix A for a specific description). The speed and size of 
this type of svstem makes the initial designing of an optimal controller fairly practical. 
Ihe same methodology can, however, be applied to a microcomputer based program. 
One system envisioned would have a microcomputer for the system controller. 
This computer would also monitor the running system, then using stochastic 
processes, generate the time averaged values of the A, B, and C matrices. These values 
could then be given to an ADSL like program that was running concurrently. New 
optimal gains would be computed, and then provided back to update the controller. 
Currently. work is in progress on the development of practical microcomputer 
sontrcllers, but significant problems stil remain. One such problem 15 themes 
obtaining sufficiently high enough sampling rates, while still being able to perform the 


large number of calculations required for a complex controller. 


B. PROGRAM DEVELOPMENT 

It was assumed that a microcomputer based integral controller, similar to that 
used in Case 2 could be developed, as well as the necessary A, B, and C matrix 
generation routines. A program structured to work like ADSL was written for the 
IBM Personal Computer (IBM/PC) using the Microsoft FORTRAN77 
(MS-FORTRAN) compiler and linker programs [Ref. 15]. 

1. Microcomputer Specifics 

The microcomputer selected to implement ADSLPC was the IBM/PC models 

XT and AT, equiped with the Intel 18087 and 180287 math co-processors respectively. 
This decision was based on local availibility, general acceptance as a standard machine, 
and speed of processing. Additionally, the degree of numerical accuracy is directly 
comparable with the IBM mainframe computer [Ref 16:p. 38]. It should be pointed 
out that significant computational time differences were found, with the mainframe 
being able to run the ADSLPC program in approximately 20 seconds while 20 minutes 


were required using the microcomputer. The ADSLPC program generated 
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approximately 380 Kbytes of executable code, and this size required the PC to have 
512 Kbytes of memory, and a hard disk drive for mass storage. 
2. ADSLPC Coding 

ADS was fully implemented using the procedures outline in the ADS 
Reference Manual [Ref. 4]. DSL was replaced with an independent FORTRAN main 
program named ADSLPC.FOR. The designation of ADSLPC was assigned to this 
program as a natural follow-on to the ADSL program for a personal microcomputer 
(PC). The five functional blocks of the ADSL program, described in Chapter 4, were 
closely followed in the design of the ADSLPC program. Control of time was added to 
the main program, as well as integration routines for solving the system of differential 
equations. The following programs from the International Mathematical and 
ЕШ ЕГІСІІ Гтотату (11 М251:) [Ref. 17] were used for determination of the differential 


equation solutions and eigenvalues: 


*  DGEAR 
° DGRCS 

“ DGRPS 

¢ DGRST 

, LEQTIB 
¢ LUDATF 
° LELMF 

е WUERTST 
° EIGRF 

° EBALAF 
° EBBCKE 
° EHESSF 
° EQRH3F 


The ADS code (Version 1.12) was obtained in a format suitable for processing 
by the IBM/PC using the Microsoft compiler. It should be noted that not all 
FORTRAN compilers can be used with this version of ADS, due to the source code 
using FORTRAN Level 66 dynamic dimensioning standards. The IMSL routines 
required some slight modifications to eliminate use of unique IBM extensions to the 
FORTRAN 77 standard, and thus make them compatible with the Microsoft compiler. 


The specific changes are noted in Appendix G. 
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The main program was developed by modifying the translated DSL output file 
(FORT FORTRAN) for the Case 2 model using the eigenvalue method for the initial 
pole to gain determinations. The DSL generated code was first converted to single 
precision expressions, to be compatible with ADS, and the integration routines were 
replaced hy a call to the IMSL subroutine DGEAR. This methology required that the 
derivative block of the program be broken out into its own subroutine, and that the 
Adams integration method be selected. Additionally, the EISPACK routine were 
replaced by the single precision | МЭ ирсоли ЕЕЕ 

One major difference, from the ADSL program, was to incorporate the 
evaluation of constraints and the integration termination criteria at the end of fixed 
intervals (for the Case 2 model this was established at 0.1 secs), rather than at the end 
of every integration step. Figure 7.1 shows the interrelationships of the program calls. 
The complete source code is listed 1n Appendix H, as is a listing of the batch file used 


to control processing of the program. 


C. COMPARATIVE RESULTS 
lable 22 shows how the mainframe and PC program results compare. The 
general conclusions that can be drawn from these results are: 
*  ADSLPC can be used to find optimal controller gains for an integral feed back 
controller, with results comparable to those found using ADSL. | 
e The difference in final objective function values can be attributed to finding 
local minima, which can be most likely traced backed to using a different, and 
less accurate integration method. 
° As with the ADSL model, further improvement of the numerical tolerances 
could result in additional reduction of the objective function. 
ө The obtained results compare favorably with the unscaled ADSL results found 
im [abies]: 
e ADSLPC could be a powerful tool for optimal, online updatn ai 


microcomputer controller processes. 
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TLZ) 
CALL ADS FOR зно VARIABLES 
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Figure 7.1 ADSL Program Segments 
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COMPARATIVE RESULTS FOR ADSL AND ADSLPC PROGRAMS 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


ЫГ CONCLUSIONS 
The following conclusions can be draw from the previous disccussions: 

]l A complicated, constrained, multiple input, multiple output (MIMO) controller 
model can be developed using time domain methods and standard Dvnamic 
Simulation Language (DSL) techniques. This model can then be optimized 
based upon a performance index of choice using the Automated Design 
Synthesis (ADS) program. 

2. the ADSL methodology provides a robust way in which state feedback gains 
can be obtained without use of matrix approximations associated with the 
classical Riccati equation solution. 

3. For MIMO systems in which sinusoidal and random noise inputs were 
simulated, optimal feedback gains could be found. In general they were such as 
to slow the system down, but closelv approximated the gains found for a 
constant noise input. 

4. The ADS optimization code, with an eigenvalue solver routine, can be 
emploved to solve the problem of finding the system poles corresponding to 


given gains for anv complex MIMO controller. 


Cn 


From the studies described herein, the following approaches were determined to 

be well suited for developing ADSL programs: 

e Use a stiff integration method 

ә For inital pole to gain determination, if an exact determination 15 
impractical, use ADS and an eigenvalue solver method as developed in this 
study. 

е For basic controller optimization, application of the Method of Feasible 
Directions strategy (047) proved to be the most practical overall, but the 
Sequential Unconstrained Minimization Techinque using an external 
penalty function strategy (133), in general, provided better optimums at the 
expense of more iterations. 

өе Order of magnitude changes of the initial poles proved to be a reasonable 
method for establishing the validity of the optimization results for the final 


gains. 
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e If the final solution has active constraints, these should be investigated 


individually, by reduction of the corresponding constraint specification. 


6. For the state feedback controller considered in this study, many local optimums 
can be found. The reasons for this appears to line in the definition of the 
performance index and the relative freedom from severe design constraints. 

7. The basic ADSL program methodology can be down scaled to run on a 
personal computer, thus opening the possibility to using optimization in on-line 


applications. 


B. RECOMMENDATIONS 
The following recommendations are offered: 

l. A simplified state feedback model should be investigated as the results obtained 
for the integral state feedback model demonstrated that the selected design 
space had many local minima. An acutal controller, with known physical 
parameters, would likely have a well behaved design space. This would allow 
for gaining a clearer understanding of the interactions investigated in this study. 
Additionally, the results of actual tests could then be applied to verify the 
ADSL predictions. 

2. А large scale controller model should be optimized. One such controller model 
has already been developed by NASA and has been documented by Merrill 
ЕКИ 

3. The methodology should be expanded to provide for online interaction, and 
recursive optimization capability. 

4. The same methodology used with ADSL can be applied to a micro-computer 
based design system using trimmed down codes. One possible system would be 
replacement of DSL by the IBM PC-Engineering Simulation Program and ADS 
With the Engineering Design Optinuzation program Microdot, which is a 


reduced version of the ADS code in a format for use on a micro-computer. 
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APPENDIA A 
гс:  УТЕВБЁАСЬЕ5РЕСТЕЇС5 


1. COMPUTER DESCRIPTION 

The ADSL program was initially developed at the Naval Postgraduate School W. 
R. Church Computer Center using a IBM 3033 Attached Processor System (16 
Megabytes) loosely coupled with a IBM 3033 Model U (16 Megabytes) and a IBM 
4381 Model M1 (8 Megabytes). Interactive computing was provided under VM/SP 


CMS, and batch-orocessing under MYS with JES3 networking. 


2, SPECIFICATIONS USED FOR COMPILING ADS (VERSION 1.19) 
MOKI RAN V5 
o. Auts-Double Precision (DBL) 
ӨГ тш 1011 Level 


d. Language Level 66 


3. DSL/VS (VERSION 1.1) IN THE VM MODE 

The resulting TEXT file from the compilation of ADS was converted to a 
FORTRAN library fle using the online library generation procedure [Ref. 19:pp. 
35-36]. The library was named MYLIB TXTLIB in order to take advantage of the 
default DSL user library specification for the VM environment. 

The default DSL executive was modified to streamline the interfacing of of ADS 
ате DSL output, and is shown in Appendix B as the ADSL EXEC. 


4. DSL/VS (VERSION 1.1) INTHE MVS MODE 

The ADS TEXT file was also placed onto the Mass Storage System (MSS) as a 
user library in order that DLS/VS could access ADS using the MVS environment. The 
standard job control statements for using DSL must then be modified to concatenate 
the ADS, LINPACK, and EISPACK libraries as shown in Appendix B. The job 


control for compiling ADS on MSS is shown below. 
5. JOB CONTROL FOR FILING ADS ON THE MASS STORAGE SYSTEM 


//FILEADS JOB (0180,9999),’FILEADS’ ,CLASS=C 
// EXEC FORTVCL,;,PARM.FORT-'LVL(C66),ADCDBL ) ,NOMAP ; OPT( 3 ) , NOS ,NOX"' 


//FORT.SYSIN OD ж 
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x3exx9336eo* REPLACE THIS LINE WITH THE ADS SOURCE CODE 39eeeeeeeeeeeee 


/% 
//LKED.SYSLMOD DD UNIT=3330V ,MSVGP=PUBSCC , 
2 05М-М55.Ғ0180.А05(А05), 


Ж DISP=(NEW,CATLG, DELETE },SPACE=(CYL,(4,4,5)) 
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APPENDIX B 
STEPS TO RUN ADSL AT NAVAL POSTGRADUATE SCHOOL 


1 VM ENVIRONMENT 
a. Logon in the VM environment (note that the user must have have previously 
| requested an increase of his default virtual machine size from 750 Kilobytes to at least 
1500 Kilobytes through the Computer Center Accounting Office). 
b. Define increased storage by using: 
р. 02208 1500К 
| Г ЄМ 
c. Linkto the optimization disk with the command: 
l I [O O0JSOP IJ OPTIM 
d. Link to the DSL and MATHPACK Disks with the command: 
ӘЛІГЕ 
This command exccutes the DSLINK executive procedure described below. 
e. Run the ADSL program using: 
ADS ЕТ | ТЕК» 


ЕРЕ (Пе user provided filename containing the DSL input program. The 
filetype of the input program must be DSL. PLOT or TEK is optional, but required if 
any plotting 15 called. TEK 1s used when TEK-618 or IBM 3279 terminal plotting 1s 
desired, while PLOT 1s used for all other plotting devices. Also note that the user can 
not have a file name MAIN TEXT on his user disk, as this file will cause DSL to load 


abnormally. 


f. The specified tabled output will be saved in a file named FN LISTING. 
g. Plotting will then be done on the specified graphic output device and printer 
proce! be concatenated to the FN LISTING file. 


2. DSLINK EXEC LISTING 


EXEC LINKTO DSL 
EXEC LINKTO GDDM40 191 GDDHMR 
EXEC LINKTO PLI5 191 PLI5 


EXEC LINKTO 0062P 191 EISPACK 
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3. ADSL EXEC LISTING 


"/^ Install REXX 1f not here. 
&CONTROL OFF NOMSG 
EXEC REX 1 
&IF SRETCODE * 0 EXEC 0 &1 82 83 86 85 86 27 88 89 810 211 812 813 816 815 816 817 818 


ЗЕХЇТ @RETCODE 


(end of REX installation) м / 


/Імммннниминнимниининимимиминнм цими мичнимимим м минмимнинимнинининимимимними мими / 


/и +++ REXX EXEC TO PERFORM COMPLETE DSL SIMULATION “/ 
/ м “/ 
/" **** DSL/VS SYSTEM ... ITEM Corporation 1986 "/ 
/" **** LATEST UPDATE: 4/03/84 “/ 


/(Мминининимпипимининнимийнийничинмийни йнминийни ийн мийн ниминийни йн иммийинн/ 

Trace value 'OFF' 

/*ARG FN FT FM DISP "('" GRAF REST»*/ 

ARG FN GRAF 

эсгэн 

TELE: Do 
Say ‘The ADSL command makes a complete DSL run in the foreground. Бу! 
Say ‘executing the following EXECS in order: SLT. DSLC. and DSLG.' 
бау ‘The plot post-processor GRAFAEL may also be invoked at the same' 
Say ‘time as an option to produce graphic output at a selected graphic’ 
Say ‘output device.’ 
Say ' ' 
SAY ‘COMMAND FORMAT : ADSL «FN» «PLOT» ' 
Say ""ЕМ" 1s the filename of a previously created DSL program file.' 
Sey ' ' 
Say 'NOTE that the filetype must be DSL.' 
SAY о 
Say '"PLOT" ís an optional parameter which specifies if graphic post-’ 
Say ‘processing 1s desired. Leaving 1t out. means no plotting is done.' 
SAY ‘Use the command "PLOT" for printer plots or “TEK” for any other’ 
SAY 'displays.' 

End . 

IETEN з "" ТНЕМ ЕЙІТ 

/м **»* invoke DSL Translator 4/ 


FT='DSL' 
FM='A' 
DISP='FILE' 


/* GRAFz'GRAFAEL' */ 
'"EXEG DSUT' fn ft fm 
lfere эз 9 then exit Fc 
COMPILE: 
/м ++++ invoke Fortran Compiler #/ 
Do 
SAY ' >>>>>>> RUNNING VS FORTRAN COMPILER <<<<<<' 
End 
'FORTVS FORT (LANGLVL (77) OPT(3) NOXREF NOMAP NOSDUMP NOTF NOTERM' 
/*5::: "FORTVS" INVOKES THE VS FORTRAN COMPILER (PROG. 5 5748-FO3) :::1«/ 
/"::: Txtlibs assumed: | VALTLIB  VFORTLIB sions 
IF RC >= 8 THEN 
Do /* Type Error Diagnostics "/ 


Say є 1 

Say бийд: инцнинимиенйнининниминмечичинмийнинмницийнининйиицчини чийнийиййниийинчиминн” 
Say ' z2£s52 ERRORCS) DETECTED IM COMPILATION PHASE паша” 

Say ' ' 


Say ' Fortran listing will be displayed by XEDIT: Enter “OQUIT™ to ехіє" 
Say * bg Sd OF Ad OE E м н мм Od Oe м” 
Say ' ' 

Queue 'PREFIX OFF' 

'XEDIT FORT LISTING' 

ЕХІТ 22725 

End 


qt 


419 820 


IF RC >=B THEN EXIT 
/4 ee Test for axistence of user's MYLIS TXTLIS file “/ 


Е ТАТЕ MYL.S ҮХТЕТВ” /“%:: МҮ(ІЗ пау be renamed :::«/ 
ГЕ ГС a Q Име: verify tne VS Fortran Library 
/*t:: names assumed: VALTLIB 82 VFORTLIB :::«/ 


/ 4 АШ20 .iNKTO ZISPACK AND LINKPACK LIBRARIES »/ 
БЕНО 22АГ anions 252518 2522 ХК5 VALTIL IB. VFORTILIBWEAMSL IB. ' 
ЗОГС Гц, ЭРЭЛ” 
else ‘GLOBAL :XTLIS MYLIG' /4::: MYLIS may ое renamed :::0W/ , 
°OSLSIM DSLSLKS VALTLIB VFORTLIB CMSLIB, ° 
'ZISPACK LINPACK' 
Z erase old "points" file 1f it exists., ”/ 
Enc (LCT DATA' 
If rc = 0 then "ERASE PLOT DATA’ 
'FILEDEF 06 DISK QUTP DATA (RECFM FBA LRECL 133 BLKSIZE 665! 
'FILEDEF 05 DISK AUX BATA’ (RECFM FB LRECL BO BLKSIZE BOO' 
'FILEDEF 12 DISK PLTC OATA (RECFM FB LRECL BO BLKSIZE 800° 
'FILEDEF 13 OISK TABL OATA (RECFM FB LRECL BO BLKSIZE B00' 
'FILSDEF I4 DISK PLOT SATA CRECFM VB BLKSIZE 1024° 
K 2ATA DATA (RECFM FB LRECL 80 BLKSIZE BOO’ 


EXSIDPEDEF 15 215 
'FILEOEF i5 2ÍISK ?RTSAV DATA (RECFM VBS LRECL 540 SLKSIZE 544" 
ENDE ot solos "UM DATA 4«ECFM 7B LREC. 89 BLXSI-CE 800” 
/* ::: 3dceitisral user *ile derinitions you may Se added here :::*/ 
Do 

SAY 


GAY ' »»»»»»» RUNNING JSL/VS SIMULATOR «««««««' 
спід 
зум; 
Биг 21 SMSGaGSE' 
*SLOBAL LOAOLiB VFLODLIB* 
‘LOAD FORT (CLEAR RESET MAIN START NOMAP* 
samme гс 
SERESET EMSG TEXT: 
If simre >= B then 


Do 
зам" ' 
СУМ n o MI MEM MN ME MN OUmmM M KL M ONRMDM DOM HO NN ORQICOR0NTO) (OG 9 н и ' 
Say ' teeen ERROR(S) DETECTED IN SIMULATION PHASE zzsuug' 
SY * WE a0 De De ad Oe De add be ad FEO Oe mm OE DOE Y MES M Oe OOO na 
Say ' ' 

End 

CONT : 


PGreaGteliisG ON’ 

ШҮЗІПГЕ эс: B then Exit 55553 

IF GRAF » 'PLOT' THEN 'GRAFAEL.' 

IF GRAF s 'TEK' THEN 'GRAFAEL' 

'ERASE 'FN ' LISTING A" 

IF DISP s 'PLOT' THEN 'COPY OUTP DATA A PPLT OATA A ' FN * LISTING A’ 
ELSE 'RENAME OUTP DATA A ' FN ' LISTING A* 

'X * FN ' LISTING ' 

Exit /* end of program w/ 


4. SAMPLE JOB CONTROL FOR MVS ENVIRONMENT 


ZZADSL JOB (0180,9999),"АО51!,С1А55:8 

КИЕТ. EXEC DSLVS,PARHM.C='NOMAP ,NOX,NOL ,NOS ,OPT( 3 )' 
Z//ZX*FORMAT PR,DDNAME -JESMSG COPIES-0 

//ZXFORMAT PR,DDNAME-SYSMSG,;COPIES-0 

/ZXFORMAT PR,DDNAME-JESJCL ,COPIES-0 

ZZBEST EXEC DSLYS,PARM.C='NOMAP ‚NO:2, NOL ,NOS,OPTI{3 )' 


// DSL. INPUT DD ж 
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SN 


%%%%* INSERT ADSL PROGRAM CODE HERE жжжж 


END 


STOP 


КАЛД Ез су сезып 


iS. 00 

7274 DD 

227 00 

БА А DD 

СА DD 

204 | 00 DISP=SHR,DSN=MSS.F0180.ADS(ADS ) 

22 DD DISP=SHR,DSN=MSS.SYS3.EISPACK. LOAD 
77 DD DISP=SHR,DSN=MSS.SYS3.LINPACK. LOAD 
/% 

22 
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APPENDIX C 
EXAMPLE SISO PROGRAM LISTING 


TITLE SISO 2ND ORDER SYSTEM-OGATA SECTION 7-3 WITH CONSTRAINTS 


3€3€3€ €3€€X9€ 3€3€ 2€3€ 9€ € 3€ 3€ € € 3€ € 9€ HE EEE IEEE IE EE EE KE FE IEE EE EE IE EEE EE HE EE ЭС ЭЄЭС ЭЄЭС ЭС ЭСЭС ЭЄЭСЭС ЭХ HEE EE EE HE HEE 


x LISTING OF ALL VARIABLES 


x 


x 


DF 

G 

H2 

I 

ІС 

IDG 
IGRAD 
INFO 
TONED 
ТОРТ 
IPRINT 
ISTRAT 
IWK 
LSTOBJ 
LSTTIM 
NCOLA 
NCON 
NDV 
NGT 
NRA 
МВІНК 
МВНК 
OBJ 

R 


TAU 


WORK MATRIX FOR ADS 


VECTOR 


VECTOR 


OF DESIGN VARIABLE GREADIENTS 


OF CONSTRAINT VALUES 


VALUE OF THE RESPONSE OSCILATION 


LOOP COUNTER 


VECTOR 


VECTOR 


SCALAR 


SCALAR 


SCALAR 


SCALAR 


SCALAR 


SCALAR 


OF CONSTRAINT GRADIENTS 

OF CONSTRAINT TYPES 

FOR INDICATING METHOD FOR OBTAINING GRADIENTS 
INDICATING OPTIMIZATION STATUS 

INDICATING TYPE OF ONE-DIMESIONAL SEACH 
INDICATING TYPE OF OPTIMIZATION 

INDICATING TYPE OF PRINT CONTROL FOR ADS 


INDICATING TYPE OF STRATEGY 


INTEGER HORK VECTOR FOR ADS 


VALUE OF THE OBJECTIVE FUNCTION AT PREVIOUS INTEGRATION STEP 


VALUE OF TIME AT THE PREVIOUS INTEGRATION STEP 


SCALAR 


NUMBER 


NUMBER 


SCALAR 


SCALAR 


SCALAR 


SCALAR 


COLUMN DIMENSION OF A MATRIX 
OF CONSTRAINTS 

OF DESIGN VARIASLES 

ADS CONTROL VARIABLE 

ROM DIMENSION OF A MATRIX 
DIMENSION OF IWK VECTOR 


DIMENSION OF WK VECTOR 


OBJECTIVE FUNCTION VALUE 


VALUE OF THE INPUT TO THE SYSTEM 


SYSTEM TIME CONSTANT 
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* TIME DSL TIMER VALUE 


x МІВ VECTOR OF DESIGN VARIABLE LONER BOUNDS 


ж МОВ VECTOR OF DESIGN VARIABLE UPPER BOUNDS 


x WK ADS WORK VECTOR 
ж X VECTOR OF DESIGN VARIABLES 
ү VALUE OF THE SYSTEM OUTPUT 


ж ҮМАХ MAXIMUM VALUE OF THE SYSTEM OUTPUT 


ж YOSC VECTOR CF SYSTEM OUTPUT VALUES AT PREVIOUS INTEGRATION STEPS 


з 36363638 3836 33636 3636 JEDE JE JEJE JE JEJE JE JE JE HEHEHE JE JE JE IE € 9€ € 9€ € 9€ JE JE HEHE EEE 3€ 9€ € 3€ 9€ 9€ 3€3€ 3€ €9€ € HEHE HIE HH HHH HHH HHH 


x PREAMBLE TO INITIALIZE THE ADS PARAMETERS 
FIXED ISTRAT,;IGPT,;,IONED,;IPRINT,INFO,IGRAD 


FIXED NDV;NCON,;IDG,;NGT;IC,;NRA,;NCOLA,;NRNK;IWK,NRIHK,I 


* DIMENSIONING REQUIREMENTS ARE FOUND IN THE ADS MANUAL TABLE 5 
X ARRAYS MUST BE ASSIGNED USING THE DIMENSION STATEMENT IN DSL 


D DIMENSION A(2;,6) 


* 


VECTORS CAN BE ASSIGNED USING THE ARRAY STATEMENT IN DSL 
ARRAY WK( 1000 ),INK(500),DF(2) 


ARRAY X(2),VLB(2),VUB(2),6(5),IDG(5),IC(5) 


* INITIALIZATION OF ADS INPUT CONSTANTS (ADS MANUAL TABLE 5) 
PARAM МБА-1 ,МСОГ1А-1 ,МКМК-1000 ,МКТИК-500 


PARAM IGRAD=0,INFO=0 »NDV=1,NCON=5 


* SETTING ALL CONSTRAINTS TO NON-LINEAR, INEQUALITY 


TABLE IDG(1-5)=0,0,0,0,0 


* INITIAL GUESS FOR DESIGN VARIABLE ZETA 


TABLE X(1J=1. 


* UPPER AND LOWER BOUNDS FOR DESIGN VARIABLE 


TABLE VLB(1)=.01 , VUB(1)=2.0 


x SET ADS OPTIMIZATION STRATEGY (ADS MANUAL PART 4.6) 


PARAM ISTRAT=1, IOPT=3, IONED=3, IPRINT=2020 


*SET ADS PRINT PARAMETER (ADS MANUAL TABLE 5) 
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PARAM IPRINT=1000 


x SET DSL INTEGRATION OUTPUT PARAMETERS (DSL MANUAL PART 4.2) 
ewig STIrr 

CONTROL “ІМТІМт85 

SAVE /.0BJ 

GRAPH (DE=TEX618) TIME,Y ,OBJ 


PRINT Y,08J 


EIEH JEJE EJE JE JE JEME JEJE JE JE JE JE JE JE JE HE JEME JEJE JE E JE JE JE JE JE JE JE JE FE HE FE JEJE JEJEJEJE E E ванню аж 
INITIAL 


x CALL TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 
ж CALLING STATEMENT SPECIFIED IN ADS MANUAL PART 6.0 
CALL ADS(INFO,ISTRAT,IOPT,;iO0NED,;iPRINT,;IGRAD,;... 
NDV ,NCCN 4, VLB , VUB,CBJ,G, IDG,NGT IC,... 


DF A; NRA;NCOLA , NK ;NRWK , INK ,NRIMK ) 


x THIS STATEMENT TURNS ON THE PRINTER OUTPUT ONLY ON THE LAST RUN 
x NHEN INFO HAS BEEN RESET TO 0 ВУ ADS TO INDICATE OPTIMIZATION 
X HAS BEEN COMPLETED 
ІҒ (ІМҒО.Е0.0) DELPRT=1.5 
*THIS STATEMENT TURNS ON THE RECOROING OF PLOTTER VALUES 


IF (INFO.EQ.0) DELPLT=.4¢ 


X RE-INITIALIZE CONSTRAINT CHECK VALUES AFTER CALL TO ADS 
YMAX=0. 
TAU=0. 
(5Т08./50. 
LSTTIM-O. 
YOSC1-0. 
YOSC220. 
YOSC3=0. 
DO 15 i-1,NCON 


15 G(I)-1.D20 


FEE IE FE FE EH HEHE IESE HEE ЭЄЭЄ ЭЄ3© ЭЄЭЄ ЭЄЭЄ IESE FEE IESE FE IE FE IE FE HE ESE HE SESE SE SE IE 56 IE IESE FE EEE € € XC 3€ € € 3€ € € CC € X 
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DERIVATIVE 


NOSORT 


* INPUT TO THE MODEL (DSL MANUAL SECTICN 5.1) 


coset eo.) 


* OUTPUT OF SECOND ORDER LAG SYSTEM (DSL MANUAL SECTION 5.4) 


YsCMPXPL(O.,0.,XC 1), 1.,R) 


ж ІТАЕ PERFORMANCE INDEX TO MINIMIZE AMPLITUDE (OGATA P. 299) 


OBJ-INTGRL(O. TIMEXABS(R-Y )) 


X393 € 3€ HE JE 3€ 9€ € X 9€ 3€ X Y ME JE 2€ 36 HEE JE JE € € 9€ 9€ € 3€ 363€ JE JE 9€ 3€ 3€ JE JE JE JE JE JE JE 9€ 3€ 36 36 3€ 3€ 3€ 9696 3€3€ 9€ 96 9€ 2€ 2€ 9€ 
DYNAMIC 


* UPDATE CONSTRAINT CHECK VALUES AT EACH INTEGRATION STEP 
1OSCL-YOSC2 
" YOSC2-YOSC3 
YOSC3=Y 
H2-ABS(YOSC2-R) 
IF(TAU.NE.O.)GOTO 20 
IF(Y.GE.(0.65212%R))TAU=sTIPFE 


20 IF(Y.GT.YMAX) YMAX=Y 


* MAX AMPLITUDE CONSTRAINT ( <1407) 


IF(Y.EQ.YMAX) G(1)72YMAX-(C1.4*R) 


x MAX DELAY TIME CONSTRAINT ( «1.5 sec) 


IPFCCY.LE.C.5*R)). AND. CY. EQ. YMAX) ) GC 2)2TIME-1.5 


x MAX RISE TIME CONSTRAINT ( «23.5 sec) 


IF(CY. LE. RO. AND. CY. EQ. YHAX)) GC S3 )ZTIME-3.5 


ж MAX PEAK TIME CONSTRAINT ( «4.5 sec) 


ІҒ(Ү.ЕС.ҮМАХ) б(4):11МЕ-4.5 


x MAX COMPLETION TIME CONSTRAINTS ( <12 sec) 
IF(G(5).LT.0.)GOTO 30 


* FOR UNDER DAMPED CASES 
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IF((H2.GT.ABS(YOSC1-R)).AND.(H2.GT.ABS ... 
I OSC2-F)) AND. (H2.LE.( .Dl1xR11) 6(5)-ТІМЕ-12. 


FOR OVER DAMPED AND CRITICALLY DAMPED CASES 


* 


IF(((Y-R).LT.(.01%R)).AND.(TIHE.GE.(5.%TAU))... 


AND.(TAU.NE.0. )J.AND.(Y.EQ.YMAX)) G(5)=TIME-12. 


INTEGRATION COMPLETE CHECK USING THE OBJECTIVE FUNCTION SLOPE 


* 


* 


CHECK FOR FIST INTEGRATION STEP 
30 IF((OBJ.LE.0.).OR.(LSTOBJ.LE.0.))GOTO 40 


DETERMINE SLOPE AND CHECK FOR PAST COMPLETION TIME - IF MET 


* 


* 


ENORUN DIRECTS PROGRAM TO THE TERMINAL SEGMENT OTHERWISE THE 


* 


NEXT INTEGRATION STEP IS TAKEN UNTIL FINTIM IS MET 
BEUCUDOBJSESTOBJJZCTIME-LSTTIM).LT.1.E^75).AND. 
(((H2.LE.(.01%R)).AND.(H2.GT.ABS(YOSC1-R)) 
mD. (HP .GT.ABS(Y0SCSR))).OR.( (TIME. 
GE.(5.xTAU)).AND.(TAU.NE.0.))))CALL ENDRUN 

40 LSTOBJ-OBJ 


LSTTIM=TIME 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX% XXX 
TERMINAL 


x CHECK TO OPTIMIZATION COMPLETE TO EXIT OSL OTHERWISE RERUN WILL 
X START THE PROGRAM OVER AT THE INITIAL SEGMENT 
IF (INFO.EQ.0) CALL ENDJOB 


CALL RERUN 


END 


STOP 
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APPENDIX D 
STATE VARIABLE PROGRAM! LISTINGS 


TITLE DC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 


x LISTING OF ALL VAPIABLES 


* 


* 


AA 


BB 


CC 

СЕ 

CI 
CON 
СЭР 
DELC 
DELU 
DELX2 


ОУ. 


FV1 


IC 
IDG 
IER1 
IERR 
IGRAD 
TI 
INFO 
TONED 


ТОРТ 


ЇРВЇМТ 


MATRIX OF STATE COEFFICIENTS 

VECTOR OF CONTROL COEFFICIENTS 

SCALAR OUTPUT VALUE 

MATRIX OF STATE OUTPUT COEFFICEINTS 

FINAL EXPECTED OUTPUT VALUE 

INITIAL OUTPUT VALUE 

VECTOR OF CONSTRAINT VALUES 

MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 

DIFFERENCE OF THE ACTUAL OUTPUT TO THE INITIAL OUTPUT 
DIFFERENCE OF THE ACTUAL CONTROL. TO THE INITIAL CONTROL 
DIFFERENCE OF THE ACTUAL X2 STATE TO THE INITIAL X2 STATE 
VECTOR OF DESIGN VARIABLE GREADIENTS 

VECTOR OF INPUT COEFFICIENTS 

EISPAC WORK VECTOR 

VECTOR OF DESIGN VARIABLE GAIN VALUES 

VECTOR OF OUTPUT COEFFICEINTS FOR THE INPUTS 

LOOP COUNTER 

VECTOR OF CONSTRAINT GRADIENTS 
VECTOR OF CONSTRAINT TYPES 
LINPAC ERROR INDICATOR 
EISPAC ERROR INDICATOR 
SCALAR FOR INDICATING METHOD FOR OBTAINING GRADIENTS 
COUNTER OF NUMBER OF FUNCTION EVALUATIONS 

SCALAR INDICATING OPTIMIZATIOM STATUS 

SCALAR INDICATING TYPE OF ONE-DIMESIONAL SEACH 
SCALAR INDICATING TYPE OF OPTIMIZATION 


SCALAR INDICATING TYPE OF PRINT CONTROL FOR ADS 
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IFVT 


ISTRAT 


LSLOPE 


LSTOBJ 


ШӘТТІМ 


МАТ2 


МЕАМ 


* NCOLA 


NCON 


MOV 


NGT 


© NM 


NN 
NRA 
МВІНК 
МЕК 
ов.) 
008Ј 
P1 

P2 

P3 


РА 


SD 
SEED 
TAU 
те 
ТІМЕ 


TSP 


UI 


LINPAC 


SCALAR 


ЕЗЭРАС 


INTEGER HORK VECTOR 
INDICATING TYPE OF STRATEGY 


INTEGER AORK VECTOR 


ІНТЕЗЕТ 


AORK “ECTCR FOR 405 
(СОР COUNTER 
OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 


SLOPE OF THE 


VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 


VALUE OF TIME AT THE PREVIOUS INTEGRATION STEP 

EISPAC SCALAR VALUE TO INDICATOR EIGENVALUE DETERMINATION 
DSL NORMAL DISTRIBUTION MEAN 

SCALAR COLUMN DIMENSION OF WA MATRIX 
NUMSER OF CCNSTRAINTS 

OF DESIGN VARIABLES 
ADS CONTROL VARIABLE 

DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 
DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 
SCALAR ROW DIMENSION OF WA MATRIX 
SCALAR DIMENSION OF IWK VECTOR 
SCALAR DIMENSION OF WK VECTOR 

OBJECTIVE FUNCTION VALUE 

OBJECTIVE FUNCTION FOR INTEGRATION 

DESIRED INITIAL COMPLEX POLE VALUE 

DESIRED INITIAL COMPLEX POLE VALUE 

DESIRED INITIAL COMPLEX POLE VALUE 

MATRIX OF INITIAL GAIN DETERMINATION COEFFICIENTS 
STATE WEIGHTING MATRIX 

CONTROL WEIGHTING VECTOR 

STANDARD DEVIATION OF DSL NORMAL FUNCTION 

RANDOM MUMBER SEED FOR DSL NORMAL FUNCTION 
SYSTEM TIME CONSTANT 

TIME WHEN COMPLETION TIME SPECIFICATION IS MET 
DSL TIMER VALUE 

MAXIMUM COMPLETION TIME 

SCALAR FEEDBACK CONTROL VALUE 


INITIAL VALUE OF THE CONTROL VECTOR 
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x UMAX MAXIMUM CONTROL VALUE 


* USP MAXIMUM VALUE OF THE CONTROL VECTOR 


x UT TRANSLATED FEEDBACK CONTROL VALUE 
к MATRIX FOR FINAL POLE DETERMINATION USING EISPACK 
4253 "20708 OF DESIGN VARIABLE LOWER BOUNDS 


* VUB VECTOR OF DESIGN VARIABLE UPPER BOUNDS 


x WI DISTURBANCE VALUE 

ж AZ REFERENGE SETTING 

ж МА. ACS WORK VECTOR 

ж НІ VECTOR OF FINAL IMAGINARY POLES 

* WL CIRCULAR FREQUENCY OF SINUSOIDAL DISTRUBANCE 
x WR VECTOR OF FINAL REAL POLES 

x Xxl MOTOR SHAFT SPEED STATE VALUE 


= KIT TRANSLATION OF Xl STATE 

ооо MOTOR CURRENT STATE VALUE 

ж X2I INITIAL VALUE OF THE X2 STATE 

* X2MAX MAXIMUM VALUE OF THE X2 STATE 

* X2SP MAXIMUM SPECIFIED VALUE OF THE X2 STATE 

ж ХХІ FUNCTION DESCRIBING THE Xl STATE FOR INTEGRATION 

ж Хх2 FUNCTION DESCRIBING THE X2 STATE FOR INTEGRATION 

x YOSC1 VALUE OF SYSTEM OUTPUT OSCILATION AT TWO PREVIOUS INTEGRATION STEPS 
* YOSC2 VALUE OF SYSTEM OUTPUT OSCILATION AT PREVIOUS INTEGRATION STEP 
* YOSC3 VALUE OF SYSTEM OUTPUT OSCILATION 

г EISPACK DUMMY SCALAR 

хх << << 5ЭЭ 55555 << HAI IEE HEHEHE ЭЄЭС-ЭС ЭЄ ЭЄ ЭСЭСЭЄ ЭЄЭС ЭС ЭС ЭС ЭСЭСЭ6 ЭЄЭСЭЕЭСЭСЭСЭ6 
%* INITIALIZE THE ADS PARAMETERS 

FIXED NRA,NGT,NCOLA ,NRAK >NRIWK ,IGRAD ,NDV »NCON,IDG ,IC , ITWK 

FIXED ISTRAT ,IOPT , TONED ,IPRINT ,INFO,IPVT 

ARRAY DV(4),VLB(4),VUB(4),CON(6),IDG(6),IC(6),DF(4),WK(1ODD) ,TINK(1D0D ) 
D DIMENSION WA(4D,4D ) 

PARAM NRA=4D,NCOLA=40 ,NRWK=1000 ,NRIWK=1DD0 »NDV=3 »NCON=6 , IGRAD=D , INFO=D 
* ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

TABLE IDG(1-6 )=3*D,2*2,D 

* LOWER BOUND FOR DESIGN VARIABLES 


TABLE VLB(1-3)-23*-1DDD. 
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% UPPER BOUND FOR DESIGN VARIABLES 

TABLE VUB(1-5 }=3%100C. 

ж А05 CPTIMIZATION METHOD SELECTION 

PARAM ISTRAT=0, 1ТО?РТ-4, IONED=7, IPRINT=1010 

3€ 3€ FEI IEEE EEE HE EEE юю З юю юю ее зе зе 
т AHUIILIALIZE THE DSL PARAMETERS 

PRINT CBJ;,C,X1,Xx2,U 

SAVE ( RUN) OBJ,X1,X2,U,C 

CONTROL DELPLT=0.005,DELPRT=0.01 

CONTROL FINTiM-54. 

METHOD STIFF 

RELERR ӨБЈ=1.=-4 

АВЗЕКК 083:1.Ё-4 

зе ее а з жа зе 36-36 з 36-36-36-36-36-36-36-36 36-36-36-36-36 ЭЭС ЭЕ ЭЭЭ ЭЕ ЭС ЭЕ ЗЕ ЗЕ Ж 36-36-36-36-36 36-36-36-3636-36-36-3-36:36 36 
x INITIALIZE THE MODEL PARAHETERS 

ЕТХЕр 1,11,1ЕВ1,ТЕВБ,ММ,ММ,ГТУ1,МАТ2,5ЕЕОр 

COMPLEX P1,P2,P3 

ARRAY ВВ(С3),Н(2),В(2),1РУТ(3 ),66(5 ),МВ(3 ),нН1(3),1У1(3 ),РУ1(3) 

р ОЇМЕМ5ЗЇОМ АА(3,3),СС(2,2),84(2,2),РА(3,3),Е(2,2),У(3,3),2(3,3) 
x ZERO ALL MATRICES 

D DATA AA,F,CC,Q/21*0.7 

TABLE 88(1-3 )=3*0. »R( 1-2 )=3%0. »H( 1-2 J=3X0. 

% SET CONSTRAINT SPECIFICATIONS 

* CISINITIAL OUTPUT VALUE; CF=FINAL EXPECTED OUTPUT VALUE 

*X CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 

*% TSP=MAXIMUM COMPLETION TIME 

* USP=MAXIMUM VALUE OF THE CONTROL VECTOR 

PARAM CI=1. ,CF=0. ,CSP=5. ,TSP=1.00,USP=1.05,X2SP=1.0 

* SET INITTAL INPUT VALUES (W1=OISTURBANCE, W2=REFERENCE SETTING} 
PARAM W1=0., W2=1., SEED=113, MEAN=1., SD=0.1, WL=157 

* SET INITIAL RUN FLAG AND EISPACK PARAMETERS 

PARAM II=0,MATZ=0 »NM=3 ,»NN=3 ,I=1,J=1,IER1L=0,IERR=0 

HHH AK LEH FAK HK MK IH HL EE EI IIE E EE EEE 
INITIAL 


x CHECK FOR INITIAL RUN 


SS 


ЇЕ(1Ї1.6Т1.0) бОТО 200 

x THE DESIGN VARIABLES ARE THE FEEDBACK GAINS 

X INITIAL SYSTEM POLES TO BASE FEEDBACK GAINS UPON (COMPLEX FORM) 
Pl-CMPLX(-10.0,10.0) 
P2=CMPLX(-10.0,-10.0) 
P3=CMPLX(-300.0,0. ) 

* SET NON-ZERO COEFFICIENTS 
АА(1,2 )=50. 
АА( 2,1 )=-200. 
АА( 2,2 )=-200. 
АА(5,1)=-1. 
ВЗ(2)-200. 
CC(1;1)--1. 
H(2)=1. 
F(1,1)--59. 
Q(1,;1)-1. 
4(2,2)21. 
(1021. 

X ROUTINE TO FIND THE SYSTEM GAINS FROM KNOMN POLES (KUO SECTION 8-6) 
G(LI)zREAL(-P1-P2-P2)*AA(1,1)*4AA( 2,2) *AACS , 2) 
Gl 2 J=REAL(P1¥P2+P1*P3+P2*P3 )-AA(1,1 J¥AA(2Z,2)-... 
АА(2,2)%ДА(5;5)-АА(1,1)%ҰАА(5,5)%АА(2,5)%АА(5,2)“АА(1,2)% ... 
АА(2,1)%АА(1,5)%АА(5,1) 
6(5)-БКЕДІ(-Р1%ҰР2%ҰР5)%АА(1,1)%АА(2,2)%ҰАА(5,5)-АА(1,19)ж... 
АА(2,5)%ЖАА(5,2)%АА(1,;2)%ҰАА(2,5)ҰАА(5,1)-АА(1,2)%ҰАА(2,1)%... 
АА(5;5)Ұ%АА(1,5)%АА(2,1)%АА(5,2)-АА(1;5)ЖАА(35,1)%ДА(2,2) 
РА(1,1)-88(1) 
РА(1,2)-ВВ(2) 
РА( 1,5 )=ВВ(5) 
РА(2,1)=(-АА( 2,2 )-АА(5,5 ) )жВВ(1 )+АА(1,2 )жВВ(2 )+АА(1,35 )жВВ(5) 
РА(2,2 )=АА(2,1 )жВВ(1) +(-АА(5,5 )-АА(1,1) )ж8В(2 )+АА(2,5 )жВВ(5) 
РА(2,5 )=АА(5,1 )ж8В(1)+АА( 5,2 )хВВ(2 )+(-АА( 2,2 )-АА(1,1) )*ВВ(5) 
РА(5,1)-(АА(2,2)жАА(5,25)-АА(2,5)ҰАА4(2,2) 3ЖВВ(1)-(АА(1,2,)% ... 
АА(5,5)-АА(1,5)%АЛ(5,2))Ж8В(2)%( АА(1,2)%АА(2,5)-АА(1,5)% ... 


ДА( 2,2 ))Ж88( 3) 
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PACS3,2)2-CAAC2,1)XAA€ 5, 3) -AAC2, $5 )XAA (3 , 1) JXBBC1) 5(AAC 1, 1) 
AACS,3)-AACY,3 )J*AACS , 1) )XBBC2)-CAA C 1 1)*AAC 2,3 )-AA (15 3 )* 
^A 2,12)*8B8( 3) 
7AC5 523 j2! AAC2 o 17064365, 2 )- AACS 41 JXAA0U 2,2 ) JXBBC1)7CAAC1, 1)X 
АА( 5,2)-аАА(1,2)%АА( 5,1) 1%ҰБВ(2)%( АА(1,1)%АД(2,2)-АА(1,2209% ... 
АА(2,1) )%ВВ(5) 
WHRITE(CO,297)P1,92,P3 

297 FORMAT('O INITIAL POLES ДЕБ 7 -35( (75Е11.6;", ,Е11,6,” )) ')) 
CALL DGEFA(PA,3,3,IPVT,IER1) 
IF(IERI.NE.O)GOTO 300 
CALL 9GESL(PA,3,3,IPVT,G,0) 
МКІТЕ{ 6,299 )611),6(2),6(5) 

299 FORMAT(*O ІМІТІАІ бАТМБ5 6(11),6(2),6(3) з !,3Е15.9) 


WRITE(S,290) 


гэ 
“0 
о 


FORMAT '9 STATE CONTROL MATRIX (B AB) = ') 
WRITE(6,298) BB(1),AA(1,1)*BB(1)4+AA(1,2)*BB( 2) 
NRITE(6,298) BB(2),AA( 2,1 )*BB(1)+AA(2,2 )*BB( 2 ) 
298 FORMAT(2E15.4) 
WRITE( 6,291) 
291 ҒОКМАТ( "0 OUTPUT CONTROL MATRIX (CB CAB) = ') 
WRITE(6,298) CC( 1,1)*XBB( 1)4CC( 1,2) *BB( 2) ,CC( 1,1)*( AAC 1, 1) 
BB(1)+AA( 1,2 )*BB( 2) )+CC(1,2 )*( AA( 2,1 )¥BB( 1 )4AA( 2,2 )*BB( 2) ) 
WRITE( 6,292 ) 
292 FORMAT('G OBSERVABILITY MATRIX (Cx AxCx¥) = ') 
WRITE(6,298) СС(1,1),АА(1>1)%СС(1,1)%АА(2,1)%СС(1,2) 
WRITE(6,298) CC(1,2),AA(1,2)*CC(1,1)+AA(2,2)*CC(1,2) 
GOTO 302 
500 WRITE(6,301) ІЕВІ 
301 FORMAT('Q POLE SELECTION ERROR IER1 = ',17) 
CALL ENDJOB 
302 CONTINUE 
DV(1)=G(1) 
0У(2)-6(2) 
рУ(3 )-6(3) 


IF(INFO.EQ.0) GOTO 200 
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CALL ADSCINFO,ISTRAT ,IOPT ; IONED; IPRINT ; IGRAD;NDV ;NCON;,DV,VLB,;... 


VUB OBJ , CON, IDG ;NGT »1C,DF »WA,NRA »NCOLA »WK ;NRHK » LYK »NRIWK ) 
IWK{ 2 )=0 


ARITE( 6,220) 


і} 


120 FORMATE '9 SCALING NOT USED") 
* CALL TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 
200 CALL ADSCINFO,ISTRAT ,IOPT ,IONED,IPRINT ; IGRAD ;,NDV;NCON,; DV, VLB;... 
VUB ; OBJ ; CON, IDG, NGT ,IC,DF WA »NRA ,NCOLA , K >NRNK , IA^K ;NRIMK ) 
Gt 1)-2DV( 1) 
G( 2)zDVt 2) 
Gt 3 )zDV( 3) 
* CHECX FCR TYPE CF RUN IN ORDER TO SAVE INITIAL AND FINAL RESULTS 
ТЕСІПУМЕ 03500255 
ЗЗА?Н !ІМІТТА(/9ММ,89т1,бЕсТЕК618 ) ТІМЕ,Х1,Хх2,0) 
GOTO 245 
235 ІР(ІМҒО.МЕ.0)60Т0 240 
DELPRT=0.01 
DELPLT=0.01 
GRAPH (FINAL/RUN,RU=*,DE=TEK618 ) TIME,X1,X%2,U 
GOTO 245 
260 DELPRT=0.0 
DELPLT=0.0 
245 11-1 
ж INITIALIZE CONSTRAINT CHECK VALUES 
* LAST THREE VALUES OF THE CUTPUT 
ҮО501-01 
YOSC2-CZI 


YOSCS-CI 


* 


MAXIMUM VALUE OF THE OUTPUT 


СМАХ-0. 


* 


THE TIME CONSTANT OF THE OUTPUT 

TAU-D. 
% TIME WHEN COMPLETION TIME SPECIFICATION IS MET 
TC-0. 


ж INITIAL VALUE OF THE CONTROL VECTOR 
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UI=0. 
x MAXIMUM VALUE OF THE CONTROL 
ОМАХ-0. 
ж INITIAL VALUE OF YdE X2 STATE 
^-I-0. 
X MAXIMUM VALUE OF THE X2 STATE 
X2MAX=0. 
* VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
LSTOBJ=0. 
x SLOPE OF THE CBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
LSLOPE-0O. 
X VALUE OF TIME AT END OF LAST INTEGRATION INTERVAL 
LSTTIM=0. 
DERIVATIVE 
* PERFORM INTEGRATIONS TO FIND NEW STATES 
Xl-INTGRL(O. ,XX1) 
X2-INTGRL(O. ,XX2) 
AXl1-zAACl;l)XX1*AAC1,2 )*X2 *BB(C1)XU-*F(1;l)XWlI*F(1,2)xW? 


ХХ25-АА(2,1 В6Х14АА(2,2 X X2 *BBC2 )XU*F(V 2,1) XH1*F(CV 252 ) xMW2 


* 


FIND NEW CONTROL VALUE 


U=-GU1 )*KX1-G( 2 J*¥X2-G(3 J¥XINTGRL(O. »C) 


Ж 


FIND NEW OUTPUT VALUE 


C=CC(1,;1)%X1+CC(1,;2)%X2+H(1)%HW1+H( 2 )xXW2 


* 


FIND NEW OBJECTIVE FUNCTION VALUE USING QUADRATIC PERFORMANCE INDEX 
X AFTER TRANSLATION OF OUTPUT AND CONTROL VALUES 
X1T=X1-1. 
UT=U-1. 
OBJ=INTGRL(0. ,OOBJ ) 
OOBJ=(X1T¥Q( 1,1 )4+X2*Q( 1,2) )XX1T+(X1T*Q( 2,1 )4+K2*Q( 2,2) )*X2+ ,... 
UT*R( 1 )*UT 
X JE JEJE E JE JE JE JEJE JE JEJE JE JE JEJE JE JE JE JE JE JE JE FE 9€ JE JE JE JE JE JE JE C JE JE JE JE JE JE JE JE FE JE ME JE JE JE JE JE JE JE JE JE € € € € MEHE € € € E E E E EE 
DYNAMIC 
* DISTURBANCE FUNCTION 


ж W1=NORMAL(SEED,MEAN,SD ) 
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М121.%0.1Ж5ІМІНІЖТІМЕ) 
CONSTRAINT CHECKS AT EACH INTEGRATION 5ТЕР 
ESTABLISH MAX OUTPUT OSCILLATION CHECKS (YOSC) 
STEP FORWARD THE LAST THREE OUTPUT VALUES 
YOSC1=YOSC2 
YOSC2=YOSC3 
YOSC3=C 
CALCULATE THE DIFFERENCE OF THE ACTUAL OUTPUTS TO THE DESIRED OUTPUTS 
H1=H2 
H2-H3 
НЗ-АВ5(Ү05С3-СР) 
CALCULATE THE DIFFERENCE OF THE ACTUAL OUTPUT TO THE INITIAL OUTPUT 
DELC-ASGS(CI-C) 
CALCULATE THE OIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 
DELU-ABES(UI-U) 
FIND MAXIMUM CONTROL VALUE 
IF(DELU.GE .UMAX JUMAX=DELU 
CALCULATE THE DIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 
DELX2=ABS( X21 -X2 } 
FIND MAXIMUM X2 STATE VALUE 
IF(DELX2.GE.X2MAX )X2MAX-DELX2 
ESTABLISH TIME CONSTANT CHECK (TAU) 
IF(TAU.NE.O.)GOTO 20 
IF(DELC.GE. ABS(0.63212*( CI-CF ) ))TAUZ TIME 
20 IF(DELC.GE.CMAX)} | CMAX-DELC 
MAX AMPLITUDE CONSTRAINT (CON(1)) 
(*) IF ABOVE SPECIFICATION (-) IF BELOW SPECIFICATION 
IF(DELC.EQ.CMAX) CON( 1)=DELC-ABS(CI-CF )*( CSP/100*1. ) 
MAX COMPLETION TIME CONSTRAINT (CON(2)) 
SKIP CHECK IF COMPLETION TIME ALREADY MET 
IFUTC.NE.O.)GOTO 150 
FOR UNDER DAMPED CASE 
CHECK FOR A OSCILLATION REVERSAL, AND IF THE PEAK MAGNITUDE IS LESS 
THAN 17 OF THE DIFrERENCE BETWEEN THE INITIAL AND FINAL OUTPUT VALUES 


IF((H2.GE.H1).AND.(H2Z.GE.H3) 
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DABNBUSUHZODESABSPUGI-CEI*.01J)) TC-TIME 

CHECK FCR THE LAST THREE OUTPUT VALUES TO BE DECREASING AND ALL LESS 

THAN 14 OF THE CIFFERENCE BETNEEN THE INITIAL AND FINAL OUTPUT VALUES 
2212:122:1:12/1 21021821 7 .Н13,ЖМ 01.11. 

aBS((CI-CF }*.01))) TC=TIME 

FOR OVER DAMPED CASE 

CHECK FOR THE CUTPUT AT 5 TIMES THE TIME CONSTANT TO BE LESS 

THAN 1/4 OF THE DIFFERENCE BETWEEN THE INITIAL AND FINAL OUTPUT VALUES 
IF((C.GT.(CF-ABS(CI-CF)%.01)).AND.(TIME.GE.(5.%TAU )) 
.AND.(TAU.NE.0. ).AND.(C.EQ.CMAX)) TC=TIME 

ШЕ СВИРВФЕТТОМ CTZME NOT ME: SET CONSTRAINT POSITIVE 
CCON( 2 )=10. 

ScT CONSTRAINT VALUE NEGATIVE ЗУ VALUE FROM TIME SPECIFICATION 
F (I CONE.D. 1 CON(2)=TC-TSP 

MAX ALLCWASLE CONTROL VECTOR CONSTRAINT (CON(3)) 

(+) IF ABOVE SPECIFICATION. (-) IF BELOW SPECIFICATION 

150 IJF(DELU.EQ.UMAX) CON(3 J=DELU-USP 

UNSTABLE CONSTRAINT CHECK (CON(4)) 

(+1) IF OUTPUT IS TEN TIMES THE MAX ALLOWABLE (-1) OTHERWISE 
CON(4 )=-1. 
IF (ABS(C).GT.ABS(CI-CF )¥10.*(CSP/100.+1. ) JICON(4 )=1. 

TERMINATE RUN EARLY IF UNSTABLE TO PREVENT "BLOWING UP" 
IF(CON(4).GT.0.) CALL ENDRUN 

NEGATIVE RESPONSE CONSTRAINT CHECK (CON(5)) 

(41) IF OUTPUT INITIALLY GOES NEGATIVE (-1) OTHERWISE 
CON(5 )=-1. 
IF( (C ABS(CF-C).GT.ABS(CF-CI ))J. AND. ( TAU. EQ. O. ))CON( 5)721. 

TERMINATE RUN EARLY IF UNSTABLE TO PREVENT "BLOWING UP" 
IF(CON(5).GT.0.) CALL ENDRUN 

MAX ALLOWABLE X2 STATE CONSTRAINT (CON(6)) 

(+) IF ABOVE SPECIFICATION (-) IF BELOW SPECIFICATION 
IF(DELX2.EQ.X2MAX) CON( 6 J=DELX2-X2SP 

INTEGRATION COMPLETE CHECK 

50 IF((OBJ.LE.0O.).OR.(LSTOBJ.LE.O.))GOTO 40 


FOR OBJECTIVE FUNCTIONS THAT APPROACH A STEADY STATE VALUE CHECK FOR 
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X A OBJECTIVE FUNCTION SLOPE OF LESS THAN .14 OTHERWISE CHECK FOR THE 
x SECOND DERIVATIVE OF THE OBJECTIVE' FUNCTION TO ВЕ LESS THAN 17 
x ALONG WITH COMPLETION TIME BEING REACHED IN EITHER CASE 
SLOPE-(2BJ-LSTOBJ)/( TIME-LSTTIM) 
25.СРЕ“ЭБОӘЕ-ГСГОРЕ ӘСӘН 
ЇЕСГС5 28Ё-1Т..0019-0КЭРю5е ОВЕ  Е220081::2724 
.AND.(TC.NE.0.)) CALL ENDRUN 
x ADYANCE THE LAST FUNCTION, TIME AND SLOPE COUNTERS FOR NEXT INTERVAL 
40 LSTOBJ-0BJ 
ESUEIMSTIME 
LSLOPE=SLOPE 
TERMINAL 
IF (INFO.NE.O) GOTO 450 
30 4809 1-1,5 
GO 410 J=1,3 
410 V(I,J)-AA(I,J)-BB(I)XG(J) 
400 CONTINUE 
МКІТЕ(6,429936(1),6(2),6(5) 
499  FORMAT('O FINAL GAINS 6(1),G(2),G( 3) 7 ',S3bE15.4) 
X CALL TO FIND THE EIGNVALES HHICH ARE IN TURN THE FINAL SYSTEM POLES 
CALL RGINM,NN,V,WR,WI ,MATZ,Z,IV1,FV1L,IERR) 
IF(IERR.NE.O )HRITE(6,695)IERR 
495  FORMAT('O ERROR FROM EISPACK',;I3,'-FINAL POLES MAYBE MWRONG') 
WRITEL6,597)CLWROID),NIUILD),I-1,3) 
497  FORMAT('O FINAL POLES ARE : '5,3('(';E11.4,',',Ell.4,' J) %)) 
CALL ENDJOB 
450 CALL RERUN 
END 


STOP 
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APPENDIX E 
PROGRAM LISTINGS FOR POLE TO GAIN APPROXIMATIONS 


1, DETERMINANT METHOD 

The following listing would replace the initialization and initial segments of the 
DSL program listed in Appendix D, for using the determinant evaluation method to 
find the initial gains corresponding to desired initial system poles. The DERIVATIVE, 
DYNAMIC, and TERMINAL sections of the state variable program would remain 


unchanged. 


TITLE OC MOTOR INTEGRAL CONTROL SYSTEM-KUO EXAMPLE 8-7 

XIX XXX XXX5 XXXXXXXXXX XIX XXX XXX XXX X XXX XXXX X XXX X XXX XX XXXAX XX XXX XX XXX XX XXX 
* INITIALIZE THE ADS PARAMETERS 

FIXED NRA,NGT »NCOLA ,»NRWK »NRIWK , IGRAD ,NDV,NCON ,IDG,IC,IWK 

FIXED 75ТБАТ ,ТОРТ ,ТОМЕО ,ТРКЇМТ ,ТМЕО 

АКВАҮ DV(7),VLB(7),VUB(7),CON(5),IDG(5),IC(5),DF(7),NK(1000),INK(1000) 
D DIMENSION WA(40,40) 

PARAM NRA=40 ,NCOLA=40 »NRWK=1000 ,NRIWK=1000 ,NDV=3 ,NCON=1 ,IGRAD=0 ,INFO=0 
* ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

TABLE IDG(1-5)-3x-]1,2x2 

* LOWER BOUND FOR DESIGN VARIABLES 

TABLE VLB( 1-6 )=3*-1000. ,3x0. 

x UPPER BOUND FOR DESIGN VARIABLES 

TABLE VUB( 1-6 )=3*1000. ,3*0. 

* ADS OPTIMIZATION METHOD SELECTION 

PARAM ISTRAT=9, IOPT=5, IONED=7, IPRINT=3140 

HHH HHH IHF EHH HEHE HK хх HHH HHH HEH HHH HEHE HIE XOCX € CX ЭСЭС ЭЭС ЭС ЭСЭС ЭС ЭСЭС ЭС-ЭСЭС ЭСЭС ЭС ЭС-ЭС-ЭС ЭС HEH HIE 
x INITIALIZE THE MODEL PARAMETERS 

* SET CONTROL PARAMETERS 

РІХЕО 1,4, К,11,ТЕР1,ТЕКБВ,ММ,ММ,ТУ1,МАТ2 ,5ЕЕр 

АККАҮ BBUS ),HC2)RC2),GC B) )HRCS JD, WICS ),IV1C 3 ), FV1(3) 


D DIMENSION AA(3,2),CC( 2,2),Q(2,2),F( 2,2), V( 3,3),2(23,3) 
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x SET INITIAL PARAMETERS FOR POLE TO GAIN OETERMINATION 
TABLE G(1-3)=-.38,.6,-6. 

0 INTEGER NP,JOB,IPVT(3),INF 

0 СОМРТЕХэ16 Р(5),РА(3,3),РМС3,3),РОЕТ(2 ),РИК(5) 
EXCLUOE NP,JOB,IPVT ,INF »P»PA,PM,POET ,РҮК 

ж ZERO ALL MATRICES 

0 OATA AA,F,CC,Q/21*0.7 

TABLE BB(1-3)=3*0. ,R( 1-2 )=2*0. ,H( 1-2 )=2*0. 


* SET CONSTRAINT SPECIFICATIONS 


* 


CI=INITIAL OUTPUT VALUE} CF=FINAL EXPECTEO OUTPUT VALUE 


ж 


CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 


* 


TSP=MAXIMUI1 COMPLETION TIME 


+ 


USP =MAXIMUM VALUE OF THE CONTROL VECTOR 
PARAM CI=1. ,CF=0. ,CSP=40. ,TSP=1.5,USP=10. 
x SET INITIAL INPUT VALUES (NI-DISTURBANCE; W2=REFERENCE SETTING ) 
PARAM W1=0., W2=1., SEEO=113, MEAN=1., SO=0.1, WL=1.6 
* SET INITIAL RUN FLAG ANO EISPACK PARAMETERS 
PARAM II=0 ,MATZ=0 ,NM=3,NN=3,I=1,J=1,K=1,IER1=0,IERR=0 
XXXX5< XXX ЖЖ ЭЖ ЭЖЖ ЖЖК ХХ ХХ << ХХХ HEHEHE IE HEHEHE HEHEHE HEHEHE HEHEHE HE HEHEHE 
* INITIALIZE THE OSL PARAMETERS 
PRINT 08J,C,X1,X2,U 
SAVE (RUN) OBJ,X1,X2,U,C 
CONTROL OELPLT=0.005,0ELPRT=0.01 
CONTROL FINTIM=4. 
METHOD STIFF 
RELERR O8J=1.E-4 
ABSERR OBJ-1.E-4 
HHH HH HHH IEE HE HEE HEH IEE HEHEHE HE IEE HEHEHE E EXE QE HEE IEEE PEPE IE HEHE DE HE HEE HE IE EXE EE E EX € EXE EE EE EE HEHE HEH HEHE H 
INITIAL 
* CHECK FOR INITIAL RUN 
IFCII.GT.0O) GOTO 200 
* TRE DESIGN VARIABLES ARE THE FEEOBACK GAINS 
* INITIAL SYSTEM POLES TO BASE FEEDBACK GAINS UPON (COMPLEX FORM) 
Р(1)-СМРЇХ(С-10.0,10.0) 


P(2)=CMPLX(-10.0,-10.0) 


Эн 


Р(3)-0МРЇХ(-200.0,0.) 
* SET NON-ZERO COEFFICIENTS 

АА(1,2 )=50. 

аг 11-50200. 

ААС2,21)--2400. 

АА(5,125-1. 

88( 2 )=290. 

221(21,135-1, 

Н(21:1, 

БО151)5-50. 

Q(1,1)=1. 

Qt 2,2121. 

mU). 


* Sci THE JESIGN VARIABLE MATRIX 


DV( 2)26(2) 
DV(3 J=G(3) 
DV(4)-Q(1,1) 
DV(5)zQ( 2,2) 
DV(6 )-R(1) 
ж ROUTINE TO FIND THE SYSTEM GAINS FROM KNOHN POLES USING ADS 
NP =š 
JOB=10 
INFO=0 


ж USE ADS TO FIND G'S THAT MAKE SUM(DET(RELATION MATRIX)) с 0. 


660 CALL ADS(INFO,;ISTRAT,IOPT,IONED,;IPRINT,;IGRAD,;NDV,;NCON;DV,VLB,;... 


YUB ,08J, CON, IDG,NGT ; IC DF ; HA NR ;NCOLA HK ;NRAK , ІНК , МВІМК ) 
IF(INFO.EQ.O)GOTO 650 
DET-Q. 

* CALCULATE THE RELATION MATRIX: (P-A+BxG } 

ж 00 600 К-1,МР 
DO 610 I-1,NP 
DO 620 J=1,NP 

620 РА(1,415(-АА(1,33) 


610 РА(І,І)=Р(І)-АА(І,І) 
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DO 630 I=1,NP 
G(I)=DV(I) 
00 640 J=1,NP 
540 PM{J,I)=PA(J,I)+BB(J)*G( 1) 
550 CONTINUE 
* ALL 70 CONDITION THE RELATION MATRIX 
CALL ZGEFA(PM,NP NP , IPVT , INF ) 
* CALL TO FIND DERTERMINATE OF RELATION MATRIX 
CALL ZGEDI(PM,NP NP ,IPVT ,PDET ,PHK , JOB ) 
ж FIND OBJECTIVE FUNCTION-LOOKING FOR MINIMUM TO BE ZERO 
CON( K )=ABS( REAL( POET (1)*10.**PDET(2))) 
OBJ=OBJ+CON( K ) 
590 CONTINUE 
GOTO 660 
* OUTPUT INITIAL POLE AND GAINS INFORMATION 
380 ARITEt9,297)P(1),PC 2), PCS) 
297 ҒОЮМАТ( "0 DESIRED POLES WERE : ',3('(',E11.6,', ',E11.6,' J) ')) 
WRITE(6,299)G(1),G( 2),G( 3) 
299 FORMAT('O INITIAL GAINS G(1),G(2),G(3) - ',3E15.4) 
* CALL TO FIND THE EIGNVALES WHICH ARE IN TURN THE INITIAL SYSTEM POLES 
DO 680 I-1,3 
00 690 2-1,5. 
690 V(I,J)=AA(1I,J)-BB(I)xG(J) 
680 СОМТІММЕ 
CALL RG(NM;NN, V,HR,HI,MATZ,Z,IVl,FVl,IERR) 
IF( IERR.NE.0 )WRITE( 65485 JIERR 
485  FORMAT('D ERROR FROM EISPACK',IS,'-INITIAL POLES MAYBE WRONG’ ) 
МВЇТЕ(6,9873(М8(13,И1(13,151,3) 
487 ҒОНМАТ( “0 INITIAL POLES ARE™: ",5("(",Е11:4/ ЕНЕ УІ 
НЕЇТЕ(6,290) 
290 FORMAT('O STATE CONTROL MATRIX (B AB) = ') 
WRITE( 6,298) BB(1),AA(1,1)*BB(1)+AA(1,2)*BBI 2) 
WRITE( 6,298) BB(2),AA(2,1)*BB(1)+AA( 2,2 )*BB( 2) 
298 FORMAT(2E15.4) 


WRITE(6,291) 


29 


291 FORMAT( 'O OUTPUT CONTROL MATRIX (CB CAB) = ‘) 


NRITE(6,298) CC(1,1)*BB( 1)4CC( 1,2 )*BB( 2) ,CCC 1,1 )XC AA CI, 1) ... 


BB(1)+AA(1,2)*BB( 2) )+CC(1,2)%*(AA( 2,1 )*BB( 1 )+AA( 2,2 )*BB(2)) 
WRITE( 65292) 

292 FORMAT('O OBSERVABILITY MATRIX (Cx A%wC%) = ') 
ИВІТЕ(6,292) СС(1,1),АА(1,10%СС(1,1)%АА(2,10%СС(1,2) 
КОЇТЕ(6,298) СС(1,2),АА(1,2)ЖСС(1,1)4АА(2,2)5СС(1,2) 

x RESET ADS PARAMETERS FOR MODEL OPTIMIZATION 

CALL ENDJOB 

МСОМ-5 

Ї5ТВАТ-0 

102Т-4 

ТОМЕр-7 

IPRINT=1010 

І06‹1 )=0 

IDG( 2)=0 

IDG(3)=O 


* CALL TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 


200 CALL ADS(INFO,ISTRAT ,IOPT ,IONED,IPRINT,IGRAD,NDV,NCON,DY,VLB,... 


YUB ,OBJ ; CON; IDG,NGT ; IC; DF ;HA,NRA,NCOLA ;WK ,NRWK , IWK ;NRIMK ) 
GC 1)-0V(1) 
G(2)=DV( 2) 
G(3)=DV(3) 
Q(1,1)=DV(4) 
Q(2,2)=DV(5) 
R(1)=DV(6 ) 
% CHECK FOR TYPE OF RUN IN ORDER TO SAVE INITIAL AND FINAL RESULTS 
IF(II.NE.0)GOTO 235 
бВАРН (1МТТТА!/БОМ,КО-1,0Е-ТЕК618 ) ТЇМЕ,Х1,52,0 
| бОТО 295 
235  IF(INFO.NE.O)GOTO 240 
DELPRT-0.01 
DELPLT=0.01 
GRAPH (FINAL/RUN;RU-X*,DE-TEK618 ) TIME,X1,X2,U 


GOTO 245 
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240 DELPRT=0.0 


DELPLT=0.0 


245 11-1 


ж 


INITIALIZE CONSTRAINT CHECK VALUES 


* 


LAST THREE VALUES OF THE OUTPUT 


YOSC1=CI 


YOSC2=CI 


YOSC3=CI 


ж 


MAXIMUM VALUE OF THE OUTPUT 


CMAX=0. 


ж 


THE TIME CONSTANT OF THE OUTPUT 


TAU=0. 


* 


TIME WHEN COMPLETION TIME SPECIFICATION IS MET 


TC=0. 


* 


INITIAL VALUE OF THE CONTROL VECTOR 


UI=0. 


* 


MAXIMUM VALUE OF THE CONTROL 

UMAX=0. 
* VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
LSTOBJ-0. 


SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 


* 


LSLOPE=0. 


* 


VALUE OF TIME AT END OF LAST INTEGRATION INTERVAL 
LSTTIM=0. 


HHH IH HH HHI HEHE HEH EEF HHH IH HEH HEHEHE HEH HHH HHH HE HHH HHH HK HHH HHH € € € €9€9€ EXC €CXEC 3€ 39€ 3€ 3€ X 


2. EIGENVALUE METHOD 

The following program listing would be used if the the eigenvalue approximation 
method was desired for finding the gains corresponding to the desired initial system 
poles. The DERIVATIVE, DYNAMIC, and TERMINAL sections of the 8 
variable program in Appendix D remain unchanged, and would follow the listing 


below. 


TITLE DC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 


HEHE PEK HEIL € € 9€ 9€ 9€ € 9€ 3€ 3€ 9€ 3€ 9€ 3€ € 9€ € 3€ € € 9€ 3€ CC 9€ € 9€ € € E € € E € € € JEJE 9€ € € € € € € € 3€ 2€ € € € XC 3€ 9€ XX 3€ C € CX € 9€ € 
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* INITIALIZE THE ADS PARAMETERS 

FIXED NRA,NGT,NCOLA ,»NRWK »NRIWK ,IGRAD,NDV,NCON,IDG,IC,IWK 

ETXE ISTRAT,IOPT,IONED.IPRINT.INFO 

ARRAY ‘3V07),VLBO7*.VUBC 7 ),CON(5),206(5),IC(5),DF(7),NK( 1000) ,INK( 1000 ) 
0 QLIMENSION WA(S0,450 ) 

PARAM NRA=40 ,NCOLA=40 ,NRNK=1000 ,NRIKK=1000 ,NDV=3 »NCON=4 , IGRAD=0 , INFO=-2 
* ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

TABLE IDG(1-5 }=4%-1,2 

* LOWER BOUND FOR DESIGN VARIABLES 

TABLE VLB(1-6)-23*-1000.,23x0. 

x UPPER BOUNO FOR DESIGN VARIABLES 

TABLE VUOBIL1-65,23*1000.,3*0. 

ENSDS OPTITMIZATIOM METHOD SELECTION 

PARAM ISTRAT-1, I0PT-3, ICONED-5, IPRINT-z1100 

XC Je Me HEE ХЕ ХХХ EEE EEE SE IEE SE EE FE EE EEE DEE IE HE IEICE HEE FE IEE IEEE € 
* INITIALIZE THE MODEL PARAMETERS 

ж SET CONTROL PARAMETERS 

ЕТХЕр 1,4,К,11,1ЕВ1 ,ТЕББ,ММ,ММ,ТУ1,МАТ2,5ЕЕОр 


ARRAY 88(3),H(2}),R(07),G6(3),WR(3),WI(3),IV1(3),FV1(3) 


D ОІМЕ“І5ІОМ АА(5,25),СС(2,2),0(2,2)5Ғ(2,2),У(5,5),2(5,3) 
0 COMPLEX*16 P(3) 
EXCLUDE P 


X SET INITIAL PARAMETERS FOR POLE TO GAIN DETERMINATION 
TABLE G(1-3)-23*-1. 

* ZERO ALL MATRICES 

D DATA AA,F,CC,Q/21*0./ 

TABLE BB( 1-3)-23*0.,R( 1-2)z2X0. ,H( 1-2)z2x0. 

* SET CONSTRAINT SPECIFICATIONS 

% CI=INITIAL OUTPUT VALUE; CF=FINAL EXPECTED OUTPUT VALUE 
% CSP=MAXIMUM PERCENT OVERSHCOT FOR THE OUTPUT 

X TSP-MAXIMUM COMPLETION TIME 

x USP-MAXIMUM VALUE OF THE CONTROL VECTOR 

PARAM CI=1l.,CF=0. ,CSP=100. ,TSP=1.9,USP=1. 

x SET INITIAL INPUT VALUES (W1=DISTURSANCE, W2=REFERENCE SETTING) 


PARAM W1=0., W2=1., SEED=113, MEAN=1., SD=0.1, WL=1.6 
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x SET INITIAL RUN FLAG AND EISPACK PARAMETERS 
РАВАМ ІІ=0,МАТ2=0,М№М=3 ,М№=3,1=1,7=1,К=1,1ЕВ1=0,ІЕВА=0 
JE IE IE HEA HEHEHE 9€3€ 3€ 93€ 3€ 3€ 3€ 9€ 3€ 3€ 9€ 3€ 9€ 3€ 3€ 3€ € 9€ 3€ 3€ 9€ 3€ 3€ 9€ 3€ 9€ € 3€ 9€ 3€ 3€.€ €. 9€ 3€ 3€ 9€ 9€ 3€ 9€ 3€ EH HHI HHH HH 
* INITIALIZE THE DSL PARAMETERS 
MINT 0B4,C,41,X2,U 
SAVE (RUN) OBJ,X1,X2,U,C 
CONTROL DELPLT=0.005,CELPRT=0.01 
CONTROL FINTIM=4. 
METHOD STIFF 
RELERR OBJ=1.E-4 
ABSERR O8J=1.E-4 
FE 3€ 16 3€363€ 3€ 3€ 29€ 3€2€ 3€ 9€ € 2e 3€ 3€ 3€ € 3€ 3€ 9€ 3€ 3€ 9€ IE IE IE EE IIE HK HEHEHE IH IE EM IIE HEHEHE ЭС ЭСЭС ЭСЭС ЭС ЭС ЭСЭС ЭЄЭС ЭЄЭСЭСЭСЭС 
INITIAL 
+ CHECK FOR INITIAL RUN 
IF(12.GT, 6) 50100000 
ж THE DESIGN VARIABLES ARE THE FEEDBACK GAINS 
x INITIAL SYSTEM POLES TO BASE FEEDBACK GAINS UPON (COMPLEX FORM) 
Р(1)-СМРЇХ(-100.0,10.0) 
Р(2)-СМРЇХ(-100.0,-10.0) 
Р(5)-СМРЇХ(-300.0,0.) 
x SET NON-ZERO COEFFICIENTS 
AA(1,2)=50. 
AA(2,1)=-200. 
AAC2,2)2-200. 
ДА(3,1)2-1. 
BB(2)=200. ` 
БӘ =з: 
H(2)=1. 
Е(1,1)2-50. 
Gi 19537 
QU2,2)z1. 
R(1)=1. 
x SET THE DESIGN VARIABLE MATRIX 
DV(1)=G(1) 


DV( 2 )=G( 2 ) 
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0У(3)-5(3) 
DV(&)=Q(1,;1) 
DV(5)=Q(2,2 ) 


DV(6)=R(1) 


CALL ADS(INFO, ISTRAT,IOPT,IONED IPRINT;IGRAD,NDV,NCON,DV ,VLB, ... 


VUB ,OBJ ,CON ,IDG,NGT ,IC ,DF >WA NRA ,МСОПТПА , УС ,МБИК ,ТНК ,МЕТНК ) 
ІМК (2 )=0 
* ROUTINE TO FIND THE SYSTEM GAINS FROM KNOWN POLES USING ADS 
66D CALL ADS( INFO, ISTRAT ,IOPT ,IONED ,IPRINT > IGRAD NDV,NCON,DV,VLB,.. 
VUB ,OBJ ,CON,IDG,NGT ,IC,DF »NWA,NRA»NCOLA HK >NRWK » INK »NRIAK ) 
IF(INFO.EQ.D)GOTO 650 
X CALCULATE THE RELATION MATRIX: (A-B*G) 
DO 7DD I-1,3 
G(I)=DV(I) 
DO 71D J=1,3 
710 V(J,I)-AACJ,I)-BB(J)X*G(I) 
700 CONTINUE 
* CALL TO FIND THE RELATION MATRIX EIGENVALUES 
CALL RG(NM NN, V HR NI ,MATZ,Z,IV1,FV1 yIERR) 
IF(IERR.NE.D)MRITE(6,)791)IERR 
791 FORMAT('O ERROR FROM EISPACK',I2,'-FINAL POLES MAYBE WRONG' ) 
ТЕ(МАХ(МВ(1),ИВ(23МБ(33).6Т.0.26(1)5-6(1) 
CON( 1)-ABS(MAXWR(1) ;HRCÉ2 ) ;NRCS ) )-MAXCREALGCP(C12),... 
REAL(P(23),REAL(P(3)))) 
CON( 2 )JZABS(MIN(GNR(O1),;NR(G2);NRC2 ) ) -MIN(CREAL(CPC1)),... 
REAL(P(2));REALCPC3)))) 
CON( 3 J=ABS(MAXCWI(1),WNI(02),WI03) J-MAXCIMAG(P(1)),... 
IMAG(P(C2)),; IMAGCP(C2 )))) 
CON( & )ZABS(MIN(GWI(GLÓ),HIG2),NICGGS))-MINCGIMAG(P(C1)),... 
ІМАЄ(Р( 2 )),ІМАЄ(Р(5 )))) 
OBJ=CON( 1 )+CON( 2 }+CON( 3 )+CON( 4} 
GOTO 660 
* RESET ADS PARAMETERS FOR PERFOMRANCE INDEX OPTIMIZATION 
65D NCON=5 


IDG(4 }=2 
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296 


297 


299 


290 


298 


291 


292 


x CALL 


200 


WRITE(6,296 )P(1),P(2),P(3) 

FORMAT(C'O DESIRED POLES ARE : *,S5( (© 382279, °;* bina J) oe 
НВЇТЕС6,297ЭНВ(1 ),М1(1)»МКС2),НЇС2 ),МВКС3),Н1(5) 

ҒОНМАТ( "0 INITIAL POLES ARE : !,35(1!(!,Е11.9,",",Е11.9," 9) ")) 


WRITE( 6,299 )G(1),6(2),6(3 ) 


FORMAT('O INITIAL GAINS G(1),G(2),G(3) = ',5Е15.4) 
ARITE(6;,290) 
CALL ENDJOB 

FORMATI 'O STATE CONTROL MATRIX (B AB) = ') 
ARITE(6,298) BB(1),AA(1,1)*BB(1)+AA(1,2)*BB( 2) 
ARITE(6,298) BB(2),AA(2,1)XBB( 1) *4AAC 2, 2)*BB( 2) 

FORMAT(2E15.94) 


ARITE(6,291) 


FORMAT('O OUTPUT CONTROL MATRIX (CB CAB) = ') 
ARITE(C6,298) CC(C 1,1)*BBC1)*CCC1, 2)*BBC2),CCC 1, 1 )*CAAC1, 1)* 
BB(1)*AAC1,2)xXBB( 2) ) 9CC( 1,2 )*(C AAC2, 1 )*BBC 1)4A4A( 2 2 )*XBBC2)) 
НВЇТЕ(6,292) 

FORMAT('O OBSERVABILITY MATRIX (C* AxwC%x) = ') 
HRITE(C6,298) CC( 1,1),AAC 1; 1)*CC( 1,1) *AAC 2, 1)*XCC( 1,2) 
WRITE(6,298 ) CC(1,2),AA(1,2)XCC(1,1)+AA(2,2)XCC(1,2) 

TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 


CALL ADS(INFO,ISTRAT,IOPT,IONED,IPRINT »>IGRAD,NDV,NCON ,DV,VLB,... 
VUB ,0BJ ,CON, IDG »NGT ; IC, DF , HA, NRA ;NCOLA ,WK , NRWK , INK »NRIWK ) 
G(1)=DV(1) 

6121-0412) 
G(3)=DV(3) 
Q(1,1)=DV(%) 
Q(2,2)=DV(5) 


ВС1)-0У(6) 


x CHECK FOR TYPE OF RUN IN ORDER TO SAVE INITIAL AND FINAL RESULTS 


IF(II.NE.0)GOTO 2355 


GRAPH (INITIAL/RUN;RU-2,DE-TEK618 ) ТЇМЕ,Х1,Х2,0 


GOTO 245 


235 IFCINFO.NE.OIGOTO 240 


DELPRT=0.01 
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DELPLT=0.01 
GRAPH (FINAL/RUN,RU=*,DE=TEK618 ) ТЇМЕ,Х1,52,0 


GOTO 245 


* INITIALIZE CONSTRAINT CHECK VALUES 

x LAS: THREE VALUES OF THE OUTPUT 
ҮОЗС1-01 
YOSC2=CI 
Ү05Сс5-СІ 

* MAXIMUM VALUE OF THE GUTPUT 
CMAX=0. 

< NE TIME CONSTANT CF THE OUTPUT 
TAU=0. 

x TIME WHEN COMPLETION TIME SPECIFICATION IS MET 
TC-0. 

* INITIAL VALUE OF THE CONTROL VECTOR 
UI=0. 

x MAXIMUM VALUE OF THE CONTROL 
UMAX=0. 

* VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
LSTOBJ=0. 

% SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
LSLOPE=0. | 

* VALUE OF TIME AT END OF LAST INTEGRATION INTERVAL 
LSTTIM=0. 


хх X 3€ 3€ HEHE з ж HE IEE DE € HE € 3€ 9€ € 3€ € € € 3€ € € € € € 3€ € € 3€ € 3€ € 9€ 9€ € 3€ 9€ € 3€ € 9€ € 9€ € € Э6-эс ЭС-ЭС 36 эс эс эс ЭС эсээс 
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АРРЕМЮ В 
CASE 2 RESULTS FORE Ня 


j: RESULTS FOR STRATEGY 057 WITH SCALING APPLIED 


DSU/VS ..... RELEASE 1 ...-. MODIFICATION LEVEL 1 ..... 1985 


>>>> DSL SIMULATION INPUT DATA <<<< 
TITLE DC MDTDR INTEGRAL CONTROL SYSTEM-NO NDISE-KUD EXAMPLE 8-7 
FIXED NRA.NGT»,NCDLA , NRZIK , NRIWK , IGRAD,NDVNCCN, ТОС, ІС, ІМК 
FIXED ISTRAT.ICPT,IONED.IPRINT.INFO,IPVT 
PARAM NRA=40,NCDLA=460 -NRWK=1000 -NRIWK=21000-NDV=3 -NCON=S. [GRAD=0. INFO=0 
ТАВЬ 20601-5)53»0-2Ч2 
"МВ2Е СЛ.8(1-6)-:25-1000.,340. 
TABLZ VUB(1-$223*1000..3»*0, 
PARAM `STRAT=O, 1О0РТ-5, 10МЕ0=7,. IPRINT=3050 
SAVE 282,41,32,4 
ОВБАРН (0рЕ-5РЗ1МТ ) ТЇМЕ,0В.,.ХЇ,.Х2,.0 
PRINT OBJ.C.XIX2,U- е 
CONTROL FINTIM=4. 
METHDD STIFF 
RELERR 0BJz1.E-4 
ABSERR OBJs1.E-4 
FIXED I II IERI IERR:NN;NM,IVI MATZ 
TABLE ВВ(1-3):3мн0..8(1-2):3м0,,Н(1-2):340. 
РАВАМ С1:1,.,.СЕ-0, ,С5Р-40, ,Т5Р-1,5,05Р-:10, 
PARAM МІсз0., Мм241, 
РАХАМ 11:0,МАТ2:0,ММ-3,ММ-3,1:1.3:1,1ЕВ1-:0,.1ЕВВ-0 
END й 
чим К.ЕМОН = 4245. KPOINT = 1775, AVAILABLE SPACE LEFT IN CDMMON/CURVAL/ = 2470 DOUBLE-WORDS www 
чим S TIPE INTEGRATION METHDD USED www 
INITIAL POLES ARE : (-0.1000E*02, 0.1000E*02 J) (-0.1000E*02,-0.1000E*02 0) (-0.3000E+03, 0.0000E+00 J) 


INITIAL GAINS G(1).G(2)-GC03) = -0.3800E*00 0.6000Е%00 =0.6000Е +01 
STATE CONTROL MATRIX (В АВ) = 
0.0000Е +00 0.1000Е+05 
0.20002 +03 -0.4000Е +05 
OUTPUT CONTROL MATRIX (CB CAB) = 
0.0000E*00 -0,1000E*05 
DBSERVABILITY MATRIX (Сх Аниси) = 
-0,1000E*01 0.0000E*00 
0.0000Е%00 -0,5000E*02 
AAAAA DDDDDD 55555 
А A D D s 
A A D D s 
AAAAAAA D D sssss 
A A D D s 
A A D D S 
A A 000000 555555 
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CONTROL PARAMETERS 
ISTRAT s 0 


1GRAD 


ПІМЕ 
. 00000Е +00 
.00000E- 02 
.00000£-02 
. 00000£- 02 
00000£- 02 
00000E- 02 
.00000E- 02 
.00000E- 02 
.00000Е-02 
00000Е-02 
.00000Е-01 
11000 
12000 
.13000 
14000 
.15000 
.16000 
0.17000 
0.18000 
0.19000 
0.20000 
0.21000 
0.22000 
.23000 
.26000 
.25000 
.26000 
.27000 
.28000 
.29000 
- 30000 
- 31000 
» 32000 
-33000 
- 34000 
‚35000 
- 36000 
‚57000 
0.:8000 
0.39000 
0.40000 
0.41000 
0,42000 
0.45000 
0.44000 


о OO му = Ul £f t4 F9) о 


о O O O O = 


o 


о о фо о O O O O O O о о о о 


AUTOMATED 


- 0 


FORTRAN 


10РТ z 5 
NDV з 3 


мим DSL OUTPUT LISTING. 


ШЕЗЗ 1 5 М 


PROGRA 


SY Noten = S 1 


үс СВ-Ө 00321. 


IONED ° 
HCON з 
GROUP 1 


7 IPRINT = 3050 


DC MOTCR INTEGRAL CCNTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 


08. 
0. 00000Е +00 
1.94922Е-02 
3-78729Ё-02 
5.5798928-02 
7.00222Е-02 
8.36120Е-02 
9.69395Е-02 
0.10467 
0.11272 
0.11926 
0.12447 
0.12854 
0.13165 
0.13397 
0.15567 
0.13688 
0.13771 
0.13827 
0.13862 
0.13884 
0.13896 
0.13903 
0.13906 
0.13908 
0.13910 
0.13911 
0.13913 
02127315 
0.13918 
0.13922 
0.13925 
0.13929 
0.13953 
0.13936 
0.13940 
0.13943 
0.13945 
0.13947 
0.13949 
0512951 
0.13952 
0.13953 
0.13954 
0.115955 
0.13955 


C 
1.0000 
0.99486 
0.97439 
0.94014 
0.89509 
0.84193 
0.78319 
0.72092 
0.65689 
0.59259 
0.52921 
0.46775 
0.40896 
0.35343 
0.30157 
0.25365 
0.20983 
0.17015 
0.13459 
0.10304 
7.53384E- 02 
5.12988E-02 
3.06886E-02 
1.32590E-02 
-1.25109E-03 
-1.31081E- 02 
-2,257868-02 
-2499229Е-02 
-2.55969Е-02 
”-3,9226:2Е-02 
-4,16677Ё-02 
-4.09059E-02 
-4.31429Е-02 
~4.25561E-02 
-4.15057Е-02 
23.-332332Е-02 
-3.73683Е-02 
-3.69191Е-02 
23-22317Ё-02 
= 55 3s77E- 02 
=2067572Е=02 
22.39965E-02 
20.13020E-02 
-1.87104Е-02 
-1:52500Ё-02 


Х1 
0.00000E+00 
5.13756E-02 
e.,56087E-02 
95: 5955 1Е=02 
0.10491 
0.15807 
0.21681 

.27908 
0.34311 
0.40741 
0.47079 
0.5322 
0.59104 
0.64657 
0.69843 
0.76635 
0.79017 
0.82985 
0.86541 
0.89696 
0.92466 
0.94870 
0.9693] 
0.98674 

1.0013 

1.0131 

1.022 

1.0299 

1.0354 

1.0392 

1.0417 

1.0429 

1.0431 

1.0426 

1.0613 

1.0395 

1.0%74 

1.0349 

1.0323 

1.0295 

1.0268 

1.0240 

150213 

1.0187 

1.0163 


X2 
0.00000E*00 
2.52547E-02 
5,58460Е-02 
8.06192Е-02 
9.91254Е-02 
0.11263 
0.12165 
0.12683 
0.12878 
0.12805 
0-12513 
0.12047 
0.11448 
0.10751 
9 .98488E-02 
9.17725Е-02 
8.36994Е-02 
7,52150Е-02 
6.70709E-02 
5.9188 7Е- 02 
5.16667Е-02 
4.45697Е-02 
3.79544Е-02 
3.18513Е-02 
2.62778Е-02 
2.12392Е-02 
1,67296Е-02 
1.27346Е-02 
21055358 С-05 
6.20011Е-03 
3.60431Е-03 
1.41584Е-03 

-4.05144Е-04 
-1.88703E-03 
-3.06614Е-03 
-3.375927Е-03 
-4.64895Е-03 
-5.11580Е-05 
-5.40541Е-03 
-5.54461Е-05 
-5.558 3552053 
-5.46898Е- 03 
-5.29686Е-03 
~5.06008E-03 
-4,77465Ё-03 
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U 
0.00000Е +00 
4.67107E-02 
9.52865Е-02 
0.15103 
0,21202 
0.27627 
0.34196 
0.40765 
0.47215 
0.53451 
0.59399 
0.65003 
0.7022 
0.75039 
0.79433 
0.83402 
0.86952 
0.90094 
0.92847 
0.95229 
0.97267 
0.98984 

1.0041 
1.0157 
1.0249 
1.0320 
1.0372 
1.0408 
1.0430 
1.0440 
1.0441 
1.0433 
1.0619 
1.0400 
1.0577 
1.0552 
1.0524 
1.0296 
1.0268 
1.0240 
осте 
1.0186 
1.0161 
1.0138 
1.0116 


. START1NG RUN 


1 


чин 


0.45000 0.13956 -1.39415E-02 1201559 -4,454455Ё-03 1,0097 
0.46000 0.13956 -1:17992Ё-02 1.0118 -6.11187E-03 1.0079 
0.47000 0.13956 -9.85171:-05 1.0098 -3.75689E-03 1.0063 
SCALAR PROGRAM PARAMETERS 
REAL PARAMETERS 
1) ALAMDZ = 0.00000E+00 20) EXTRAP = 0.50000Е +01 
2) BETAMC = 0.00000€+00 ӨЛӘЗІЕШЕН з 0.10000Ё-01 
SJ CT з -0.50000Е-01 22) FOCHM = 0.10000Е-02 
4) CTL з -0.50000E-02 23) GMULTZ = 0.10000E+02 
S) CTLMIN = 0.10000E-02 24) РЗА12 = 0.95000Е%00 
6) CTMIN = Q.10000E-01 25) RMULT = 0.50900Е%01 
7) ОАВА5ШР з 0.10000Е-03 26) RMVLMZ = .20000Е%00 
B) DABDBJ = 0.129562-03 27) RP = 0.10000E+02 
9) DABDSM = 0.27913E-03 28) RPMAX з 0.10000Е%11 
10) DABSTR *  0.153586Е-05 29) RMULT з 0.20000Е-00 
11) DELALP = 0.50000E-02 30) ЕРРМІМ з 0.10000Е-09 
12) DELDBJ г 0.10000Е-02 31) RPPRIM * 0.10000E*03 
13) DELD8M = 0.10000E-01 32) SCFD з 0.10000Е+01 
14) DELSTR = 0.10000Е-02 33) SCLMIN = 0.10000Е-02 
15) 00831 = 0.10000Е +00 26) STOL = 0.10000Е-02 
16) DLOBJ2 = 0.10000Е+04 35) ТНЕТА2 = 0.10000E+00 
1777041 = 0.10000с<-01 36) XMULT = O.261805+01 
1B) DX2 = 0.20000E+00 57) 280 = 0.10000E-04 
19) EPSPEN з -0.50000Е-01 5B) PMLT = 0.10000Е +02 
INTEGER PARAMETERS 
1) ICNDIR = á &) ITRMDP = 6) JONED = 7 
27 15СА. - 1 5) ITRMST = 7) JTMAX = 20 
3) ITMAX = 40 
ARRAY STDRAGE REQUIREMENTS 
DIMENSIONED REQUIRED 
ARRAY 512Е 512Е 
WK 1000 299 
IWK 1000 196 
SCALING INFORMATION CALCULATED BY ADS 
OBJECTIVE FUNCTIDN SCALE ҒАСТОВ, SCFD = 0.46685Е +02 
DESIGN VARIASLE SCALE FACTDRS 
1 0.65528Е+01 0.16842Е +01 0.3747ВЕ +00 
CONSTRAINT SCALE FACTORS 
CONSTRAINT NUMBERS 
1 5 á 2 3 


SCALE FACTORS 
0.68265Е%01 0.10000Е%02 0.10000Е%02 0.10167Е%00 0.73710Е%01 
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]DPT = 5: MDDIFIED METHDD DF FEASI8LE DIRECTIONS 


INITIAL DESIGN 
ОВ) = 0.13956Е +00 


DECISIDN VARIABLES (X-VECTOR) 


1) -0.38000Е-00 0.60000Е%00 -0.60000Е%01 


LDHER 8DUNDS DN THE DECISION VARIABLES (VLB-VECTDR) 
1) -0.10000Е%06 -0.10000Е%06 -0.10000Е%04 
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UPPER BOUNOS ON THE DECISION VARIABLES (VUB-VECTOR) 
1) 0.10000Е%04 0.10000Е%04 0.10000Е%06 


CONSTRAINT VALUES (G-VECTOR) 
12 -39.26251E*01 -6.1031e0BE*00 -0.66014E+02 -9.10000Е%02 -0.10000Е%02 
-- EEGIM i ERATION NUMBER 
THERE ARE Q aCTIVE CONSTRAINTS AND G VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIOE CONSTRAINTS 


GRADIENT OF "HE OBJECTIVE FUNCTION (OF-VECTOR) 
1) 0.67675Е%00 0.1262BE+01 .22360Е +01 


SEARCH DIRECTION (S-VECTOR) 
1) -0.30266Е%00 -0.56476Е%00 -0.10000Е%01 


2ROPOSED ALPHA = :20000£*00 
SALCULATEO ALPHA s 9.22652Е%01 
OBJECTIVE = 9.57096Е-01 


SEC. STON J/ARIABLES (Х-МЕСТОВ) 
12 -9.55081820С6 -1.4948ВЕ%00  -0.14712Є»02 


SUNSTRAINT VALUES (2-МЕСТОН) 
1) -0.104573Е01 -0.1244ВЕ%00 -0.64694Е%02 -0.10000E+02 -0.10000Е%02 
~~ BEGIN ITERATION NUMBER 2 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRAOIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) 20,695517»00 0.55507Є6200 -0.19154Е-01 


SEARC" DIRECTION (S-VECTOR) 
1) 0.70557Е%00 -0.10000Е%01 -0.23626Е%00 


PROPOSED ALPHA = 0.17326Е%01 
CALCULATEO ALPHA = 0. 25660Е-01 
OBJECTIVE = 0.56591E-01 


OECISION VARIABLES (X-VECTOR) 
1) -0.52В05Е%00 -0.51012Е%00 -0.14729E+02 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.16636E*01 -0.12500E*00 -0.64511Е%02 -0.10000E*02 -0.10000E+02 
-- BEGIN ITERATION NUMBER 2 
ТНЕВЕ АВЕ 0 ACTIVE CONSTRAINTS АМО 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) -0.48796E*00 0.59645Е%00 -0.22670Е-01 


SEARCH DIRECTION (S-VECTOR) 
1) 0.75708E*00 -0.10000E*01 -0.11070Е%00 
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PROPOSEO ALPHA з 0.16455Е%01 
CALCULATEO ALPHA з 0, 69В95Е-02 
ОВОЄСТІУЄ з 0.56456с-01 


2:21210М JARIABLES (Х-УЕСТСЕ) 
or -0,221462440 -0.5140276»00  «0,14731:»02 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.16745Е%01 -0.12506Е%00 -0.64519Е%02 -0.10000Е%02 -0.10000E+02 
-- ЗЕСІМ ITERATION NUMSER 4 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) =0.4B500E+00 0.52880Е%00 -0.23835Е-01 


SEARCH ЗІРЕСТІОМ (5-МЕСТОЯ)|) 
1) 0.30312E+00 -0.10000E+01 -0.61035€-01 


SSOP OS=e0 ALPHA я 0,1635252- 01 
ZALCULATED ALPHA = 0.22228Е +00 
JBJECTIVE : 0.55864Е-01 


OECISION VARIABLES (X-VECTOR) 
12 -0.49983Е%00 -0.64625Е%00 -0.14767Е+02 


CONSTRAINT VALUES (G=VECTOR ) 
1) -0.175461E*01 -0.12181E+00 =-0.64550E*02 -0.10000E*02 -0.10000E+02 
-- 8EGIN ITERATION NUMBER 5 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATEO CONSTRAINTS ^ 
THERE ARE 0 ACTIVE SIOE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) -0.47829Е%00 -0.33528Е>00 -0.87296Е-01 


SEARCH OIRECTION (S-VECTOR) 
1) 0.100008 +01 0.70100Е%00 0.18252Е%00 


PROPCSEO ALPHA = 0.11464Е%00 
CALCULATEO ALPHA = 0.17350Е%00 
ОВЛЕСТІМЕ = 0.52701E-01 


OECISION VARIABLES (X-VECTOR) 
1) -0.47335E+00 -0.57406Е%00 -0.146826»02 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.20609€*01 -0.12899E*00 -0.65109Е%02 -0.10000Е%02 -0.10000Е%02 
-- BEGIN ITERATION NUM8ER 6 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 


THERE ARE 0 ACTIVE SIDE CONSTRAINTS 
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GRADIENT OF THE O8JECT1VE FUNCTION (DF-VECTOR) 
12 -0.24605Е%00 0.58072Е%00 0,44972Е-01 


SEARCH DIRECTION (S-VECTOR) 
13 0.100002%01 -0.34160Е%00 0.13126Е-01 


PROPCSEO ALPHA = 0.19789= +00 
CALCULATED ALPHA = 0.54227E+00 
OBJECTIVE = 0.50417E-01 


DECISION VARIABLES (X-VECTOR) 
13 -0.39060E*00 -0.68402E*00 -0.14663E*02 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.26250Е%01 -0.12582Е%00 -0.65743Е%02 -0.10000Е»02 
-- BEGIN ITERATION NUMBER 7 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE S1DE CONSTRAINTS 


GRADIENT OF THE OSJECTIVE FUNCTION (DF-VECTOR) 
1) -0.52639Е-01 -0.95469Е-01 0.10113Е%00 


SEARCH OIRECTION (S-VECTOR) | 
1) 0.10000E+01 0.69751E*00 -0.91949E+00 


PROPOSED ALPHA = 0.35788E+00 
CALCULATED ALPHA = 0.99532Е-01 
OSJECTIVE з 0.50195Е-01 


DECISION VARIABLES (Х-УЕСТОВ) 
1) -0.37561Е%00 -0.64280E*00 -0.14908Е%02 


CONSTRAINT VALUES (С-МЕСТОВ) 

1) -0.26885E*01 -0.1344á3E*00 -0.65935E*02 -0.10000E+02 
-- BEGIN ITERATION NUMBER 8 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE Q0 ACTIVE SIDE CONSTRAINTS 


GRAOIENT OF THE OSJECTIVE FUNCTION (OF-VECTOR) 
1) -0.19434Е-02 0.22563Е%00 0.14040Е%00 


SEARCH OIRECTION (S-VECTOR) 
1) 0.76057E*00 0.38162Е-01 -0.10000Е%01 


PROPOSED ALPHA = 0.32090Е +00 
CALCULATED ALPHA = 0.27000Е%03 
OBJECTIVE = 0.25264E-01 


OECISION VARIABLES (X-VECTOR) 
1) 0.30963Е +02 0256150Е01 -0.73533E+03 


CONSTRAINT VALUES (G-VECTOR) 
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-0.10000E+02 


-0.10000E+02 


-0.10000E+02 


-- ВЕОЇЧ ITERATION NUMBER 9 
THERE ARE 0 АСТІУЕ CONSTRAINTS ANO Q0 VIOLATEO CONSTRAINTS 
THERE ARE 0 ACTIVE SIOE CONSTRAINTS 


GRAOIENT OF THE OBJECTIVE FUNCTION (OF~VECTOR) 
1) 0.70317Е-02 0.13976E-01 0.59360E-02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.50312E+00 -0.10000Е201 -0.42831Е-00 


PROPOSEO ALPHA = 0.13505E+03 
CALCULATES ALPHA = 0. 00000Е •00 
OBJECTIVE = 0.352646E-01 


OECISION VARIABLES (X-VECTOR) 
1) 0.30963E*02 0.54750Е+01 =0.73533Е +03 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.27515E+01 -0.12889E+00 -0.66352E+02 -0.19000Е%02 -0.10000E+02 


+ 


-- BEGIN ITERATION NUMBER 10 
THERE ARE Q ACTIVE CONSTRAINTS AND 0 VIOLATEO CONSTRAINTS 
THERE ARE Q ACTIVE SIOE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) 0,/20217Е-02 0.13976E-01 0.59860Е-02 


SEARCH DIRECTION (S-VECTOR) 
l) =0.50312E+00 -0.10000Е%01 -0.42831E+00 


PROPOSEO ALPHA = 0.13500E+03 
CALCULATEO ALPHA = 0.00000E+90 
OBJECTIVE = 0.352646E-01 


OECISION VARIABLES (X-VECTOR) 
19 0.30963E*02 0.56750Е%01 -0.73533E+03 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.27515Е»01 -0.12889Е»00 -0.66352Е%02 -0.10000Е%02 -0.10000Е%02 
-- BEGIN [ITERATION NUMBER 11 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIOE CONSTRAINTS 


GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) 0.70317E-02 0.13976E-01 0.59B60E-02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.50312E+00 -0.10000Е01 -0.42831Е%00 


PROPOSED ALPHA =  0.10000Е-02 
CALCULATEO ALPHA = 0.57766Е%01 


OBJECTIVE = 0.33286E-01 
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DECISION VARIABLES (X-VECTOR) 
i) 0.30520E*02 0.20453E+01 -0.74193E+03 


CONSTRAINT VALUES (G-VECTOR) 

1) -2:276005-01 -0.12966Е%00 -0.66357Е%02 -0.10000Е%02 
-= ZZGIN ITERATION HUMBER с 
THERE ARE Q ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


СААОТЕМТ SF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) 0.72761Е-02 0.58022Е-01 0.56003E-02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.45297E+00 -0.10000€+01 -0.38651Е%00 


PROPOSED ALPHA =  0.28882Е»01 
CALCULATED “ALPHA = 90.50260Е%00 
OBJECTiVE з 2.23270E-01 


DECISION VARIABLES (X-VECTOR) 
m) J. 39485-92 2.176093Е01 -9.74265E*03 


2GNSTRAINT VALUES :0-УЁСТОВ) 

1) -9.27534E+91 -0.12857E+00 -0.66285Е+02 -0.10000Е+02 
-- 3EGIN ITERATION NUMBER 13 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) 0.70757Е-02 0.66616Е-02 0.54342Е-02 


SEARCH DIRECTION (S-VECTOR ) 
1) -0.10000E*01 -0.63056Е%00 -0.76801Е%00 


PROPOSED ALPHA = 0.31395Е%01 
CALCULATED ALPHA = 0.10690Е +0] 
OBJECTIVE s 0.33101E-OI 


DECISION VARIABLES (X-VECTOR) 
1) 0.30322Е%02 0.13466Е%01 -0.74464E+03 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.27598Е%01 -0.12961Е%00 -0.65688Е>02 -0.10000Е%02 


-- 8EGIN ITERATION NUMBER 14 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
19 0.72210Е-02 -0.10101Е-01 0.59956E-02 


SEARCH OIRECTION (S-VECTOR) 
1) -0.10000E*01 0.75990E-01 -0.78970E+00 
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-0.10000E*02 


-0.10000E+02 


-0.10000Е»02 


PROPOSED ALPHA = 0.7В5В1Е+00 
CALCULATED ALPHA = 0.58461E+02 
OBJECTIVE = 0.24306E-01 


DETISICN VARIABLES (Х-МЕСТОВ) 
2 1:21ё40025902 0.39863E*01 -0.B67B2E*03 


CONSTRAINT VALUES (G~VECTOR ) 

1) -0.27575Е%01 -0.14014Е+00 -0.66347E+02 -0.10000Е%02 
-- SEGIN ITERATION NUMBER 15 
THERE АВЕ g ACT1VE CONSTRAINTS AND OQ VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) 0.22096Е-02 0.10575Є6-01 0.14623Е-02 


SEARCH 01РЕСТ10М (5-ЧЕСТОЯ) 
1) -0.10000Е-01 -0.63193Е»00 -0.77107E+00 


PRCPOSED ALPHA = 0.29765Е-02 
CALCULATED ALPHA = 0. 00090Е +00 
Ј8ЈЕСТІУЕ = 0.24306Е-01 


DECISION VARIABLES (Х-УЄСТОВ ) 
1) 0.21400Е +02 0.39843Е+01 -0.86782Е +03 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.27575Е%01 -0.14014Е%00 -0.66347Е%02 -0.10000Е%02 
=- BEGIN ITERATION NUMBER 16 
THERE ARE 0 ACT1VE CONSTRAINTS AND G VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) 0.22096E-02 0.10575E-01 0.14623E- 02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.20895E*00 -0.10000E*01 -0.13828Е%00 


PROPOSED ALPHA = 0.29230€+02 
CALCULATED ALPHA = 0.00000E+*00 
OBJECTIVE = 9.26306E-01 


DECISION VARIABLES (X-VECTOR) 
1 0.21400E*02 0.39843E+01 -0.86782E+03 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.27575E+01 ~0.14014E+00 -0.66347E+02 -0.10000E+02 
-- BEGIN ITERATION NUMBER 17 
ТНЕВЕ АВЕ Q ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 


THERE ARE 0 ACTIVE SIDE CONSTRAINTS 
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-0.10000Е»02 


-0.10000Е%02 


-0.10000E*02 


GRAOIENT OF THE OS8JECTIVE FUNCTION (OF-VECTOR) 
I) 0.22096E-02 0.10575E- 0I 0.14623Е-02 


SEARCH DIRECTION (S-VECTOR) 
I) -0.20895E*00 -90.10000Е%01 -0.13828<%00 


PROPOSED ALPHA = 0.10000E-02 
CALCULATED ALPHA = 0.22054Е%01 
OBJECTIVE з 0.25914Є-01 


OECISION VARIASLES (X-VECTOR) 
I) 0.21330E+02 0.26749E+01 -0.86863E+03 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.27529Е01 -0.13901Е%00 -0.65241Е%02 -0.10000Е%402 
=- BEGIN ITERATION NUMBER 18 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CCNSTRAINTS 


GRAOIENT OF THE OS8JECTIVE FUNCTION (OF-VECTOR) 


1) 029995502 0.48979Е-02 0.11110£-02 
/ 


SEARCH OIRECTION (S-VECTOR) 
I) -0.38877E*00  -0.10000Є6»01 -0.17775Е%00 


PROPOSEO ALPHA =  0.11027E*0I 
CALCULATEO ALPHA =  0.95684Е-0Ї 
O8JECTIVE з 0.23913Е-01 


OECISION VARIASLES (X-VECTOR) 
1) 0.21324E*02 0.26180E*01 -0.86868Е%03 


CONSTRAINT VALUES (G-VECTOR) 

I) -0.27632E+01 -0.13983E+00 -0.65154Е%02 -0.10000Е%02 
-- BEGIN ITERATION NUMBER 19 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATEO CONSTRAINTS 
THERE ARE 0 ACTIVE SIOE CONSTRAINTS 


GRAQIENT OF THE O8JECTIVE FUNCTION (DF-VECTOR) 
1) 0-25915Е-02 0.25879Е-01 0.13133Е-02 


SEARCH OIRECTION (S-VECTOR) 
1) -0.10014Е%00 -0.10000Е%01 -0.50749Е-01 


PROPOSEO ALPHA = 0.11505E+01 
CALCULATEO ALPHA = 0. 00000Е +00 
OSJECTIVE = 0.25913Е-01 


DECISION VARIABLES (X-VECTOR) 
1) 0.21324E+02 .26180E+01 -0.86868E+03 


CONSTRAINT VALUES (G-VECTOR) 
12 -0.27632Е01 -0.13983Е%00 -0.65156Е%02 -0.10000Е%02 
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-0.10000E+02 


-0.10000E+02 


-0.10000E+02 


-- BEGIN ITERATION NUMBER 20 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT CF THE OBJECTIVE FUNCTION (LF-VECTOR) 
1) 0.25915Е-92 0,25879Е-01 0.15152Е-02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.10014Е%00 -0.1І0000Е%01 -0.50749Е-01 


PROPOSED ALPHA =s  0.47842Е-01 
CALCULATED ALPHA = 0.97728E-02 
OBJECTIVE * 0.23913E-OI 


DECISION VARIABLES (X-VECTOR) 
1) 0.21324Е +02 0.26122E*01 -0.86868Е +03 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.27635E*0I]I -0.13985E*00 -0.65146E*02 -0.10000Е%02 


FINAL OPTIMIZATION RESULTS 
NUMBER OF ITERATIONS = 20 
OBJECTIVE = 0,22917Е-0Т 


DECIS'IN VARIABLES (X-VECTOR) 
1) .21526Е»02 0.26122E+01 -0.86868E+03 


CONSTRAINT VALUES (G-VECTOR) 
I) -0.27635E*0I -0.I2985E*00 -0.65I66E*02 -0.10000E+02 


CONSTRAINT TOLERANCE, CT = -0.30000E-01 CTL = -0.50000E-02 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


TERMINATION CRITERIA 


-0,10000Е»02 


-0.10000E*02 


RELATIVE CONVERGENCE CRITERION WAS MET FOR 3 CONSECUTIVE ITERATIONS 


ABSOLUTE CONVERGENCE CRITERION WAS MET FOR 3 CONSECUTIVE ITERATIONS 


MAXIMUM К-Т RESIDUAL = 0.00000E+00 IS LESS THAN 0.10000Е-02 


cun Gun am ee сор ee ee — — с=ш= шә = = 


OBJECTIVE FUNCTION VALUE 0-259315Е-01 


DESIGN VARIABLES 


LOWER UPPER 
VARTAS8LE 8OUND VALUE BOUND 
-0.10000E+04 0.21324E+02 0.10000E+04 
-0.10000E+04 0;25122E*0] 0.10000Е +04 
3 -0.10000E+04 -0.86868E+03 0.10000Е»04 
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DESIGN CONSTRAINTS | 


11 -0.49068Е-00 -0.13278Е-01 -0.8838Е-01 -0.1000Е-01 -2.1000Е»01 
SUNCTION Z'/ALUATIONS : 173 
зө. cow 20ТРЫТ UL23TiNG. SROUP -. 
» STARTING RUN 175 www 
7 IME 08.) C XI X2 9 

9.30000Е-40 0.00000Е-00 1,3000 0.00000:200 0.00000Е%00 90.00000Е%00 
1.00000Е-02 1.197442-02 0.75301! 0.246099 0.69497 05842111 
2.20000Е-02 1.91402Е-02 0;47626 0.52376 0.42929 0.85409 
3.20090Е-02 2.19869Е-02 9.50355 0.62645 0.27368 0.90859 
G.Q00COE-02 2.31382Е-002 0.19342 0.80658 0.17440 0.94170 
5.00000Е-02  2.35021E-02 0.12329 0.87671 0.11115 0.96280 
6.90000Е-02 2.37926Е-02 7.85797Е-02 0.92142 7.08448Е-02 0.97630 
7.90000Е-02 2.38697Е-02 5.00829Е-02 0.94992 6.51536Е-02 0.98489 
8.00000Е-02 2.39009Е-02 3.19248Е-02 0.96808 2.87808Е-02 0.99037 
9.00000Ё-02 2.39136č-02 2.03445E-02 0.97966 1.85616Е-02 0.995386 

.20000Е-01 2.291882--02 1.296202-02 0.98704 1:215853Е-02- 0.29509 
9.11090 2,222092-02 2.25822Е-03 0.93917/6 7.64523с-03 0.99751 
FINAL GAINS 652).6(2),5(2) = Ө.с132Е02 0.2612E*01 -0.8687E*03 
FINAL POLES ARE : (-0.5387£*03, 9.2795Е%03 J} (-0.2287E*03,-0.2793E*03 J) (-0.4508E+02, 0.0000E+00 91 

SUMMARY: =NGING RUN NO. = 175, NO. SF RERUNS REQUESTED = 175 


чию SUN “ERMENATED ЗМ "250208" чч4 


RESULTS FOR STRATEGY 057 WITHOUT SCALING APPLIED 


2: 
DSL/VS ..... RELEASE 1 ..... MODIFICATION LEVEL 1 ..... 1985 
>>>> DSL SIMULATION INPUT DATA <<<< 

TITLE DC MOTOR INTEGRAL CONTROL SYSTEM=NO NOISE=KUO EXAMPLE 8-7 

FIXED NRA,NGT,NCOLA , NRHK , NR IWK , 1GRAD. NDV;NCON, IDG, IC. IMK 

FIXED ISTRAT.IOPT:1O0NED:; IPRINT; INFO, IPVT 

РАКАМ МКА=&0, МСО А240 , МАИК = 1 000, МАІМК = 1000, МОМ = 3, МСОМ=5, 1І6КАр=0, ІМЕО=-2 

TABLE I1DG(1-5)=3*0,2"2 

ТАВЬЕ М.8(1-6):38-1000.,3 0. 

TASLE VU8(1-6)=3"1000..3*0. 

PARAM 1STRAT=O, IOPT=5, IONED=7, IPRINT=1050 

ЗАУЕ 08.,Х1.Х2,0 

GRAPH (DE=SPRINT ) ТІМЕ,08Ј,Х1,х2, 0 

РВЇМТ О8),С.ХІ,Х2,Ц 

CONTROL FINTIM=4. 

METHOD STIFF 

RELERR O8J=1.E-4 

A8SERR O8J=1.E=-4 

ЕЇХЕП 1,.11.1ЕВ1,1ЕВВ.ММ.ММ,1У1,МАТ2 

TABLE 88(1-3)23*0., RCI-2)23*0., HC1-2)22 3*0. 

РАКАМ С1-1.,СР-0..С5Р-40.,Т5Р-1,5,05Р-10. 

PARAM М1:0., W2sI. 

PARAM I1=0.MATZ=0.,NM=3.NN=3,1=1,J=1.1ER1=0,I1ERR=0 

ENO 
ими KLENGH = 4245, КРОІМТ = 1775, AVAILA8LE SPACE LEFT IN COMMON/CURVAL/ = 2470 DOU8LE-WORDS www 
иик STIFF INTEGRATION METHOO USED «мн 


INITIAL POLES ARE : (-0.1000Е%02, 
INITIAL GAINS G(1).G(2),G(3) 
STATE CONTROL MATRIX (8 АБ) 

0.0000Е%00 0.1000Е%05 


0.1000E*02 3) 
-0.5800Е%00 


(-0.1000Е»02,-0.1000Е%02 J) 
0.6000E*00 -0.6000E*01 


(-0.3000E*03, 0.0000E*00 J) 
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0.2000Е%03 =0.4000Е +05 
OUTPUT CDNTRDL MATRIX (C8 CA8) = 


0. 0000Е +00 -0.1000Е +05 
О035ЕВУАВ11.1ТҮ МАТА1Х (04 АЯС) = 
=0.1000Е +01 0. 0000Е +00 

2.2000Е»00 -0.5900E*02 


AAAAA DDDDDD ssssss 
А А D D s 
А А D D s 
ААААААА D D SSSSS 
A A D D s 
A A р р s 
А А DDDDOD 555555 


FDRTRAN PRDGRAM 


ДОЦ ТОМУ Ч” ЕС DES si) GN SYNTH Ges Is 


VERSIDN 1.10 


CONTRDL PARAMETERS 
ISTRAT = 0 IDPT - 5 IDNED = ? 
IGRAD = 0 NDV з 5 NCON z 5 
SCALING NDT USED 
мия DSL OUTPUT LISTING»; GROUP 1 


IPRINT = 10S0 


» STARTING RUN | ммм 
DC MOTDR INTEGRAL CDNTROL SYSTEM-ND NOISE-KUD EXAMPLE 8-7 
TIME DBJ C X1 X2 U 


0.00000E*00 0.00000E+00 1.0000 0.00000E*00 0.00000Е%00 0.00000Е%00 
1.00000E-02 1.94922E-02 0.99486 5.13756E-03 2.52547E-02 9.67107Е-02 
2.00000E-02 .3.78729E-02 0.97439 2.56087E-02 5.58440Е-02  9.52865Є6- 02 
3.00000E-02 5.47989E-02 0.94014 5.98551Е-02 8.04192Е-02 0.151093 
4.00000Е-02 7.00222Е-02 0.89509 0.10491 9.91256E-02 0.21202 
5.00000E-02 8.34120E-02 0.84193 0.15807 0.11263 0.27627 
6.00000Е-02 9.49395E-02 0.78319 0.21681 0.12165 0.34196 
7.00000Е-02 0.10467 0.72092 0.27908 0.12683 0.40765 
8.00000Е-02 0.11272 0.65689 0.34311 0.12878 0.47215 
9.00000E-02 9711726 0559252 0.40741 0.12805 0.53451 
1.00000E-01 0.12447 0.52921 0.47079 0.12513 0.59399 
0.11000 0.12854 0.46775 0.5322 0.12047 0.65003 
0.12000 0.13155 0.40896 0.59104 0.11448 0.7022 
0.13000 0.13397 0.35343 0.64657 0.10751 0. 5057 
0.14000 0.13567 0.30157 0.69843 9.98688E-02 0.79433 
0.15000 0.13688 0.25365 0.74635 9.17725E-02 0.83402 
0.16000 ОРЕ ЛИИ 0.20983 0.79017 8.34995Е-02 0.86952 
0.17000 0.13827 0.17015 0.82985 7.52150E-02 0.90094 
0.18000 0.13862 0.13459 0.86541 6.70709E-02 0.92847 
0.19000 0.13884 0.10306 5.9188/E-02 0.95229 


0.89696 
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0.20000 0.13896 7.53384E-02 0.92466 5.16647Е-02 0.97267 
0.21000 0.13903 5.12988Е-02 0.94870 6.65697Е-02 0.98984 
0.22000 0.13906 3.06886E-02 0.96931 3.79566Е-02 1.0041 
0.23000 0.13908 1.32590Е-02 0.98674 3.13513Е-02 1.0157 
2.24000 0.13910 -1.25109Е-03 1.0013 2.62778E-02 1.0269 
0.25000 0.13911 -1.31081Е-22 1.0131 2.12392Е-02 1.0320 
0.26002 0.12913 -2.25784E-02 1.022 1 57296Е-02 1.0572 
0.27000 0.12915 -2.99229E-02 1.0299 1.27366Е-02 1.0608 
0.28000 0.13918 -3.53949Е-02 1.0356 .23358E-03 1.0430 
0.29000 0.13922 -3.92342Е-02 1.0392 6.20011E-03 1.0640 
0.30000 0.13925 -6.16677Е-02 1.0617 3.60631E-03 1.0641 
0.31000 0.13929 -4.29059Е-02 1.0629 1.61386E-03 1.0633 
0.32000 0.13923 -4.31428Е-02 1.0431 -6.05164Е-06 1.0619 
0.33200 0.13936 -4.25561Е-02 1.0626 -1.88703E-03 1.0420 
0.34000 0.13940 -6.13057E-02 1.0613 -3.06616Е-03 1.0377 
0.35000 0.13943 -3.95339Е-02 1.0395 -3.97597E-03 1.0352 

. 26000 0.13945 -3.73683Е-02 1.0374 -4.64895Е-05 1.0324 
0.37000 0.13947 -3.69191E-02 1.0349 -5.11580E-03 1.0296 
0.38000 0.13949 -3.22817E-02 1.0323 -5.40561Е-03 1.0268 
0.39000 0.13951 -2.95377Е-02 1.0295 -5.564661E-03 1.0240 
0.40000 0.13952 -2.675722-02 1.0268 -5.55832E-03 1.0212 
0.41000 0.13953 -2.399652-02 1.0260 -5.66898E-03 1.0186 
0.42000 0.13954 -2.13020E-02 1.0213 -5.29686E-03 1.0161 
0.43600 0.13955 -1.871066-02 1.0187 -5.060035-03 1.0138 
0.44000 0.13955 -1.62500Е-02 1.0163 -6.77665E-03 1.0116 
0.45000 0.13956 -1.39415Е-02 1.0139 -6.45455Е-03 1.0097 
0.46000 0.13956 -1.17992Е-02 1.0118 -6.11187E-03 1.0079 
0.47000 0.13956 -9.83171E-63 1.0098 -3.75689E-03 1.0063 


— — — — — «шә «шы «ше «шә — — — — — — чай чай чай чай чай чай чай "— чир чаб чию чаб чно чаб НЬ Оз чш сш чш x s чи? uu ЄНЬ ЄНЬ чи? чи? чиэ чи å cmm unb ub unD cb cub cmm aum cum uum cb MD — тар сы» қыр тыр ee зар «ыр — 


1ОРТ * 5; MODIFIED METHOD OF FEASISLE DIRECTIONS 


чий A өш» "до чий се с се си чий «шә - си» чай чий s s чир чий h s == сыз чыз бшш чи? чиЭ ӨНЬ чиэ чиэ чиэ чи? р чи вэ чи? ЯНЬ чиэ чи? сню ЧАО чао чаб ЕЮ чаб Gib umm s UND ч к=з ош эз чыз єз — ча? «шо = з шш» uu cum m 


-- INITLAL DESIGN 


08) = 0.13956Е +00 


DECISION VARIASLES (X-VECTOR) 


I) -0.38000E+00 0.60000E*00 -0.60000E+01 


LOWER BOUNDS ON THE DECISION VARIABLES (VL8-VECTOR) 
1) -0.10000Е%06 -0.10000Е%04 ~0.10000E+04 


UPPER S80UNDS ON THE DECISION VARIASLES (VU8-VECTOR) 
1) 0.10000Е +04 0.10000Е +04 0.10000E+04 


CONSTRAINT VALUES (G-VECTOR) 


1) -0.35685E*00 -0.10217Е%01 -0.89558Е%01 -0.10000Е%01 -0.10000Е%01 
-- BEGIN ITERATION NUMBER Ї 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE O8JECTIVE FUNCTION (DF-VECTOR ) 
1 0.94990Е-01 0.45556Е-01 4:-1219351Е-01 


SEARCH DIRECTION (5-МЕСТОВ) 


1) -0.10900Е%01 -0.47959Е%00 -0.18897Е>00 
PROPOSED ALPHA s 0.11608Е%00 
CALCULATED ALPHA = 0.22788E+00 


OBJECTIVE = 0.12996E+00 
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DECISION VARIABLES (X-VECTOR) 
1) -0.607В8Е%00 0.69071Е%00 -0.60431Е%01 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.19685Е%00 -0.72B29E*00 -0.8/7893Е%01 -0.10000Е%01 
-= BEGIN ITERATION NUMBER 2 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRAOLENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) -0.59658E-01 0.81346Е-01 0.11749Е-01 


SEARCH DIRECTION (S-VECTOR) 
1) -0.21230E*00 -0.10000Е%01 -0.22827Е%00 


PROPOSED ALPHA = 0.17198E+00 
CALCULATED ALPHA = 0.90049E+00 
OBJECTIVE = 0.10507E+00 


DECISION VARIABLES (X-VECTOR) 
1) -0.79995E*00 -0.60977E*00 -0.62686E*01 


CONSTRAINT VALUES (G=VECTOR ) 
1) -0.69459Е-02 -0.78792Е%00 -0.85804Е+01 -0.10000E+01 
=- BEGIN ITERATION NUMBER 3 
THERE ARE 1 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
CONSTRAINT NUMBERS 
Ї 


THERE ARE 0 ACT1VE S1DE CONSTRAINTS 


GRADIENT OF THE O8JECTIVE FUNCTION (OF-VECTOR) 
1) -0,43352Е-00 0.78925Е-01 -0.22729Е-02 


GRADIENT OF CONSTRAINT NUMBER l 
1) -0.24212Е%01 0.16891Е%00 -0.53977Е-01 


K-T PARAMETERS. BETA = 0.50816Е%00 MAX. RESIDUAL = 0.491B4E+00 


SEARCH DIRECTION (5-УЄСТОВ ) 
1) 0.10000E*01 -0.18206Е%00 0.52629Е-02 


PROPOSED ALPHA = 0.56418Е +00 
CALCULATED ALPHA = 0.15068Е +00 
ОВЈЕСТІМЕ = 0.80689Е-01 


DECISION VARIABLES (X-VECTOR) 
1) -0.64927E+00 -0.43721E*00 -0.6247BE+01 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.32129Е%00 -0.11902Е%01 -0.В9122Е%01 -0.10000Е%01 


-- ВЕСІМ ITERATION NUMBER 4 


THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 
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-0.10000Е%01 


-0.10000Е%01 


-0.10000E*01 


THERE ARE Q ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE O8JECTIVE FUNCTION (DF-VECTOR) 
1) 9.12485E-01 0.52862Е-01 0.67057Е-02 


SEARCH OLRECTICN (5-УЕСТОВ) 
17 -0.22518ЕҰ92  еї.100006-01 -0.12685с%00 


“ROPOSED ALPHA = 2.5455 9Е +00 
CALCULATED ALPHA = 0.52268E-02 
OBJECTIVE = 0.80397E-O1 


CECISION VARIABLES ЇХ-чЕСТОВ) 
17 =C.eS0SCE S00 -0.44241Е-00 -0.62485Е01 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.52116Е%00 -09.11920Е%01 -0.89119E*0I -0.10000Е%01 
-- BEGIN ITERATION NUMBER Э 
THERE ARE Q ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
шыр е АРЕ J ACTIVE SIDE CCNSTRAINTS 


ЗВАСЇЕМТ ОҒ “HE OSJECT:VE FUNCTION (DF-VECTOR) 
23 3-413068-21 0.52925Е-01 9.66478Е-02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.23055с%00 -0.10000Е%01  Á-0.I2623E*00 


PROPOSED ALPHA з 0.77955Е-01 
CALCULATED ALPHA з 0.20409E+00 
OBJECTIVE = 0.70223E-01 


DECISION VARIA8LES (X-VECTOR} 
1) -0.69756Е%00 -0.64652Е%00 -0.62763Е%?01 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.50233Е%00 -0.12122Е%01 -0.88751Е%01 -0.10000Е%01 
-- 8ЕО1М ITERATION NUM8ER 6 
ТНЕВЕ АКЕ 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT O* THE O8JECTIVE FUNCTICN (DF-VECTOR) 
ID 0. SII21IE-01 0.446775Е-01 0.554992-02 


SEARCH DIRECTION (S-VECTOR) 
1) 0.45047Е-01 -0.10000Е%01 --0.II225E*00 


PROPOSED ALPHA = 0.10466€£+00 
CALCULATED ALPHA = 0.10466€+00 
OBJECTIVE = 0.67007E-O1 


DECISION VARIAS8LES (X-VECTOR) 
1) -0.69286E*00 -0.75118Е%00 -0.62860E*0I 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.22042Е%00 -0.115341Е%01 -0.88657Е%01 -0.10000Е%01 
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-0.10000Е%01 


-0.10000E*0I 


-0.10000Е%01 


-- BEGIN ITERATION NUMBER 7 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE S1DE CONSTRAINTS 


SRADIENT СЕ THE OBJECTIVE FUNCTICN (DF-VECTOR) 
-4,64291М:-41 27-576 366-02 2,1882 3E-02 


SEARC.: DIRECTION (S-VECTOR ) 
1) 0.45896E+00 -0.10000Е +01 -0.12388E+00 


PROPOSED ALPHA = 0.15437E+00 
CALCULATED ALPHA 3 0.20074E-01 
OBJECTIVE = 0.66636E-01 


DECISION VARIASLES (X-VECTOR) 
1) =0.68363Е+00 -0.77125E+00 -0.62885E+01 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.34051E+00 -0.11226E+01 -0.88804E+01 -0.10000Е%01 
-= SEGIN ITERATION NUMBER 8 
THERE ARE Q ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE O8JECTIVE FUNCTION (DF-VECTOR) 
1) -0.36555Е-01 -0.67581Е-02 0.19196E-02 


SEARCH DIRECTION (S-VECTOR) 
1) 0.10000E*01 0.18688Е%00 -0.52513Е-01 


PROPOSED ALPHA a 0.62366Е-01 
CALCULATED ALPHA = .24117Е-01 
O8JECT1VE = 0.66174Е-01 


DECISION VARIABLES (X-VECTOR) 
1) -0.65951E+00 -0.76680Е%00 -0.62898Е%01 


CONSTRAINT VALUES (G-VECTOR ) 

1) -0.37897E*00 -0.11891E*01 -0.89462E*01 -0.10000E*01 
-- ВЕСІМ ITERATION NUMSER 9 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE S1DE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUHCTION (DF-VECTOR) 
1) 0.82166E-02 0.94493E-02 0.38751E-02 


SEARCH DIRECTION (S-VECTOR) 
1) -0.42778E+00 -0.10000Е%01 -0.47765Е%00 


PROPOSED ALPHA =  0.22035Е-01 
CALCULATED ALPHA =  0.67373Е-02 


OBJECTIVE з 0.66144Е-01 


128 


-0.10000Е%01 


-0.10000E*01 


ВЕС1510М VAR1A8LES (X-VECTOR ) 
1) -0.66239E*00 -0.77353E+00 -0.62930Е+01 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.37421Е%00 -0.11589Е%01 -0.89255Е%01 -0.10000Е%»01 
-- 8EGIN ITERATION NUMBER 10 
THERE ARE 0 ACTIVE CONSTRAINTS AND OQ VIOLATED CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) 0.78501E-03 0.27141E-02 0.34708Е-02 


SEARCH DIRECTION (S-VECTOR ) 
1) -0.22617Е%00 -0.78196Е%00 -0.10000Е%01 


PROPOSED ALPHA =  0.15427Е-01 
CALCULATED ALPHA =- 0.21057E-02 
О8ЈЕСТІУЕ = 0.66142Е-01 


DECISION VARIABLES (X-VECTOR ) 
1) -0.662872%00 -0.77518Е%00 -0.62951Е%01 


CONSTRAINT VALUES (G-VECTOR ) 

1) -0.27204Е%00 -0.11599Е%01 -0.89222Е%01 -0.10000Е%01 
-- BEGIN ITERATION NUMBER 11 
ТНЕВЕ АВЕ 0 ACTIVE CONSTRAINTS AND Q VIOLATED CONSTRAINTS 
THERE ARE Q ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 
1) -0.86474Е-02 -0.80495Е-01 0,33776Е-02 


SEARCH DIRECTION (S-VECTOR) 
1) 0.10743Е-01 0.10000Е%01 -0.41960Е-01 


PROPOSED ALPHA = 0.44215E-02 
CALCULATED ALPHA =  0.11576Е-01 
OBJECTIVE з-: 0.65590Е-01 


DECISION VARIAS8LES (X-VECTOR) 
1) -0.66275Е%00 -0.76260Е%00 -0.62956Е%01 


CONSTRAINT VALUES (G-VECTOR) 
1) -0.29095Е%00 -0.12925Е%01 -0.89288Е%01 -0.10000Е%01 


е 


г» 8EGIN ITERATION NUMBER 12 
ТНЕНЕ АВЕ 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
THERE ARE Q ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE O8JECTIVE FUNCTION (DF-VECTOR) 
1) 0.92074E-01 0.21756Ё-01 0.12965Е-01 


SEARCH DIRECTION (S-VECTOR) 
1) -0.99496Е%00 0.10000E*01 -0.19292Е%00 
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-0.10000ЕФ01 


-0.10000E+01 


-0.10000E*01 


PROPOSEO ALPHA = 0.68407Е-02 
CALCULATED ALPHA = 0.00000E+00 
OBJECTIVE 2 0.65590E-01 


DECISION VARIABLES (X-VECTOR) 
13 -0.66275E*00 -0.76360E*00 -0.62956E+01 


CONSTRAINT VALUES (G-VECTOR) 

1) -0.39095E*00 -0.12935E*01 -0.8938BE*01 -0.10000E+01 
-- BEGIN ITERATI1CN NUMBER 13 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 
THERE ARE 0 ACTIVE SIDE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
D 0.92074Е-01 0.21756Е-01 0.12965Е-01 


SEARCH DIRECTION (S-VECTCR) 
1) -0.16000E+01 -0.23629E+00 -0.16081Е%00 


PROPOSEO ALPHA = 0.57878E-02 
CALCULATEO ALPHA = 0.10867Е-03 
OBJECTIVE = 0.65589Е-01 ` 


DECISION VARIABLES (X-VECTOR) 
1) -0.66285Е%00 -0,763263E*00 -0.62956E+01 


CONSTRAINT VALUES {£G-VECTOR) 

1) -0.39071E+00 -0.12935E*01 -0.89385E*01 -0.10000Е%01 
=- BEGIN ITERATION NUMBER 15 
THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 
THERE ARE 0 ACTIVE SIOE CONSTRAINTS 


GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 
1) 0.92007E-01 0.21717E-01 0.129 70E-0I 


SEARCH DIRECTION (S-VECTOR) 
1) -0.,10000Е-01 -0.23616Е»00 -0.16089Е»00 


PROPOSEO ALPHA с 0.10000Е-02 
CALCULATED ALPHA = 0.72687Е- 04 
OBJECTIVE з 0.65589Е-01 


DECISION VARIABLES (X-VECTOR) 
1) -0.66293E*00 -0.76365E+00 -0.62956E+01 


CONSTRAINT VALUES (G-VECTOR) 
17 -0.39054E+00 -0.13935Е+01 -0.89383Е+01 -0.10000Е +01 


FINAL OPTIMIZATION RESULTS 


NUMBER OF ITERATIONS + 14 
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=0.10000Е +01 


=0.10000Е +01 


-0.10000E*01 


OBJECTIVE 3s 0.65589E-01 


DECISION VARIABLES (X-VECTCR) 


1) -0.66293Е-00 


CONSTRAINT JALUES СЕЧЕСТОН 2 
7 RESSE Sn 


93.230668.60 - 
CCNSTRAINT TCLZRANCE, CT 
THERE ARE 0 ACTIVE 
THERE ARE 0 ACTIVE 


TERMINATION CRITERIA 


-0.7583565Ё-00 


CONSTRAINTS AND 


-0.62958Е%01 


-02.89383Е-01 


= -0.20000Е-01 CTL = 


SIDE CONSTRAINTS 


-0.10000Е%01 


-0.10000€+01 


-0.50000Ё-02 


0 VIOLATED CONSTRAINTS 


RELATIVE CONVERGENCE CRITERION WAS MET FOR 


ABSOLUTE CONVERGENCE CRITERION WAS MET FOR 


3 CONSECUTIVE ITERATIONS 


3 CONSECUTIVE ITERATIONS 


— — = = 


ee = < тг. 


28ЈЕСТ:УЕ FUNCTION VALUE 0.65589Ё-01 
DESIGN VARIABLES 
LOWER UPPER 
VARIABLE BOUND VALUE BOUND 
1 -0.10000Е404 -0,66293Е»00 0.10000Е +04 
2 -0.10000€*06  -0.76265E*00 0.10000Е%04 
5 -0.10009Е%04 -0.62956Е%01 0.10000Е +04 


DESIGN CONSTRAINTS 


1) -0.5905Е>00 -0.1396Е>01 -0.8238Е»01 -0.1000Ё-01 -0.1000Е»01 
FUNCTION EVALUATIONS = 90 
wee DSL OUTPUT LISTING. GROUP 1 
TIME OBJ С X1 х2 U 

0.00000Е%00 0.00000Ё”00 1.0000 0.00000Е%00 0.00000Е%00 0.00000Е%00 
1.00000€-02 1.904651E-02 0.99076 9.23799Е-03 5.24768E-02 0.10900 
2.00000Е-02 2.50669Е-02 0.93704 6.29613=Е-02 0.16774 0.29362 
3.Q00000E-02 4.70222E-02 0.82258 0.17742 0.28635 0.51578 
4.00000£-02 5.51022E-02 0.65769 0.36251 0.36476 0.73190 
5.00000Е-02 6.007775-02 0.46817 0.55183 0.38167 0.90580 
6.00000Е-02 6.28622Е-02 0.28607 0.71593 0.35750 1.0164 
7.00000£-02 6.42180E-02 0.13799 0.86201 0.24934 1.0604 
8.00000E-02 6.47266Е-02 3.96980E-02 0.96050 0.14261 1.0494 
9.000090Е-02 6.683096Е-02 -5.97585Е-02 1.0060 4.28098E-02 1.0042 
1.00000Е-01 6.68637Е-02 -7.61740E-03 1.0076 -5.09017Е-02 0.94839 
0.11000 6.49237E-02 1.89452E-C2 0.98105 -09-3250528-02 0.90176 
0.12000 6.51313Е-02 S§.S8203E-G2 0.94418 -7.29285E-02 0.87684 
0.13000 6.55820Е-02 8.75761Е-02 0.91262 -9.06556tE-02 0.87737 
FINAL GAINS G(1),G(2),G(3) = =0.6629Е +00 -20.7636E*00 -0:52968 401 


FINAL POLES ARE : 


(-0.1245Е +02, 


0.5156E*02 J) 


(20.1265E*0257025156E*02 J) 


(-0.2257Е%02», 


SUMMARY: 


ENDING RUN NO. = 2, NO. OF RERUNS REQUESTED = 
wes RUN TERMINATED BY 'ENDJOS' 
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» STARTING RUN 92 www 


0.0000E*00 J) 


APPENDIX G 
ADSLPC IMSL PROGRAM CHANGES 


Changes were done to the following IMSL routines to make them compatible 
with MS-FORTRAN as follows: 


Routine Line Number Change 


DGEAR 4260 ЭРР5т-792987690601 498 
5980 oHDGEAR to DGEAR ' 
3990 6GHDGEAR to DGEAR ' 
CEO 22. OHLEOTIB to 17Е ШИ 
ЛЕКТІ) define CHARACTER*6 NAME 
550 define CHARACTER IEQ 
560 define CHARACTER*6 NAMSET, NAMUPK 
570 NAMSET - 'UERSET" 
580 МАМЕО з ' 
990 ЇЇ! се LLL 
620 NAMUPK = NAME 
640 IOUNIT = ó 
880 deleted 
: 890 IF(NAMUPK.NE.NANPISEPIGG IO 25 
1050 6A1 to 1A6 (TWICE) 
1070 6A] to IA6 (TWICE) 
1090 6Al to LAG (TWICE) 
ІТ) 6A1l to 1A6 (TWICE) 
1170 deleted 
1180 NAMEQ= NAMUPK 


ү 


Been 1020 
2740 
2770 


EQRH3F 400 
3700 


ЕЛШІБЕК - 1525575906Е-05 
I SEO ONE INT ” 
GHEIGR EAJ EIGRE 


RDELP = 1.52878906E-05 
BOREDE RO EORR3F 


Numerical conversions were calculated as follows: 
e Original IBM data statement: DATA SEPS /Z3C100000/ 
e Find exponent: ЗС - 40 = -3 


e The resulting hexadecimal number is: 0.1 x 10° 


3 


е Converting to decimal the number is: 16.0 x 1074 


* Converting to exponential form the value 1s: 1.525878906E-05 
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1. 


APPENDIX H 
ADSLPC PROGRAM LISTING 


ADSLPC MS-FOTRAN77 PROGRAM 


$nofloatcalls 


$large 


$storage:2 


E 


DC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 


(Q3 X x 3€ x € Xx * € E JE € «Xx 3x € x x ** 33x Xx X € 30 FE FE x * 9o 


C 


LISTING OF ALL VARIABLES 


AA 


BB 


CC 


CF 


CI 


CON 


CSP 


DELC 


DELT 


DELU 


DELX2 


DV 


БУТ? 


II 


IC 


IDG 


MATRIX OF STATE COEFFICIENTS 

VECTOR CF CONTROL COEFFICIENTS 

SCALAR OUTPUT VALUE 

MATRIX OF STATE OUTPUT COEFFICEINTS 

FINAL EXPECTED OUTPUT VALUE 

INITIAL OUTPUT VALUE 

VECTOR OF CONSTRAINT VALUES 

MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 

DIFFERENCE OF THE ACTUAL OUTPUT TO THE INITIAL OUTPUT 
INITIAL INTEGRATION STEP SIZE 

DIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 
DIFFERENCE OF THE ACTUAL X2 STATE TO THE INITIAL X2 STATE 
VECTOR OF DESIGN VARIABLE GREADIENTS 

VECTOR OF INPUT COEFFICIENTS 

EIGENVALUE WORK VECTOR 

VECTOR OF DESIGN VARIABLE GAIN VALUES 

VECTOR OF OUTPUT COEFFICEINTS FOR THE INPUTS 
LOOP COUNTER 

EVALUATION COUNTER 

VECTOR OF CONSTRAINT GRADIENTS 


VECTOR OF CONSTRAINT TYPES 
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IERR 
IGRAD 
INFO 
TONED 
TODT 
IPRINT 
IFVT 
ISTRAT 
EI 
ІНК 

4 

JOB 
LSLCPE 
LSTOBJ 
LSTTZM 
MATZ 
MEAN 
NCOLA 
NCON 
NDV 
NGT 

NM 

NN 

NRA 
МАТИК 
NRHK 


OBJ 


PA 


SD 
SEED 
STEP 


TAU 


EIGENVALUE cRROR INDICATOR 
SLALAR FOR INDICATING METHOD FOR OBTAINING GRADIENTS 


SCALAR INDICATING OPTIMIZATION STATUS 


SCALAR INDICATING TYPE OF CNE-DIMESIONAL SEACH 


INDICATING 


TIPE OFCOSPISEHIZATION 


SCALAR INDICATING TYPE OF PRINT CONTROL FOR ADS 


LINPAC INTEGER WORK VECTOR 
SCALAR INDICATING TYPE OF STRATEGY 
EIGENVALUE INTEGER WORK VECTOR 

INTEGER WORK VECTOR FOR ADS 

LOOP COUNTER 

STATUS INDICATOR OF OPTIMIZATION PROBLEM 
THE 


SEUES OF OBJECTIVE FUNCTION FOR THE LAST INTEGRATION 


VALUE JF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION 
VALUE OF TIME AT THE PREVIOUS INTEGRATION STEP 
EIGENVALUE SCALAR VALUE TO INDICATOR EIGENVALUE DETERMINATION 
DSL NORMAL DISTRIBUTION MEAN 

SCALAR COLUMN DIMENSION OF WA MATRIX 

NUMBER OF CONSTRAINTS 

NUMBER OF DESIGN VARIABLES 

SCALAR ADS CONTROL VARIABLE 

DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 
DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 

SCALAR ROW DIMENSION OF WA MATRIX 

SCALAR DIMENSION OF IWK VECTOR 

SCALAR DIMENSION OF WK VECTOR 

OBJECTIVE FUNCTION VALUE 

VECTOR OF DESIRED INITIAL COMPLEX POLE VALUES 

MATRIX OF INITIAL GAIN DETERMINATION COEFFICIENTS 

STATE WEIGHTING MATRIX 

CONTROL WEIGHTING VECTOR 

STANDARD DEVIATION OF DSL NORMAL FUNCTION 

RANDOM MUMBER SEED FOR DSL NORMAL FUNCTION 


EVALUATION INTERVAL FOR CONSTRAINTS 


SYSTEM TIME CONSTANT 
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С 


TC 


ТЇМЕМО 


TOL 


TSP 


У 


UI 


UMAX 


USP 


UT 


VLB 


VUB 


М1 


42 


WA 


AL 


WO 


Х1 


ХІТ 


хг 


X2I 


Х2МАХ 


X2SP 


Y 


Ү05С1 


YOSC2 


YOSC3 


TIME WHEN COMPLETION TIME SPECIFICATION IS MET 

TIMER VALUE AT EACH EVALUATION INTERVAL 

INITIAL TOLERANCE FOR DGEAR 

MAXIMUM COMPLETION TIME 

SCALAR FEEDBACK CONTROL VALUE 

INITIAL VALUE OF THE CONTROL VECTOR 

MAXIMUM CONTROL VALUE 

MAXIMUM VALUE OF THE CONTROL VECTOR 

TRANSLATED FEEDBACK CONTROL VALUE 

MATRIX FOR FINAL POLE DETERMINATION USING EIGENVALUE 
VECTOR OF DESIGN VARIABLE LOWER BOUNDS 

VECTOR OF DESIGN VARIABLE UPPER BOUNDS 

DISTURBANCE VALUE 

REFERENCE SETTING 

ADS WORK VECTOR 

CIRCULAR FREQUENCY OF SINUSOIDAL DISTRUBANCE 

COMPLEX VECTOR OF FINAL POLES 

MOTOR SHAFT SPEED STATE VALUE 

TRANSLATION OF X1 STATE 

MOTOR CURRENT STATE VALUE 

INITIAL VALUE OF THE X2 STATE 

MAXIMUM VALUE OF THE X2 STATE 

MAXIMUM SPECIFIED VALUE OF THE X2 STATE 

VECTOR OF INTEGRATION FUNCTION VALUES 

VALUE OF SYSTEM OUTPUT OSCILATION AT TWO PREVIOUS INTEGRATIONS 
VALUE OF SYSTEM OUTPUT OSCILATION AT PREVIOUS INTEGRATION STEP 


VALUE OF SYSTEM OUTPUT OSCILATION 


YPRIME VECTOR OF INTEGRATION FUNCTION DERIVATIVES 


2 


EIGENVALUE DUMMY MATRIX 


C 3€3€ €3€3€9€9€ HEHEHE 9€ € 9€ € 3€ € HEHE HEHE HEHE HEHE HEHEHE HEHEHE IE JE JE HEHEHE JE JE JE JE 9€ € 9€ € 9€ € 9€ JE JE JE JE JEE HE JE JE JE JE JE JE ME JE HEHE FEDE NE HEH 


С 


СОММОМ АА(3,3),88(3),СС(3,3),Н(2),В(2),4(2,2),Е(2»2),6(3), 


5 М1,М2,0 


INTEGER IDG(5),IC( 20), INK( 400) 


REAL DV( 3), VLBCZ ), VUBCS ), CONCB5), DFC2),HAC20,20 ), NK( 400) 
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REAL 1,5ТТЇМ,Ї5ТОВ.),1,51ОРЕ 
COMPLEX P(3),N0(3),2(3,3) 
INTEGER IGWK(4) 
REAL FV1(3),GNK(69),V(3,3),Y(4) 
EXTERNAL DUMMY ,DERIVA 
OPEN(6,FILE='OUTPUT.DAT' } 
OPEN( 7,FILE='CON: ' ) 
БОГ Re eoe INITIALIZATION SEGMENT HH KH HH HH HH ж 
C INITIALIZE THE ADS PARAMETERS 
NRA=20 
NCOLA=20 
NR¥IK=4600 
NRIWK=400 
NDV=3 
МО0М55 
IGRAD=0 
INFO--2 


C ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 


IDG(1)--1 
106(2)5--1 
106(3)--1 
106(9)--1 
106(5)-2 


C LONER BOUND FOR DESIGN VARIABLES 
VU8(1)=1000. 

VUB(2)=1000. 
VUB( 3 )-1000. . 

C UPPER BOUND FOR DESIGN VARIABLES 
VLB(1)--1000. 
VLB(2)--1000. 
VLB(3)--1000. 

C ADS OPTIMIZATION METHOD SELECTION 
ISTRAT-1 
IOPT-3 


IONED-3 
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IPRINT=1000 


(C3€3€ 39€X€ 9€ € 9€ 9€3€ 3€ 3€3€ 3€ € 9€ MEIE € 9X € 9€ 9€ 9€ 3€ 3€ € 29€ 9€ € 9€ 9€ 9€ € 9€ 9€ € 3€ € € € 9€ 3€ 3€ € € 9€ € E CX € 3€ 9COCÓ€ € 2€ 9€ € € 9€ € 3€ 


C INITIALIZE THE INTEGRATION PARAMETERS 


FINTIM-1.0 


STEP-0.1 


TOL=1.E-3 


С.ЭЄЗЄЭЄЭСЭСЭСЭСЭСЭСЭС ЭСЭС ЭСЭС ЭЄЭСЭС ЭС ЭСЭСЭСЭС ЭС-ЭСЭ6ЭСЭСЭ6 FE IEEE IE IE IE IE SEE DE EE IE HEE IE IE IE € € 9€ 9€ € 3€ € 9€ € € 9€ € 9€ 9€ X 33€ X 44 HE IEE HEE 


SET CONSTRAINT SPECIFICATIONS 

CI-INITIAL OUTPUT VALUE; CF=FINAL EXPECTED OUTPUT VALUE 
CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 
TSP=MAXIMUM COMPLETION TIME 


USP-MAXIMUSM VALUE OF THE CONTROL VECTOR 


CSP=100. 


TSP=1.5 


USP=10.0 


Х25Р-10.0 


С SET INITIAL INPUT VALUES (H1=DISTURBANCE, W2=REFERENCE SETTING) 


W1=0. 
W2-1. 
SEED=113. 
ME AN=1 
5р-0.1 


WL=157 


C SET INITIAL RUN FLAG AND EIGENVALUE PARAMETERS 


C 


МАТ2-0 
NM=3 
NN=3 
11-0 


жом ж ж ж жє ж ж INITIAL SEGMENT ж ж X X % % % з 


C SET ALL MATRIX AND VECTOR INPUTS 


P{1)=CMPLX(-10.0,10.0) 
P( 2 J=CMPLX(-10.0,-10.0) 


P( 3 )ZCMPLX(-300.0,0.) 
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АА(1,11-0.Е0 
АА(1,215:50.ЕО 
ААС1,31:0.Е0 


АСС,ү)5-200.20 


14 


ААС(2,2)5-200.:0 
АА( 2,5 )=0.Е0 
AACS,1)--1.E0 
АА(5,2)0.Е0 
АА( 5,35 )=0.Е0 
BB(1)-0.E0 


88(21)-200.ЕО 


2С11,5120.20 
CC(2,1)=0.E0 
CC(2,2)=0.E0 
CC(2,3)=0.E0 
CC(3,1)=0.E0 
СС(3,21-0.ЕО 
СС(3,31):-0.ЁО 
НІ) )20.Е0 
H(2)21.E0 
F(1,1)=-50.E0 
F(1,2)=0.E0 
F(2,1)-0.E0 
Ғ(2,2)0.Е0 
0(1,1)-1.Е0 
Q(1,2)=0.E0 
8(2,1)-0.ЕО 
Q(2,2)-1.E0 
RC1)-1.E0 
R(2)=0.E0O 
DV(1)=-1. 


DV(2)=-1. 
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DV(3)==1. 

DO 610 I=1,NCON 
610 СОМ(І )=0.Е0 

C POLE TO GAIN ROUTINE 

TALL ADS( INFO, ISTRAT ,ZOPT,IONED,IPRINT ,IGRAD ,NDV ,NCON, DV VLB 
$ VUB,OBJ,CON, IDG ,»NGT ,IC ,DF ; WA »NRA ,NCOLA ,WK ,NRAK , INK »NRINK ) 

C TURN SCALING OFF 

IWK(2)=0 
660 CALL ADS(INFO,ISTRAT,IOPT ,IONED ,IPRINT ,IGRAD NOV ,NCON,DV,VLB 
$ VUB ,0BJ ;CON, IDG,NGT , IC ;DF ;4A ,NRA ;NCOLA , HK NRHK , IHK ;NRIHK) 
IF(INFO.EQ.0)GOTO650 
DO 700 I=1,3 
G(I)=DV(I) 
9710 JS 152 
710 V(J,I)=zAA(J,I)-BB(J)xG(I) 
700 CONTINUE 

C FIND EIGENVALUES FOR PROPOSED GAINS 
CALL EIGRF(V,NM,NN,MATZ,WO,Z,NM,FV1,IERR) 
RP1=REAL(P(1)) 
RP2-REAL(P(2)) 
RP3=REAL(P(3)) 
RIP1=AIMAG(P(1)) 
RIP2-AIMAG(P(2)) 
RIPZ-AIMAG(P(Z)) 
WRISREAL(AO(CIL)) 
WR2-REAL(AO( 2) ) 
WR3=REAL(WO( 3 )) 
WI 1=AIMAG(WO(1)) 
WI2=AIMAG(WO(2)) 
WI3=AIMAG(WO(3)) 
IF(AMAX1(MR1,MR2,MR3).GT.0.E0)G(1)=-G(1) 
CON(1)=ABS(AMAX1(MR1,MR2,MR3)-AMAX1(RP1,RP2,RP3)) 
CONI 2 )=ABS( AMINI WR1,WR2,WR3 )-AMIN1(RP1,RP2,RP3)) 
CONI 3 )=ABSt AMAXL(WIl,HI2,WI3)-AMAXLIRIP1,RIP2,RIP3)) 


CON(4)=ABS(AMIN1(WIl,WI2,WI3)-AMIN1(RIP1,RIP2,RIP3)) 
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OBJ=CON(1)+CON( 2)+CON( 3 )+CON( & ) 


бОТО 660 


C PRINTOUT INITIAL POLES AND GAINS 


650 


296 


297 


277 


290 


298 


291 


292 


WRITE( 6,296 IP(1),P(2),P(3) 

FORMA 0 DeEstReEDerOLes:  55(°(°;5610.9,',°,E10.49,° J) ')) 
WRITE( 6,297 )WR1,WI1,WR2 »WI2,WR3 WIS 

FORMAT('O INITIAL POLES: ',3('(',E10.6,',',E10.4,' J) ')) 


WRITE(6,299)G6(1),GL2),GCSS ) 


FORMAT('O INITIAL GAINS G(1),G(2),G(3) г ",5Е11.4) 
ARITE(6,290) 
FORMAT('Q STATE CONTROL MATRIX (В АВ) = ') 


WRITE( 6,298 JBB(1),AA(1,1)*BB(1)+AA(1,2)*BB( 2) 
МКІТЕ(6,298)В8(12),АА(2>1)%В0(1)%АА(2,2)%ББ(2) 
ҒОВМАТ( 2Е15.4) 
АКІТЕ(6,291) 
FORMAT('0 OUTPUT CONTROL MATRIX (СВ САВ) = ') 
HRITE(6,)298)CC(1,1)%BB(1)+CC(1,)2)%BB(2),CC(1,1)%(AA(1,1)%BB(1)+ 
$ AA(1,2)*BB( 2) )4CC(1,2 XC AAC 2,1 )*¥BB(1)+AA(2,2)*BB( 2) ) 
МКІТЕ( 6,292) 
ҒОЯМАТ( "0 OBSERVABILITY MATRIX (C¥ AxwC%x) = ') 
МКІТЕ(6,298)0СС(1,1),АА(1,1)%СС(1,1)%АА(2,1)Ж%СС(1,2) 


МКІТЕ(6,298)0СС(1,2);АА(1,2)%СС(1,1)%АА(2,20%СС(1,2) 


C CONTROLLER OPTIMIZATION ROUTINE 


C RESET OPTIMIZATION PARAMETERS 


200 


IWK (2 )=1 

ІЅТКАТ =0 

ТОРТ-4 

IONED-7 

IPRINT=3030 
IDG(1)=0 

106(2)-0 

IDG(3)=0 

IDG(4)=2 

CALL AOSCINFO,ISTRAT,IOPT ,IONED,IPRINT ,IGRAD »NDV ,NCON,DV,VLB, 


$ VUB,OBJ,CON,IDG,NGT,IC ,DF »NA,NRA »NCOLA ,WK »NRAK , LAK »NRIWK ) 
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IF(INFO.EQ.O)HRITE(C6,2300) 
300 FORMAT(' TIME Kl хг 0", 
s OBJ E 
C RESET CONSTRAINT AND RUN VALUES 
GC 1)-20V(1) 
GC 2)-0V(C2) 
G( 3 )-D0V( 3) 
II-II*1 
JOB -0 
YOSC1=CI 
YOSC2-CI 
ҮО563-01 
СМАХ-0.ЕО 
TAU-O.EO 
TC=0.E0 
U=0.E0 
UI=0.E0 
UMAX-0.EO 
LSTOBJ=0.E0 
LSLOPE=0.E0 
І5ТТІМ-0.Е0 
С INITIALIZE THE OGEAR CONSTANTS FOR THE INTEGRATION TYPE DESIRED 
ТЇМЕ-0.0ЕО 
TIMEND=0.E0 
DELT=0.001 
NDE =4 
METH=1 
MITER=0 
INDEX=1 
DO 201 I=1,NDE 
201 Y(I)-O.EO 
C INCREMENT THE INTEGRATION INTERVAL 
20 TIMEND=TIMEND+STEP 
WRITE(7,455)II,TIME 


G55 FORMAT(' EVAL #',I3," TIME=',F5.3) 
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C CALL DIFFERENTIAL EQUATION SOLVER ROUTINE 


CALL DGEAR(NDE ;DERIVA DUMMY ; TIME ; DELT ; Y , TIMEND,TOL METH, 


$ MITER,INDE::,IGHWK ,GHK ; IER) 


C UPDATE 


CZ ATE AND CONTROL VALUES 

х= (1) 

X2=Y(2) 

OBJ=Y(%) 

CzCC(l;l)X*X1*CC(1,2)0*X2*H( l)*Wl*H(2)*MW2 
U-z-G(1)0*X1-G(2)*X2-G( 5)*Y( 3) 


жо ж ож ж ж ж ж ж DYNAMIC SEGMENT 


PRINT STATES AND OBJECTIVE - LAST RUN ONLY 


IF(INFO.EQ.0)MHRITECO6,310)TIMEND, X1, X2,U,OBJ,C 


310 FORMAT( 6612.4) 


c UPDATES 


CONSTRAINT CHECK VALUES 
YOSC1-YOSC2 

YOSC2=YOSC3 

YOSC3=C 

Н1-Н2 

H2=H3 

Н5-АВ5( Ү05С3-СЕ) 
DELC=ABS(CI-C ) 


DELU=ABS(UI-U) 


C EVALUATE CONSTRAINTS 


IF(DELU.GE.UMAX )UMAX=DELU 
IF(TAU.NE.0.E0)GOTO 15 


IF(DELC.GE.ABS(0.63212E0%(CI-CF ) ) JTAU=TIMEND 


15 IF(DELC.GE .CMAX JCMAX=DELC 


MH HM Эё ЭЕ 3 3 ях 


IF(DELC.EQ.CMAX)CON(1)ZDELC-ABS(CI-CF )*(CSP/100*1.E0) 


IF(TC.NE.O.EO)GOTO 150 


IF((H2.GE. HI). AND. HQ. GE. H3). AND. (H2. LE. ABSCCCI-CF )*.01E0))) 


$ TC-TIMEND 


IF(CHS3. LT. H2). AND. HA. LT. HL). AND. HL. LT. ABSCCCI-CF )*. 01E0))) 


$ TC-TIMEND 


IF(CC.GT. (CF-ABSCCI-CF)*.01E0)).AND. TIMEND.GE.( 5. EOXTAU)).AND. 


$ (TAU.NE.O.E0). AND. (C. EQ. CMAX ) )TCZ TIMEND 
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CON( 2)=10.E0 

IF (TC.NE.0.E0 JCON( 2 )=TC-TSP 
150 IF(OELU. EQ. UMAX )CON( Z ) EOELU-USP 

ССМ( 4 )=-1.Е0 
ІР(А455:С).0Т.АВЅСІ-СЕ )%10.Е0%(С5Р/100.Е0+1.Е0 ) )СОМ( Ф )=21.Е0 
IF(CON(4).GT.0.E0 )JOB=1 
CON(5)=-1.£0 
IF((ABS(CF-C).GT. ABSUCF-CI )). ANO. ( TAU. EQ. 0. EO) JCON( 5)21. EO 
IF(CON(5).GT.0.E0)JOB-1 

C EVALUATE FOR ENO OF RUN 
IF(TIMEND.LE.STEPXZ, )GOTO 40 
SLOPE -ABS( ( OBJ-LSTOBJ )/( TIMEND-LSTTIM)) 
О5 ОРЕ «4851 ! SLOPE-LSLOPE )/( TIMENO-LSTTIM)) 
IF(CCUSLOPE. LT. .001E0 ). OR. (OSLOPE . LE. .001E0 ) ). ANO. ( TC. NE. O.E0)) 

5 20851 
40 LSTOBJ-OBJ 

LSTTIM-TIMENO 
LSLOPE=SLOPE 

C CHECK FOR END OF RUN CONOITION 
IF (JOB.EQ.1)GOTO 30 

C CHECK FOR MAX TIME 
IF(TIMENO.GE.FINTIM)GOTO 30 
GOTO 20 

C жє жє ж x ж ж ж ж TERMINAL SEGMENT 9 жє ж ж ж ж Ж ж 

C CHECK FOR ENO OF JOB 

30 IF(INFO.NE.O)GOTO 200 
00 400 I=1,3 
00 410 J=1,3 
410 V(I,J)sAA(I,J)-BBCIDXG(J) 
400 CONTINUE 

C CALCULATE FINAL POLES 
CALL ELGRF(V,NM,NN,MATZ,WO,Z,NM,FV1,IERR) 
WRITE (65497 WO 1) ,HO( 2) ,HO(Z ) 
WRITE 65499 )G(1),G6(2),6(3) 


497 FORMAT('O FINAL POLES: '53('(',E10.4,',',E10.5,' J) ')) 


139 


499  FORMAT('O FINAL GAINS G(1)5,G(2),G( 3) = ',3E11.4%) 

STOP 
END 

С ж ж ж жж ж ж ж DERIVATIVE SEGMENT ж ж ж ж ж ж * x 
SUBROUTINE DERIVA(NDE,TIME,Y,YPRIME ) 
COMMON AA(3,3),BB(3),00(3,3),H(2),R02),Q02,2),F(2,2),613), 

581,Н2,0 

REAL Y(NDE ), YPRIME(NDE ) 
Х1-Ү(1) 
Х2-Ү(2) 

С DIFFERENTIAL EQUATIONS FOR CONTROLLER 
XIDOT=AAC 1,1 )5Х14АА(1,2 *X24+BB( 1 )¥XU+F (1,1 XW1+F(1,2 )*N2 
AeDOTzAA(2,1)*X41*AA( 2, 2) &X2*BB( 2 JXU*F( 2,1)XH1*F( 2,2)X*W2 
YSDOT-CC( 1,1)*X1*CC( 1, 2)*X2*H(C1 )*W14H( 2)*MW2 
U--G(1)*X1-G( 2)*X2-G( 3)*Y( 3) 
UTzU-1.EQ 
YZYXT-X1-1.EO0 
OBJDOTz(X1T*XQ(1,1)4X2*XQ( 1,2) J*X1T*( X1TXQ( 21)*X2XQ( 2,2) )*X2* 

$ (UTXR(1)xUT) 

YPRIME( 1 )=X1D0T 
YPRIME( 2 )=X2DOT 
YPRIME(3 }=Y3DOT 
YPRIME (4 )=OBJDOT 
RETURN 
END 

C DUMMY ROUTINE REQUIRED BY DGEAR 
SUBROUTINE DUMMY (N,X,Y,PD) 
INTEGER N 
REAL Y(N),PD(N;N),X 
RETURN 


END 
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2 ADSLPC BATCH FILE 


echo off 

rem: ASK 15 а program to query the user for a yes or no response 
rem: KEDIT 1s the line editor program of choice 
rem: PUTIME displays the current clock time 
rem: PUTIMEB displays the current clock time and sounds the system bell 
rem: BLANK.DATis an empty file 

rem: SETPRN sends escapecodes to the printer 
echo off 

path adsl util; adsl program: fortran; ; utiltys 
cd adsl program 

s tant 

ask Edit & Compile ADSLPC.FDR? 

1Ғ errorlevel 1 goto start2 

:besgin 

del ADSLPC.DBJ 

sstart3 

KEDIT ADSLPC.FDR 

putime 

FORI ADSLPC.ADSLPC,NUL ;, NUL 

putime 

РА52 

1f ex1st ADSLPC.DBJ goto 110Кк 

putlmeb 

ask Re-edit & Compile ADSLPC.FDR? 

if errorlevel 1 goto start2 

goto start3 

:start2 

ask Link and Run the ADSLPC program? 

1f errorlevel 1 goto linke 

:link 

putime 

LINK ADSLPC*ADSIA*ADSIB*ADS2A*ADS2B*ADSSA*ADS3B* ADS3C* IMSLI* IMSL2* IMSLS,MDDEL87,NUL, FORTRAN 8087.LIB* FORTRAN FDRTRAN.LIB 
putime 

goto runl2 

:11nk2 

ask Run the ADSLPC model? 

if errorlevel 1 goto print 

:run]2 

copy blank.dat output.dat 

putime 

MDCELB7 

putime 

pause 

BROWSE output.dat 

sprint 

ask Printout results (if yes--verify printer on, paper aligned): 
1f errorlevel 1 goto donel 

setern [27 IN(B} 

setprn [15] 

PRINT OUTPUT.DAT 

ASK PRINT ADSLPC.FDR? 

IF ERRDRLEVEL=1 GOTD NEXT2 

PRINT ADSLPC.FDR 

:NEXT2 

:donel 

ask Edit 8 Compile ADSLPC.FDR? 

1f errorlevel 1 goto cone 

goto begin 

:done 

echo on 
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